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Preface 


This volume contains the Proceedings of The East Asia Joint Symposium on 
Fields and Strings 2021, which was held in hybrid form at the Media Center 
of Osaka City University on November 22-27, 2021. About 160 physicists 
from all over East Asia attended or stayed online for this symposium and 
more than 50 researchers presented their results in the invited lectures, 
the short talks or the poster session. Quantum field theory and string 
theory in the context of several exciting developments were discussed, which 
include frontiers of supersymmetric gauge theory, anomalies and higher 
form symmetries and several issues on quantum gravity and black holes. 

We thank all of the speakers and the participants of this symposium for 
their stimulating lectures and intensive discussions. It is our sincerest hope 
that this volume will not only help to advance our field by presenting the 
latest developments in research, but also to serve to inspire new generations 
of physicists. 
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The fluctuation effect of the order parameter in a holographic superconduc- 
tor model was analyzed under holographic prescription for Schwinger-Keldysh 
closed time contour for non-equilibrium system in the paper. The time- 
dependent Ginzburg-Landau effective action, which governs the dynamics of 
the fluctuating order parameter near the critical point, is derived. The time- 
dependent Ginzburg-Landau action is computed up to quartic order of the fluc- 
tuating order parameter and first order in time derivative in a semi-analytical 
approach. 


Keywords: Black holes, Gauge-gravity correspondence, Holography and con- 
densed matter physics 


1. Introduction 


Critical exponents describe scaling behavior of observables near a critical 
point of continuous phase transitions. The characteristic properties of criti- 
cal exponents are symmetry, dimension, and properties of order parameters. 
It is believed to be independent of the details of interaction. As a result, it 
is important to examine the critical region. The results of weak and strong 
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coupling can be compared. Holography will make it possible to analyze 
dynamics (e.g. changing temperature) in the critical region: the critical 
exponent and the specific heat in theories at a quantum critical point. The 
gauge/gravity correspondence will give insights of cuprate superconductors, 
which appear in phase structure of superconductivity vs quantum criticality 
in interesting strongly correlated systems. 

Holographic superconductors have become a driving force for establish- 
ing AdS/CMT (condensed matter theory). Holographic s-wave supercon- 
ductors are referred to Abelian Higgs models %3, where s-wave, p-wave, and 
so on refer to the spatial component of the Cooper pair’s wave function. 
The phase transition of an s-wave superconductor is from a black hole with 
no hair (normal phase) to one with scalar hair. The gap of the AC con- 
ductivity at T = 0 represents the frequency w ~ 87,. This is consistent 
with the experimental results of cuprate superconductors. Because it is 
different from the BCS prediction w ~ 3.57%, the holographic result implies 
the strong coupling limit. The paper* constructed the general off-shell 
Ginzburg-Landau formulation (analytic expressions) of a holographic su- 
perconductor. An inhomogeneous order parameter was used to calculate 
the scaling coefficient and gradient term. Holographic entanglement en- 
tropy has become a nice probe of the holographic superconductor phase 
transition first realized by papers? ’. 

In this work, we apply the holographic Schwinger-Keldysh approach® 
to a holographic superconductor model’. The time-dependent physics of 
a fluctuating order parameter is investigated. This model is applicable 
to nonlinear problems. The main goal is to introduce time-dependence 
in the Landau-Ginzburg action and to consider a non-equilibrium QFT 
framework. 


2. The time-dependent Ginzburg-Landau effective action 


In this section, we derive the time-dependent Ginzburg-Landau effective 
action from a holographic superconductor model. We analyze it in the 
probe limit without considering the backreaction: the scalar QED in AdS 
Schwarzschild geometry !°. Consider an analytic holographic superconduc- 
tor, the lagrangian of which also becomes the scalar QED as follows: 


So = f Paya -3ft - (Dy Do- mv], 0 


where fmn = VuGy — Vrap and D, = Vy — iga,. U(1) gauge symmetry 
in the bulk corresponds to U(1) global symmetry on the AdS boundary. 
The formation of a charged scalar hair causes the spontaneous symmetry 
breaking of U(1). The incoming Eddington-Finkelstein coordinate is used 
to parameterize the AdS Schwarzschild metric z“ = (r,v, 2"). 


Im(r) 


Fig. 1. The radial coordinate is complexified in the gravity dual. The picture is taken 
from the paper !? 


To analyze both fluctuation and dissipation, the boundary theory must 
be placed on the Schwinger-Keldysh time contour!!. A holographic dual 
has the complexified radial coordinate®. It circles the black hole horizon 
from infinity 1 to infinity 2. See Fig. 1. 

The solution of the gauge field is in the leading order 


ay = (1 - $). (3) 


If the regular boundary condition at the black hole horizon is required, it is 
known that ODE can be solved analytically at the special point uo = 2rp. 
This unique value corresponds to the critical point at which the phase 
transition occurs. 

The scalar field has mass m = —4 which saturates the Breitenlohner- 
Freedman bound. The AdS boundary expansion of 7, w~* becomes 


logr , 
U(r — sz) = be -3 , 
r r 
‘ - logr 9, 
U(r + s) = bos ~z eae Ta (4) 


Because y, ~* are not independent, ps (®;) are not always complex con- 
jugate of wp, (®,). The two modes have conformal dimension two and are 
normalizable. The two modes are canonical conjugates of one another!. 


We fix ®,, identified as the expectation value of condensate, by adding a 
proper boundary terms at the AdS boundary. 

The EOMs are nonlinear partial differential equations (PDEs). It is not 
easy to obtain an analytic solution for PDEs. Therefore, we consider the 
following perturbation scheme to simplify the problem in the high temper- 
ature phase where temperature T is slightly above T,. We consider (1) 
off-equilibrium fluctuation of ® while keeping 6?, 6+ close to the critical 
point. (2) A slight deviation from the critical point fg is described by a 
chemical potential perturbation O(d), (3) We require the hydrodynamic 
limit: the system grows slowly with increase of time. The hydrodynamic 
limit can be represented as the low ono y expansion O(0,). Each a field 
is described by the triple expansion APOY YOM, gO), 

We compute the partially on-shell action, which leads to the boundary 
effective action Seff. UV divergences are removed by counter-term action 
added at the AdS boundaries. Due to the triple expansion of the bulk fields, 
the effective action Serr = f d‘xL is expanded as 


L = LABDA 4 LABDA 4 pOAO) 4 cA) a (5) 
Effective action in each order becomes 
LODO — 21 = 92g psy _ 2 org, (6) 
Qin 20 


LODD = ôu 2 (05 - B})(Hz — D1) — (034, ¥9,) 
T 


£5) So; Da + (G8; + 8:9,)], (7) 


log 2 
4T 


D* 3u Pa, 
(8) 

LOMO — — 9,000129006i( 6, P*)? + 0.004666880, P* (P* P, + &,6*) 
— 0.000263406i [(P} $)? + (©, ®*)*] — 0.00105363i®, $, 07 &* 
+ 0.0208333(®* P% 6? + &,0,.6*), (9) 


a rT r 


1 1 
LD0 — -301 81) BAD, + F(1 + 31) OFA Pa + 


where we have introduced difference and average combinations as the (r, a)- 
basis: 


1 
b,=0—%2, Br = 5(1 + 02). (10) 


Equation (6) is the noise term. The coefficient should be a positive imag- 
inary number because of unitarity. The chemical potential correction is 


represented by equation (7). This also includes corrections of 676,. Equa- 
tion (8) is the first order of hydrodynamic corrections. These terms con- 
tribute the symmetric and retarded two-point correlations, which obey the 
standard fluctuation-dissipation theorem up to dj at lowest order in w. 
Equation (9) includes general 8 potential terms. 

Effective action has Z-reflection symmetry ©, + Pz and P$ 4 P3 with 
sign flips of effective action. This implies that effective action has complex 
coefficients according to!*. 
by Z2 symmetry. 


Higher-point correlations will be constrained 


3. Comparison with weakly coupled results 


We compare gravity dual with weakly coupled results. The time-dependent 
Ginzburg-Landau action Sg, (T > T.)1' becomes 


Sou = av f dta [DE (LHO + PU (L71) DA + OE (L7* BE], (11) 


where Oc = (8A, PL)T are the classical part and quantum fluctuation of 
the scalar condensate (the order parameter). Matrix elements become 


7¢(3) 


T ON e — c 
(ERO = [a + D(V, + 2ieAf)? -ral - EP], (12) 
Ww = z iT 
(E>) = coth & [173P (w) - (E-)4(u)] = E, (13) 


where signs in the time-derivative term change in R(A) and TeL = 1/[8(T— 
T.)]. The symbol “K” means the K-gauge: the time-component of the 
external gauge potential is equal to zero. The equation (13) is obtained 
from Kubo-Martin-Schwinger symmetry. 

We have two fundamental parameters in the action: »—the density of 
states, and D—the diffusion constant. 

We have a similarity between the holographic model and weakly coupled 
theory. 


e The imaginary part of PP$ and eq. (6) are positive. This is a 
requirement coming from unitarity. In the gravity dual, the same 
term is corrected by eq. (7), which does not exist in the field theory 
side. 

e the coefficient of the OO}. (and O(/6%) term contains inverse 
relaxation time, which shows Te, ~ ep” ~ ee Actually, both 
models are in model A of the Hohenberg-Halperin classification for 
dynamic universality class. Static and dynamic critical exponents 


v = 1/2 and z = 2. The effective lagrangian eq. (7) in the gravity 
dual shows the critical exponent of the same dynamical universality 
class ĝu ~ u — uo ~ Te — T. Critical exponents are calculated for 
holographic s-wave and p-wave superconductors +1. 

e Only two quartic terms are present in eq. (11). These cover only 
arrr-terms. Eq. (9), on the other hand, has 8 general quartic 
potential terms. The gravity dual predicts more general results of 
holographic superconductors. 


4. Discussion 


We have outlined our obtained results'. We showed that holography pro- 
vided a possible description of strongly coupled superconductors. In par- 
ticular, time dependent Ginzburg-Landau effective action was derived from 
a holographic superconductor. Effective action is up to the quartic order 
in the fluctuating scalar condensate. This is an effective action for charge 
degrees of freedom such as currents and charge-condensate coupling. The 
holographic model was comparable with weak coupled BCS superconduc- 
tor’! and demonstrated the same dynamic universality class. Our work 
shows that the holographic Schwinger-Keldysh contour approach for non- 
equilibrium physics is relevant to nonlinear problems in the gravity dual. 

The analysis in this work is for spatially homogeneous case, in which 
the charge diffusion part is decoupled. It would be interesting to consider 
spatially non-homogeneous case. 

We are interested in how to incorporate Kibble-Zurek scaling. We 
must analyze the influence of the noise on the time dependent dynamics in 
quenching physics at that time. 
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The IIB matrix model has been proposed as a non-perturbative definition of super- 
string theory since 1996. We study a simplified model that describes the late time 
behavior of the IIB matrix model non-perturbatively using Monte Carlo methods, 
and we use the complex Langevin method to overcome the sign problem. We 
investigate a scenario where the space-time signature changes dynamically from 
Euclidean at early times to Lorentzian at late times.We discuss the possibility of 
the emergence of the (3+1)D expanding universe. 


Keywords: Type IIB matrix model; Emergent space-time; Complex Langevin sim- 
ulation 
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1. Introduction 


Superstring theory is the most promising candidate for a unified theory of all 
interactions, including quantum gravity. The theory is consistently defined in 
ten-dimensional space-time, leading to the compacting of the extra dimensions 
into small compact internal spaces. These scenarios have been investigated 
perturbatively on D-brane backgrounds and result in a vast number of vacua, 
leading to the so-called string landscape. It is, therefore, interesting to see 
what happens when one includes non-perturbative effects and whether these 
play an essential role in determining the true vacuum of the theory. The type 
IIB matrix model [1] has been proposed as a non-perturbative definition of 
superstring theory and provides a promising context to study such problems. 

The type IIB matrix model is formally obtained by the dimensional reduc- 
tion of ten-dimensional M = 1 Super Yang-Mills (SYM) to zero dimensions. 
The theory has maximal M = 2 supersymmetry (SUSY), where translations 
are realized by the shifts A, + A, +a,1, u = 0,...,9. The eigenvalues of 
the bosonic matrices A, can therefore be interpreted as coordinates of space- 
time. Thus, in this model, space-time appears dynamically from the degrees of 
freedom of matrices. In the Euclidean version of the model, the Spontaneous 
Symmetry Breaking (SSB) of the SO(10) rotational symmetry down to SO(3) 
occurs, which implies the emergence of a three-dimensional space [2-8]. 

By Monte Carlo simulation [9], it was found that a continuous time emerges 
dynamically, and a three-dimensional space expands. In Refs. [10, 11], it turned 
out that the expanding behavior of the space obeys the exponential law at early 
times and the power-law at late times. In Ref. [12], however, it was shown that 
SSB comes from singular configurations associated with the Pauli matrices, in 
which only two eigenvalues are large. This problem has been attributed to an 
approximation used to avoid the sign problem, which turned out later to be 
unjustifiable. 

In Refs. [13-15], the Complex Langevin Method (CLM) [16, 17] was used to 
overcome the sign problem without the above mentioned approximation. When 
one applies this method, one should apply the criterion for correct convergence 
of the CLM [18-24]. In Ref. [15], we found a new phase in which the structure 
of space is continuous by applying the CLM to the Lorentzian type IIB matrix 
model. See also Refs. [25-27] for other related works. 

In this work, we study the bosonic version of the type IIB matrix model 
by using the CLM. We show the equivalence between the Lorentzian and Eu- 
clidean models, which implies that the space-time in the Lorentzian model is 
Euclidean. To realize the possibility of the dynamical change of signature from 
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Euclidean to Lorentzian, we introduce a Lorentz-invariant mass term in the 
action that breaks the equivalence. We find some evidence that the signature 
of space-time changes from Euclidean at early times to Lorentzian at later 
times. 


2. The type IIB matrix model 
2.1. Definition 


The action of the type IIB matrix model is given as follows: S = Sp + Sf, 

Ss -5 Tr([A", A"[Ay, Av), Se= -z Tr (H(CI)[A,,,¥]) , (1) 
where A,, (u =0,...,9) and Y are N x N Hermitian matrices, and I’ and C 
are 10-dimensional gamma matrices and the charge conjugation matrix, respec- 
tively, which are obtained after the Weyl projection. The A, and Y transform 
as vectors and Majorana-Weyl spinors under SO(9,1) transformations. In this 
study, we omit S to reduce the computational cost. 

The partition function is given by Z = f dAe*>», Due to the phase factor 
eiS», the model is not well-defined as it is, and in this work, we define it by 
deforming the integration contour. When we rewrite the partition function as 
Zo. dAe~§ , the action of the Lorentzian model is given as 


cee -in|-2 Tr(Fy:)? + Tr(Fis)”| l (2) 


where g? = 1/N and Fy, = i{A,, Av]. According to Cauchy’s theorem, one 
can rotate the Lorentzian matrices A, to the Euclidean ones A, since the 
integration contour of A, can be deformed keeping the real part of Š positive. 
The relationship between A,, and Au is 


Ao = e 8 Ay ; A; = e's A; : (3) 
Then, the Euclidean action is given by 

s Ï ` ` 

Š= IN [2TA +Tr(Ē;)?] , (4) 


which is positive-definite. Here we have defined Fuy = i[A,,, Av]. 


2.2. Equivalence between the Euclidean and Lorentzian models 


By using Eq. (3), one can derive the relationship between the expectation 
values of Tr AZ and Tr A? in the two models: 


1 -3r 1 ~ 1 i 1 z 
a) =e a (504) 3 (=m?) = ez (522) » (5) 
(x */ NOS N i N E 
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Fig. 1. We plot the expectation values of wy Tr AZ (Left) and wir A? (Right). Those 
of the Euclidean and Lorentzian models are represented by the triangles and the squares, 
respectively. The angles between the Lorentzian and Euclidean models of (q#TrA2) and 


(TrA?) are —37/4 and 7/4, which agree with Eq. (5). 


where (-),;, and (-), denote the expectation values in the Lorentzian and 
Euclidean models, respectively. In Fig. 1 (Left), ($ TrAZ) is shown, and the 
angle between (+ TrA?)p and (4TrAĝ)i is —37/4. In Fig. 1 (Right), (4 TrA?) 
is shown, and the angle between (+ TrA?)p and (}TrA?), is 7/4. These angles 
are in agreement with Eq. (5). 

These results are consistent with the fact that the Lorentzian and the Eu- 
clidean models are equivalent. Expectation values in the Lorentzian model can 
be obtained by simply rotating the phase of those in the Euclidean model. In 
particular, (7TrA@), and (+ TrA?),, are complex and the emergent space-time 
should be interpreted as Euclidean. 


2.3. Lorentz-invariant mass term 
To realize real time and space, we introduce a Lorentz-invariant mass term in 
the action. For the Lorentzian model, the action is 


= -ŻN[-2 Tr(Fy:)? + Te(Fy)’] -— EN [Te Ao)” = Tr(Ai)"] (6) 


with y > 0. Using Eq. (3), we find that the action for the corresponding 
Euclidean model becomes 


c= ÍN [2 Tr(a)? + Tr(Fis)?| 4 Ny eft [Te(ão)? + Tr(A;)?] » 


where the real part of the mass term is negative. If y < 0, the real part of the 
mass term in the Euclidean model becomes positive, and then the matrices can 
be rotated from the Lorentzian to the Euclidean, which implies the equivalence 
between the two models. 
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The same mass term was used to study classical solutions of the Lorentzian 
type IIB matrix model [28]: 


[A’, [Av Au]] = yAu = 0 - (8) 


For y > 0, one can obtain classical solutions with smooth space and expanding 
behavior. Classical solutions with Hermitian A, make the time and space real. 
For y = 0, the classical solutions are given by simultaneously diagonalizable 
Apn, which do not necessarily have an expanding behavior. For y < 0, there do 
not exist classical solutions with expanding behavior. 


2.4. The time evolution 


As mentioned in Sec. 1, time does not exist a priori, and we define it as 
follows. We choose a basis in which Ag is diagonal and its eigenvalues are 
in the ascending order: Ap = diag(ai,a2,...,an), G1 < ag S<- < an. 
Then, we define @ as Ap = + 77, agi, and the time t, as 


p 
to = J |ūry — Ox . (9) 
k=l 


Here, we introduce the n x n matrices A; (t) as (Ai) 4 (¢) = (Ai)k4a k4p Which 
represent the space at the time t. 


3. Complex Langevin method 


The complex Langevin method (CLM) [16, 17] can be applied successfully to 
many systems with a complex action problem. One writes down stochastic 
differential equations for the complexified degrees of freedom, which can be 
used to compute expectation values under certain conditions. Consider a model 
given by the partition function Z = f dx w(x), where x € R” and w(x) isa 
complex-valued function. In the CLM, we complexify the variables x € R” —> 
z € C”, and solve the complex Langevin equation with the Langevin time ø: 
dz, 1 Ouw(z) 
Te ae Oe 
The first term of the right-hand side of Eq. (10) is the drift term, and the 
second one is the real Gaussian noise with the probability distribution 


P(n(o)) x e71 S do Xalino)? (11) 


+nk(o) - (10) 


To confirm that the CLM gives correct solutions, we use the criterion that the 
probability distribution of the drift term should be exponentially suppressed 
for large values [23]. 
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3.1. Application of the CLM to the type ITB matrix model 


To apply the CLM to the type IIB matrix model, we make a change of variables 
[13]: ag = 0, a = D e™ for 2 < i < N, where we introduce new real 
variables T,. In this way, the ordering of a; is automatically realized. Initially, 
Q; are real, and A; are Hermitian matrices. To apply the CLM, we complexify 
Tk and take A; to be SL(N,C) matrices. The complex Langevin equations are 
given by 

dT Sef 


do ~~ On AO. 


UADu ISen 
do O(Ai)tk 

where Seg is obtained from $' in Eq. (6) by adding a term associated with the 

gauge fixing and the Jacobian term associated with the change of variables. 


+ (ni)e(o) , (12) 


4. Results 


In the following, we introduce a parameter € in the mass term: 
= -4N|-2 Tr(Fy;)? + Tr(Fis)?| = 5 Ny |e! Tr(4o)? — e~" Tr(4:)?] (13) 


to shift coefficients of Tr(Ao)” and Tr(A;)? slightly from pure imaginary, and 
set € = 7/10. 


4.1. Expectation value of the time coordinate 


Fig. 2. Expectation values of the eigenvalues a; of Ag for N = 32, y = 3 are plotted. The 
solid line corresponds to the Euclidean model, where the complex phase of the expectation 
values (a;);, is exp(—i37/8). From this plot, 0s tends to become 0 at late times (at both 
ends of the distribution). 


3r 


When y = 0, Eq. (3) holds, and we expect that (a;);, = e7'® (Qi) _. This is 
true because of the equivalence between the Euclidean and Lorentzian models, 
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and time is regarded as the Euclidean one. We measure the time differences 
(Aai) = (aii) — (ai), x e**. The emergent time is real if 0, = 0. 

In Fig. 2, we plot the expectation values of the time coordinates (ai), 
on the complex plane for N = 32,7 = 3. The solid line corresponds to the 
Euclidean model, where the complex phase of (a;);, is exp(—i37/8). From the 
plot, 0, tends to become 0 at late times (at both ends of the distribution). 


4.2. Time evolution of space 


[ ewecetttttresssescccssce | 


oe | 02+ 4 


0 L 1 f 0 L L L L f f 1 
—08 —0.6 —04 —0.2 0 0.2 0.4 0.6 0.8 —0.8 -06 -04 —0.2 0 0.2 0.4 0.6 0.8 
t t 


Fig. 3. (Left) 0s(t) is plotted against t for y = 3. All values of s(t) are about 0.2 and below 
the s(t) = 7/8 line, which corresponds to the Euclidean space. (Right) |R?(t)| is plotted 
against t for y = 3. We can see that the space is expanding slightly with the time t. 


The time evolution of the extent of space is given by R(t) = 
(Air (Ai(t))”) = e7%(4)| R2(t)|. Since the matrices A; are complex, R?(t) 
is also complex. The time t is defined in Eq. (9). From Eq. (5), we obtain the 
Euclidean space when 0;(t) ~ 7/8, and the real space in the Lorentzian model 
when 6,(t) ~ 0. Therefore, the signature of space-time can change dynamically 
in this model. 

In Fig. 3 (Left) and (Right), 6,(t) and |R?(t)| are plotted against t for 
N = 32,y = 3, respectively. All values of 0,(¢) are about 0.2 and below the 
6;(t) = 7/8 line, which corresponds to the Euclidean space. We can see that 
the space is expanding slightly with the time t from the plot of |R?(t)]. 


5. Conclusions 


In this work, the CLM was applied to the bosonic type IIB matrix model 
in order to overcome the sign problem. We showed that the Lorentzian and 
Euclidean models are equivalent and that expectation values in the two models 
are related to each other by some complex phase rotation. The expectation 
values (5) in the Lorentzian model are complex, and space-time is Euclidean. 
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We introduced the model with a Lorentz-invariant mass term, which is a 
promising way to realize real time and expanding space. Then, the Euclidean 
and Lorentzian models are not equivalent anymore for y > 0. We found that 
the time, which is extracted from the expectation values of the eigenvalues 
of Ap in the Lorentzian model, may be real at late times although they are 
complex near the origin. We also studied the evolution of the extent of space 
with time. We have seen some tendency that the space becomes closer to real 
than the original model. 

To obtain a three-dimensional expanding space, we expect that supersym- 
metry will play an essential role. We are currently investigating its effect, 
which we will report in the near future. 
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Recently computational resource of quantum computers sounds growing well. 
In this article, we discuss how we can apply this development to numerically 
simulate quantum field theories. In contrast to the conventional approach by 
(Marlov chain) Monte Carlo method suffering from the infamous sign problem, 
we work in Hamilton formalism and adopt quantum algorithms which do not 
rely on Monte Carlo sampling. After brief discussion on how to put quantum 
field theories on quantum computers, we present our recent numerical results 
on the charge-g Schwinger model, where q is an electric charge of a Dirac 
fermion. We observe an exotic phenomena such as negative string tension 
behavior in potential between heavy charged particles which essentially come 
from presense of non-small 6-angle. 


Keywords: Quantum computation, Gauge theory, Schwinger model, Lattice 
gauge theory, Confinement 


1. Introduction 


It seems that resource of quantum computers has been recently growing 
well. In particular, readers would remember the news in 2019 that Google 
claimed to achieve quantum supremacy. While the news gave contraversy 
on whether or not it is really quantum supremacy, today quantum computer 
is one of hot science topics in public media. It is now possible for anyone to 
use (small-scale) quantum computer in the cloud for free. Although public 
news on quantum computers usually mention industrial applications such 
sequrity, new medicine and so on, in this article, we would like to consider 
how these developments can help us to understand physical systems. In 
particular, we discuss applications of quantum computation to numerical 
simulations of quantum field theories (QFT), and present some of our recent 
results in this context !~°. Since QFT is a common language in various fields 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
use, distribution and reproduction in any medium, even commercially, provided that the 
original work is properly cited. 
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such as high energy physics, nuclear physics, cosmology and condensed 
matters etc., it is expected to induce developments in various fields. 


2. Conventional numerical approach to quantum field 
theory 


For many purposes, one of the main motivations to use a quantum computer 
would be to perform fast numerical calculations. The purpose of this article, 
the application of quantum field theory to numerical simulations, is not 
an exception but in this case the motivation is more specific. To explain 
this, we discuss how numerical simulations of quantum field theoriesa have 
been usually performed. In the conventional approach, we use Lagrangian 
(path integral) formalism and the expectation value of the operator O for 
Euclidean case is written as 


e- Sla] 
(0(0) = F LO () 


where S[®] is the action and “integral domain” is over all possible values 
of the field ® at each point in the spacetime. Of course we cannot directly 
perform a numerical evaluation of the infinite dimensional integral and we 
need some regularization. The most standard way to do this is to cut 
the spacetime into a lattice with finite size (lattice regularization). To 
reproduce to the original theory, we need to take the continuous limit where 
the lattice spacing goes to zero*. Numerical evaluation of the regularized 
integral is usually done by an algorithm called the (Markov chain) Monte 
Carlo method, where we regard the Boltzmann weight (x e~*[®]) as the 
probability that realizes the field configuration®. Then we approximate the 
integral by using the average over the generated samples: 


(O(®)) ~ ——_ E oa). (2) 


(samples) 


jEsamples 


So far, the conventional approch has been successful to some extent. One of 
the significant successes is the derivation of the nuclear force from the first 
principle by lattice QCD (quantum chromodynamics)*. However, if the 
Boltzmann weights are not positive real numbers, the probability interpre- 
tation cannot be applied directly and we need some tricks”. In particular, 


“Tf the original spacetime has infinite volume, then we should also take a infinite volume 
limit. 

bThere are various efforts within the framework of the path integral formalism, which 
have limited success as sign problem becomes stronger (see e.g. the paper®). 
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when the integrand is highly oscillating, it is likely hard to do efficient sam- 
pling. This infamous sign problem often occurs physically, for example, in 
the presence of topological terms, chemical potentials and real time. All 
these situations are closely related to very important physical problems. 

What would happen if we work in the Hamiltonian formalism rather 
than the Lagrangian formalism? In the Hamiltonian formalism, the sign 
problem does not exist from the beginning since the problem to be solved 
is not a (path) integral. However, there is a reason why the Hamiltonian 
formalism has not been used much in numerical simulations of quantum 
field theory: since state space of QFT is typically infinite, we need to regu- 
larize it but the dimension of the state space after regularization typically 
increases exponentially with the number of “degrees of freedom”. In other 
words, naively, computers have to memorize exponentially large vectors 
corresponding to the states and multiply huge matrices corresponding to 
operators. This is the case for classical computer. What about quantum 
computers? For quantum computers, at least for some problems in quan- 
tum field theory, there are known algorithms with exponential improvement 
in computational complexity compared to classical computation. Here we 
consider what will be possible when the resources of quantum computers 
are expected to be sufficiently developed in the future. 


3. Quantum field theory as qubits 


Let us see how we can put quantum field theories on quantum computers. In 
gate-type quantum computer, the basic unit is quantum bit (qubit) which is 
a quantum system with two dimensional Hilbert space. In order to directly 
apply quantum algorithms to QFT, it is convenient to rewrite QFT as a 
spin system. Here we discuss how we can achieve this in terms of simple 
examples. 


3.1. Fermion field 


It is easiest for fermions because state space associated with fermions be- 
comes finite dimensional by simply cutting the spatial directions by a finite 
lattice. This is essentially because of the Pauli’s exclusion principle. Let 
us explicitly see this in a free Dirac fermion in 1 + 1 dimensions. We cut 
the one dimensional space by a lattice with N sites and lattice spacing a. If 


we choose staggered fermion® as a lattice fermion, the Hamiltonian of the 
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lattice theory is given by 


; N-2 N-1 
H = Ta Dixan an Lae +m b> (-1)"xh xn, (3) 
n=0 n=0 


where m is fermion mass and the lattice fermion operators (xn, X},) satisfy 


In this case the dimension of the state space is finite in dimension 2%. 

To rewrite the system in a form that we can directly apply quantum 
algorithms, we map it to a spin system. To do this, we need a spin operator 
that satisfies the canonical anti-commutation relation (4). While such a spin 
operator is not unique, the most traditional one is the so-called Jordan- 
Wigner transformation’: 


n-1 
Xn — 1Yy ; 
Xn = a II (—iZ;), (5) 
j=l 


where (Xn, Yn, Zn) is the Pauli matrix (01, 02,03) located at site n. Then 
we find 


122 Hi 


This is the same as the Hamiltonian of the XY-model with (space depen- 
dent) longitudinal magnetic field. 


3.2. Scalar field 


Next let us consider scalar field. The large difference from the fermionic case 
is that the state space is still infinite dimensional even on lattice and fur- 
ther regularization is required. This is because lattice scalar field theory is 
technically equivalent to multi-particle quantum mechanics with appropri- 
ate interactions which have infinite dimensional Hlibert space. To explain 
how to put lattice scalar field theory on quantum computers, let us consider 
a single particle quantum mechanics for simplicity: 


1, w? 2 
H(z, p) = 5p a + V(2), (7) 


where V(x) is a potential and 


(x, p] =i. (8) 
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In principle there are many ways to truncate the Hilbert space but here 
we use the most naive truncation by harmonic oscillator basis®. To do 
this, let us introduce annhilation and creation operators for the harmonic 


oscillator: 
1 w 
t= a+a'), = i,/—(a—a'). 9 
geti), pai Sa-al) (9) 
The annihilation operator in the number operator basis n is 
a= S>VnFI\n)(n +I. (10) 
n=0 
Let us truncate the Hilbert space such that it is spanned by |0),|1),--- ,|A— 
1) and replace the annihilation operator by 
A-2 
a > ay:= >) VnFI\n)(n +1). (11) 
n=0 


Then we consider the Hamiltonian 


H(z,p) > Hy, := H(xa,pa), (12) 
where 
1 w 
t . t 
lA = aa + Gay), =í a an). 13 
A Jaa! A A) PA 5 | A A) (13) 


We can map the truncated system to a spin system as follows. First we 
consider a binary representation of n: 


n= by_12*—-1 + by_92* 7 +--+ + b92° (K =log A), (14) 
and express the basis |n) as 
In) = |bo) -+ bx). (15) 
Then each component of the annihilation operator is given by 
|r) (n + 1] = 840 (|b) (bel), (16) 
which can be mapped to a spin system via 
1-Z 142 
o=, was, 
X+iY X—iY 
oa- mo- (17) 


This is a qubit description of the single particle quantum mechanics. It is 
straightforward to generalize this method to multi-particle case. 


°See e.g. the papers® 9 for other basis. 
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3.3. Gauge field 


The case of gauge fields is more complicated and depends on spacetime di- 
mensions and and boundary conditions. This is realted to the fact that in 
gauge theory, there are not only physical states but also unphysical states, 
and whether or not the dimension of the physical state space is infinite 
depends on situations. To see this concretely, let us consider a 1+ 1 dimen- 
sional quantum electrodynamics coupled to charge-q Dirac fermion. We 
will call this theory the charge-q Schwinger model or simply the Schwinger 
model!!:!?, When the space is cut into a lattice, the Hamiltonian of this 
theory is given by'8 


i a 
H eae lee ices = XE a 


N-1 ga N= pN? 
Tat 
+m y (—1) XnXn + D 5 (1, + x] ; (18) 


where g is gauge coupling and @ is theta angle. The lattice gauge field 
operators (Øn, Ln) satisfy the canonical commutation relation 


(dm; Ln] = imn. (19) 
Now the state space associated with the gauge field is infinite dimensional 
since Øn is bosonic. However, this is the state space including unphysical 
states and the physical states are restricted by the Gauss law: 


Ln — Ln- = 4 (xh = ar (20) 


This relates L,, to its neighbours and the fermionic operator. When we 
take the open boundary condition, we can rewrite L» purely in terms of 
the fermionic operator and ¢, can also be absorbed by a gauge transfor- 
mation. Therefore the Schwinger model with the Gauss law is rewritten as 
a fermionic system and we can further map it to a spin system by e.g. the 
Jordan-Wigner transformation. Thus we can directly apply the quantum 
algorithm to the Schwinger model. 


4. Quantum algorithm to prepare vacuum 


Here we would like to construct the ground state by a quantum algorithm. 
While there are several known algorithms, here we use an algorithm called 
adiabatic state preparation. Suppose that we would like to constract the 
ground state of the Hamiltonian Htarget of the target system. First, we 
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prepare an initial Hamiltonian Hp whose ground state |GSpo) is known and 
non-degenerate. Therefore Ho is practically taken to be a Hamiltonian of 
a simple system. Next, we introduce a time-dependent Hamiltonian H4 (t) 
satisfying 


H(0) = Ho, Hy(T) = target: (21) 


and we take H,(t) to change more slowly for larger T. Then we use the 
adiabatic theorem. That is, if the ground state of H(t) has non-degenerate 
for any t, then the ground state |GS) of Htarget is obtained by the following 
time evolution: 


T 
IGS) = Jim Texp € f anat )1GSo) (22) 


In actual simulations, the time evolution operator is approximated by tak- 
ing finite T and discretizing the integral in the exponent. This gives errors 
due to the fact that the constructed state in pratical simulation is not ex- 
actly the ground state. Using the approximate ground state constructed 
in this way, we can calculate the expectation value under the ground state 
approximately. 


5. Recent simulation results 


Using the adiabatic approximation, we can construct the ground state of the 
Schwinger model and calculate various physical quantities. Here we present 
a part of numerical results of our simulations! 3 (see e.g. the papers 1418 
for other simulations of the Schwinger model). The Schwinger model has 
been analyzed by various methods. Since it is 1+ 1 dimensional, a powerful 
anaytic approach by bosonization is available. In particular, it is known 
that we can exaclty solve the massless case (m = 0) 19 22. When the fermion 
mass m is nonzero, there is no known exact solution but we can get a good 
approximation by mass perturbation theory?*:?* for small m. Regarding 
the conventional Monte Carlo approach, it is known to be difficult when 0 
is not small due to the sign problem while small 0 region is accessible. 
There is one thing before presenting the results: real quantum computer 
has errors due to interactions with einvironment. Therefore we need to take 
error correction into account in order to obtain error-free results. However, 
it is known that this requires huge computational resources and this is a 
major obstacle in development on the technology side. To test quantum 
algorithms and estimate the computational resources required, people often 
use a tool called a (classical) simulator. A simulator is a tool to simulate 
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Fig. 1. Chiral condensate in the q = 1 Schwinger model after thermodynamic and 
continuum limits against the gauge coupling g for the massless case (m = 0). The red 
straight line represents the exact result. 


a quantum computer by a classical computer and we can use almost the 
same code to run a real quantum computer. Here we use a simulator called 
“qasm simulator” provided by IBM. 


5.1. The q = 1 Schwinger model without probe charges 


First, we focus on the q = 1 case, which has been investigated well. Here 
we consider a quantity called the chiral condensate!. In the continuum 
theory, this quantity is defined as the expectation value of the fermion 
mass operator y(x) under the groundstate: 


(GS|h(x)|GS). (23) 


We compute the chiral condensate for various values of the parameters 
(a, N), and take the infinite volume limit N — oo and the continuous limit 
a — 0. In fig. 1, we plot the chiral condensate for the massless case against 
g and find agreement with the exact result. In Fig. 2, we plot the chiral 
condensate? at g = 1 against m for 6 = 0 and @ = 32/5. When m is 
small, the simulation and mass perturbation results agree, whereas when 
m is large, they deviate from each other. This is intepreted as violation 
of mass petruatation theory for non-small m. In the context of quantum 
simulations of field theory, this is the first result in which the continuous 
limit is taken seriously. 


dWe regularize it by subtracting the free result. 
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Fig. 2. The chiral condensate in the q = 1 Schwinger model for massive case at g = 1 
for 0 = 0 (blue x’s) and @ = 37/5 (green circles). The lines denote the result obtained 
by the mass perturbation theory. 


5.2. With probe charges 


Next let us consider the q # 1 case. We also put two probe charges +q, 
and consider the potential between them. According to the mass pertur- 
bation theory ?ř, the potential changes its qualitative bahavior as changing 
the parameters: 


V(é)=o8 for gl >1, (24) 
where £ is the distance between the probes and ø is the string tension 


eqg 


0 + 2Tqp 0 
273/2 cos 


q q 


+ O(m?). (25) 


At first sight, this seems to show a confinement behavior but note that the 
tension ø is not always positive. In particular, we find ø = 0 when qp is an 
integer multiple of g. This implies screening of the charges. We can also see 
that ø can be sometimes negative when qp/q is not an integer. The reason 
why this can happen is closely related to the generalized global symmetry, 
but we will not go into details here (see e.g. the papers*:?"), 

To calculate the potential, we measure the difference of the ground 
state energies in the theories with and without the probe charges. In the 
(1 + 1)-dimensional U(1) gauge theory, the effect of the probe charges is 
taken into account by replacing 0 by a space-dependent 6, such that its 
value suddenly changes by 27q,. In this way we can compute the potential. 
Figure 3 shows the potential for q = 3 and qp = —1 for various values of 
0. The simulation results show a linear behavior for all values of 6. while 
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Ve)/g 


Fig. 3. The potential V(¢)/g between the probe charges (qp = —1) is plotted in the 
q = 3 Schwinger model for various values of 0 = 09. Here we take m = 0.15, N = 25 and 
ga = 0.40. 


signs of its slopes change as changing 0 at the same timing as the result of 
the mass perturbation theory. 


6. Discussion 


In this article, we discussed applications of quantum computation to nu- 
merical simulations of quantum field theories. This topic is still in its 
infancy and there are still many things to explore. There would be future 
developments both in hardware and software aspects. While we presented 
recent numerical results obtained by a classical simulator, we need use real 
quantum computers to have speed up of quantum algorithms. 

Regarding the hardware perspective, although an ideal quantum com- 
puter is a large machine with error correction, machines available in near 
future are NISQ (Noisy Intermediate-Scale Quantum device) ?°, where er- 
rors are not negligible. For such a machine, it is desirable to apply an 
algorithm that uses as few quantum gates as possible. While the adiabatic 
state preparation used in this article is guaranteed to be correct under the 
assumptions of the adiabatic theorem, it is not suitable for NISQ devise 
because it uses many quantum gates to approximate the time evolution. 
The algorithm expected to be suitable for NISQ is the one that uses only 
a small number of gates such as classical-quantum hybrid algorithm. An 
algorithm to construct the ground state along such a direction is based on 
the variational method. 


29 


Regarding the software perspective, we have not yet established an effi- 
cient way to put most field theories on a quantum computer. In particular, 
in other gauge theories above 2+ 1 dimensions such as quantum chromody- 
namics, it is hard to establish to put theories on quantum computer while 
respecting gauge symmetry. So far, a simple lattice regularization has been 
applied but this may not be best for some theories®°. In addition, there are 
many interesting states other than the ground state in QFT, but this as- 
pect has not been explored well yet in the context of quantum computation. 
One of the most interesting states in high energy physics are ones treated 
in the scattering problem, in which multiple particles are moving with a 
fixed momentum. If we can construct such a state, then we would be able 
to understand what is “actually” happening in the accelerator experiments. 
In this article, we have focused on QFT, but it is also interesting to con- 
sider applications to quantum gravity such as string theory?’. I hope that 
we will be able to literally see time evolution of our universe on quantum 
computers in the future. 
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Multi-Soliton scattering of the Anti-Self-Dual Yang-Mills 
Equations in 4-dimensional split signature 
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We construct the ASDYM 1-solitons and multi-solitons for split signature and 
interpret them as soliton walls. We show that the gauge group is G = SU(2) 
for the entire intersecting soliton walls, and SU(3) for each soliton walls in the 
asymptotic region. This is joint research partially with Masashi Hamanaka, 
Claire R. Gilson, and Jonathan Nimmo. 


Keywords: Anti-self-dual Yang-Mills equation; Soliton scattering; Exact solv- 
able model; N=2 Strings. 


1. Introduction 


The anti-self-dual Yang-Mills (ASDYM) system is notable as an exactly 
solvable model and play an indispensable role in many different fields of 
theoretical and mathematical physics, such as quantum field theory, twistor 
theory®, and integrable systems®. Especially the ASDYM equations reveal 
the essence of the lower-dimensional integrable systems because almost all 
of the soliton equations are the dimensionally reduced equations of ASDYM 
according to the Ward conjecture!?. Among these lower-dimensional soli- 
ton equations, dimensionally reduced from the 4-dimensional split signature 
(+, +, —, —) (or called the Ultrahyperbolic® spacetime) are the majority. 
On the other hand, the ASDYM equations in the 4-dimensional split sig- 
nature is the EOM of effective action for N = 2 open string theories!! and 
seem to have good compatibility with the lower-dimensional integrable sys- 
tems. Although most physicists didn’t pay too much attention on N = 2 
strings and there has been little progress in this direction for a long time, we 
still believe that the connections between the lower-dimensional integrable 
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systems and N = 2 open string theories are worth studying and not just for 
seeking the physical reality. 

In (1+1)-dimensional integrable systems, the KdV equation —4u; + 
Gut, + Urre = 0 might be the most notable one among the exactly solvable 
models and the typical KdV 1-soliton is in the form of 


2 


o 
u(a,t) = 25-2 (log coshX) = 2n*sech?X, X:=Ko+n7t+6 (1) 


which is a travelling wave solution with a constant velocity —K? 


and a 
constant amplitude 27, and hence it preserves its shape over time. In 
addition, u(x,t) is an even function symmetric with respect to X = 0 and 
U, Ux, Ure, Urre > 0 as x — too. Therefore, the distribution of u(x,t) 
is localized within a region centered on X = 0. On the other hand, 
the KdV equation is also known for having infinitely conserved quanti- 


ties®, more specifically, the mass f udx, the momentum f u?dz, the energy 
J [2u3 — (uz)?] da, ..., and so on. In particular, the energy density 


2u? — (uz)? = 16° (2sech°X —sech*X), X:i=Kat+nt+d (2) 


possesses the behavior of 1-soliton as well. Inspired by this, we simply 
consider the pure Yang-Mills action density TrF,,F"” as an analogue of 
energy density, and try to figure out whether the 1-solitonic behavior like 
(2) exists for the ASDYM equations in 4-dimensions or not. Fortunately, 
the answer is yes and our result* is 


Tr F FH” œ (2sech? X — 3sech*X) , (3) 


where X is a nonhomogeneous linear function of real coordinates z1, x”, x’, 
xt. Especially, the gauge group is G = SU(2)* for the split signature and 
hence it means our ASDYM 1-soliton could be some candidate of physical 
object in N = 2 open string theories. 

One more well-known fact is that the KdV multi-soliton 


Ta 
u(a,t) = 255 = =, (log mn) (4) 
can be represented elegantly as the Wronskian determinant 
f,O foo coe TAS 
(1) (1) (1) 
fi to a fn 
= Wr( fi, f2,- fn) = : : és : , (5) 
ford ped S pe- 
(m) _ O fi 3 
where f; : fi = coshX;, Xi = kiz + Ket + ôi (6) 


Oa? 
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Now a natural question to ask is whether the Wronskian type multi-solitons 
exist for the ASDYM equations or not. If the answer is yes and the gauge 
group can be verified to be unitary, then what is interpretation for such 
ASDYM multi-solitons in N = 2 string theories. 


2. J-matrix formulation of ASDYM and Darboux 
transformation 


Before going ahead to the main topic, we need to introduce some prior 
knowledge. Firstly, the ASDYM equations on various real spaces can be 
unified into 4-dimensional complex flat spacetime with the metric ds? = 


2(dzdz — dwdw) and we can easily get the split signature (+, +, —, —) by 
imposing some conditions on the complex coordinates, for instance 
1 Li eB p2 4 
oa es oo. to , £! r? r’, xt ER. (7) 
wz V2 \—(r? — r4) rt- r 


Here we set the gauge group to be G = GL(N, C) in general. The complex 
representation of ASDYM equations are 
Faw = 0; Fxg =0, Les wa = 0 (8) 


which can be cast in a gauge independent formulation, called the Yang 


equation +13 


ðz[(3-J) J] — dz{(0.J)I~"] = 0, (9) 
where J is an N x N matrix called Yang’s J-matrix. An advantage of this 
formulation is that the anti-self-dual (ASD) gauge fields can be reformulated 
by the decomposition of J = h~th as 
A,=—(dzh)h7', Ay =—(Owh)h-', Ag=—(Ozh)h, Ag =—(Ozh)A-+. (10) 
Sometimes we use a convenient gauge h=1to simplify the gauge fields as 

Ay==(6,3) I~; Aw = — (3w J)JT!, Az =0, Ag = 0. (11) 
Now we introduce a novel Lax representation !? of ASDYM equations : 
{ L(ġ) := [Ow — (Qu) I~*]6 — (0z4)¢ = 0 
M($) := [0 — (0.J)I~*]6 — (3z) = 0’ 
where the spectral parameter ¢ is any right-action N x N constant matrix 
rather than scalar. Under the compatible condition L(M(¢)) — M(L(¢)) = 


0, the ASDYM equations (Yang equation) can be derived from the linear 
system (12). Furthermore, if we define a Darboux Transformation !° as 


$=- pA, J=—-bAd71J, (13) 


(G=GL(N,C)) (12) 
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where $(¢) denotes the general solution of (12) and we use a new notation 
(A) to denote a specified solution of (12). Then the linear system (12) is 
form invariant under the Darboux transformation (13). In other words, 


EO) = Ow — Oy DING — (Beoie = 0 aa 


M(8) = [ð; — (8-F)I~"]9 — (zde =0 ` 
Now we can choose a seed solution J = Jı of the Yang equation (9) and 
substitute it into (12) to solve a specified solution 41 (Ai). After 1-iteration 
of the Darboux transformation (13), we get a new solution J = Jj. By 


6 


repeating the same process’, we can obtain a series of J-matrices. 


Seed solution 


a `~ D D D D D D 
Aj She SA SS hpa OS. 


and the J-matrix J;41 can be written concisely in terms of the Wronskian® 
type quasideterminant (Cf: (5)) : 


pı ors Yn 1 
pidi PA Pnn 0 


In41 = : : . : : Ji, (15) 
WAP d2AZ- vn Ar [0] 

where w;(A;) denote n specified solutions of (12) with matrix size N x N. 

We use the term quasi-Wronskian to call it for short. In fact, since the 

elements of J; are noncommutative, the quasiderminant can be considered 

roughly as a noncommutative version of determinant. By the definition? , 

we can decompose (15) into 4 blocks as 


AnnxnN BnnxN 


which shows that J; are all N x N matrices and G = GL(N,C). 
If we choose a seed solution J with det(J) is constant, by the Darboux 


transformation (13) we have det(J) = det(—A) det(J) is also a constant. 
Applying the Jacobi’s formula 


=D-CA'B (16) 


“det A(t) =Tr baay = det A(t) Tr(a a7) 


to J and comparing with (11), we find that the gauge fields are all traceless. 
That is, any J-matrices (15) generated by a seed solution Jı with constant 
determinant are belonging to G = SL(N, C) gauge theory. In fact, we will 
always set the seed solution Jı to be identity matrix and only consider 
N = 2,3 cases in the next sections. 
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3. ASDYM 1-Soliton and Multi-Soliton for G = SU(2) 


Now we set the seed solution J; to be 2 x 2 identity matrix (G = SL(2, C)) 
so that the linear system (12) reduces to 


L(¢) = Owh — (0z0)¢ =0 (18) 
M(4) = 0.6 — (3z) = 0 ` 
For the split signature (+,+,—,—), we can take the reality condition (7) 


and find a lovely solution of (18) as 


—be—" Te 


È ead 
y= ( Oe ee I ) w.r.t. a specific spectral parameter A = iG 2) , (19) 


where L = Z [Aa + B)at + (AB — a)x? + (Aa — B)a? + (AB + a)zt], 
a, b, a, 8, AEC. 

After 1 iteration of the Darboux transformation, we obtain a candidate 
of 1-soliton solution : 


= -vAy7! 


J2 = 2 
2= yA [0] 
ai (ATA MT Yao 
det (w) (A — Ajabe h) Alaj Zer tE + Albe EtE) 


The resulting action density4 of (20) is 


TrF uu F" = 8 |(aB — Tp) (A — d)] ? (2sech” X — 3sech* X) 


X =L+TL +log(|a| /|b|) ey 


We find that the distribution of action density is localized on a 3- 
dimensional hyperplane X = 0. Therefore, this kind of solitons can be 
interpreted as codimensional 1 soliton and we use the term soliton wall to 
distinguish them from the domain wall. 

Now we can prepare n different solutions w;(A;) of (18) as follows : 

Li Be hi : 
pi = ae “Oh w.r.t. spectral parameters A; = (y 2) : 
Li= a [Ass + Bi)" + (As Bi — 04) 2? + (Aiai — Bi)aP + (Aib: + on) a") , 
Qi, bi, Qi, Bi, Ai E C, 7=1,2,---,n. 


(22) 


After n-iterations of the Darboux transformation, we obtain a candidate of 
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n-soliton solution : 


pı Yo cs Un 1 
pidi PA = UnAn 0 
an ae (23) 


Jn+1 = . . . . . 
WA? p2A3 -+> YrAt [0] 
which satisfies the property ®® 


sing =] APL, 1:2 2 identity matrix. (24) 
i=l 


Ina ged) 


n 


This fact implies that the gauge fields in 4-dimensional split signature 


ay = ay =F [3 Jn) Int + (OsIn41) Spa] 
AŞTD = AÇHÐ = Z [(ðzJny1) iy + (Onda de | 


n> 


(25) 


are all anti-hermitian and therefore the gauge group is G = SU(2). 

Now a natural question to ask is whether the n-soliton solution (23) 
gives rise to n intersecting soliton walls or not. Let us use a quite similar 
technique as mentioned in reference 7 to discuss the asymptotic action 
density of (23) rather than calculating the action density directly by (25). 
First of all, we fix an I € {1,2,...,n} and consider a comoving frame related 
to the J-th 1-soliton solution : 


Ly T =br 
D YDA (gO). yD | we” bre 
J. = A ’ f S Es 26 
2 Vn (by ) Y ta Tr el! ( ) 
whose action density is 
Tef ppd 0) = [(arBy _ arBr)(Ar = Ar) 5 (2sech° Xr = 3sech* X7) (27) 
Xr = Lr + Lr + log(|ar| / |b7|) 
More precisely, we define r := y/ (£1)? + (x2)? + (a3)? + (x4)? and consider 
the asymptotic limit r + oo such that 
X, is a finite real number (28) 
Xiizr > £00 (ie. TrF F C#D — 0) 


By some mathematical techniques®® of the quasideterminant, we find that 


Too 


hag SF -TOAT ODO, n> 9, (29) 
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where DP is a constant matrix and doesn’t affect the gauge fields, and 


n n 


I] Qr-Adare™ TT Ar- Ai) bret 
(i) ie let 
= I] (Ar — Ài) bre" II (Àr — i) Gre z 
i=1 tA i=1,iAI 
TO. as Xj i447 > +00 (30) 
I] (Ar — ji) are”? II (Àr — di) by e`% 
(ii) fare rcs E 
= I] (Ar = ài) bre! II (Ar — Ai) ar ell 
i=1iAl i=1,iAl 


as Xi iZl =? CO 


Comparing the asymptotic n-solion solution (29), (30) with -th 1-soliton 
solution (26), we find that A is in the same form as YẸ up to constant 
factors. These factors lead to a position shift from the principal peak of 
action density (27), called phase shift. More specifically, the action density 
of n-soliton solution in the asymptotic region 


TrF es [(arbr — Tbr) (Àr — xl (2sech? X; = 3sech" X ) (31) 


behaves like the action density of J-th 1-soliton (27), where Xr := X,;+Ay 
and the phase shift 


fa =+, Xr = Ho (32) 


t 
eb? := —-1, Xiizr > —0 


is real-valued. Since J is any positive integer from 1 to n, and for every I the 
n-soliton solution gives rise to a soliton wall in the asymptotic region, we 
can conclude that the n-soliton solution can be interpreted as n intersecting 
soliton walls in the entire region. Furthermore, the n intersecting soliton 
walls can be embedded into G = SU(2) gauge theory in 4-dimensional split 
signature and therefore they could be interpreted as n intersecting branes 
in open N = 2 string theories. 


4. Reduction to (1+1)-dimensional real space 


To make the discussion more clearly, let us take 2-dimensional spacetime 


and 3-soliton scattering for instance. We can impose the condition x! = t, 


x? = zt = 0, £? = x on the spacetime coordinates such that 


Li = + (Aca + bijt + (Aiai — Bix), i= 1,2,3 (33) 
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and X; = Li + Li + log |a;/b;|. Then we fix an I € {1,2,3} and choose a 
complex number ¢; such that Ly = ¢; which implies 
Nia — Bi ) (a — nauk) ' 
L = | 2 jg p v2 | ZA AT le iI. 34 
(3 = a Arar — Br 7 d 
Now the setup of the comoving frame related to the J-th 1-soliton is com- 
pleted because 


when t — +00 (35) 


Xr = finite 
XiiAl — roo or F 


This setup implies that the action density of i-th 1-soliton (i 4 I) behaves 
as Tr Fy) t GAD ~ 2sech? X; — 3sech* X; — 0 as t > +00. On the other 
hand, the action density of 3-soliton in the asymptotic region behaves as 


= a se LA) 
Par ~ e A n X,—4 * A iat 5) 
Xr+ as t+ —oo 
In fact, Aw) = ZAG and it depends on 2371 = 4 choices of asymptotic 
regions. For example one of the choices is 
AW) = -A = — log t= à — log Ar As ; (37) 
Ay = A2 Ai = À3 
£ = — A2 — 
AS) = -AP = + log àA log | às (38) 
A2 = Al A3 ai A3 
y = A3 — A A3 — À 
AD = -AP = + log B=] + log] 2—2]. 39 
3 3 8 NEN 8 Ie he (39) 


5. An example of ASDYM 1-Soliton for G = SU(3) and 
Multi-Soliton scattering 


If we consider the linear system (18) for G = SL(3,C), these is a special 
class of solutions 7);(A;) that give rise to soliton walls® and behave like the 
G = SU(2) cases as well. Firstly, the candidate of G = SL(3,C) 1-soliton 
solution can be constructed by 1-iteration of the Darboux transformation 


aje’ bie T: 0 
pi := | -bjeli Tje ceti 
a = -L; Li 
d = piip t, ` ; 0 B ne f 
Av= {0% 0 (40) 
0 0 A; 


Ly= A [ioe + Bia? + (Ai Bi — ai)? + (Aiai — Bia? + (Aibi + 04) 24] , 
Qi, biy Ci, 4, Bi, Ai € C, (= 1,2, sasa (Us 
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By direct calculation, we obtain the action density of i-th 1-soliton as 
Trk pd” @) — 8 [(aiB; = Tibi) (Ai = z) (2sech? X; = 3sech*X;) 


7 , (41 
X; = Ly +T; + 4log [as /(bl? + les!) em 


which can be interpreted as i-th soliton wall as discussed in G = SU(2) case. 
Substituting ~; of (40) into (23), we get a version of n-soliton Jn41 for 
G = SL(3, C). By some mathematical techniques® of the quasideterminant 
and considering a comoving frame related to the J-th 1-soliton, we find that 


Inga 2S FD, =: GA (GO) 1D, n> 2, (42) 


Q) 


where D is a constant matrix and vu! is in the form of 


A=) be ie 
_ i=1,iŻI 
Are: Bre™™ 0 a (F) 
Y = | -Bre Are! Cre |, a i (43) 
0 =C re ™ Are” Cr = ji (Ar — A) CI 


=1if¢l 
al) ) = (Ais Ai) 


Comparing (43) with (40) and (41), we can conclude that the action density 
of n-soliton in the asymptotic region 


TrF,,F* "2S 8 [(arBy -Trb Ad] (2sech? X;, = 3sech" X,) (44) 


behaves like the action density of J-th 1-soliton, where X 1 := Xr + ^z and 
the phase shift is 


TL [ar-a dor? + ler) 
SORE en (45) 
i ar-a or? + TE, ar-a] er 
i=1,iŻI i=1,iŻI 
which is real-valued and depends on 2”~1 choices of AR AD?) := (Ai, Ad) 
Therefore, the n-soliton can be interpreted as n-intersecting soliton walls 
for G = SL(3,C). On the other hand, the asymptotic form TAI of the 


n-soliton solution Jn+41 satisfies the property ê 


1 
Ar = 5 log 


SCE FU I vs 
I aS ITIP i (46) 


which implies that the gauge fields given by JE wA are all anti-hermitian. 
Therefore for each single soliton wall in the asymptotic region, the gauge 
group is in fact G = SU(3). 
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6. Conclusion 


In this paper, we found a class of ASDYM 1-solitons in 4-dimensional split 
signature for G = SU(2) and SU(3), respectively. The resulting action den- 
sities are in the same form as TrF,,F"" œ (2sech? X — 3sech*X ) which 
can be interpreted as the soliton walls. After n-iterations of the Darboux 
transformation, we obtain the ASDYM n-soliton and interpret it as n inter- 
secting soliton walls for G = SU(2) and SL(3,C), respectively. This fact is 
a well-known feature for the KdV multi-solitons, but a new insight for the 
ASDYM multi-solitons. On the other hand, the n intersecting soliton walls 
can be embedded into G = SU(2) gauge theory and hence they could be 
interpreted as n intersecting branes in N = 2 open string theories. There- 
fore, to understand the role of such physical objects that play in N = 2 
open string theories would be an interesting future work, and the relation- 
ship between N = 2 open string theories and lower-dimensional integrable 
systems is also worth studying. 
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We first study the WKB analysis of the third order ODE, which can be regarded 
as the quantized Seiberg-Witten curve of the (A2, An )-type Argyres-Douglas 
theory in the Nekrasov-Shatashvili limit of Omega background. We then derive 
thermodynamic Bethe ansatz (TBA) equations satisfied by the Y-functions 
from the solutions of the ODE, and identify the Y-function with the WKB 
period. For the (A2, A2)-type ODE, we study the process of wall-crossing of 
the TBA equation from the minimal chamber to the maximal chamber. 


Keywords: TBA equations, WKB period, Wall-crossing, ODE/IM correspon- 
dence. 


1. Introduction 


The WKB period of the Schrödinger equation is expanded as the formal 
asymptotic series in the Plank constant, which thus needs to be Borel re- 
summed. Recently, a connection between the exact WKB period and the 
thermodynamic Bethe ansatz (TBA) equations of the quantum integral 
model has been noticed for the Schrödinger equation with arbitrary poly- 
nomial potential!:?. In particular, the Y-function of the TBA equations is 
identified to the exponential of the Borel resummed WKB period, which 
share the same asymptotic behaviors and the discontinuity on the com- 
plex fA-plane. The TBA equations thus provide a solution to the Voros’ 
Riemann-Hilbert problem of the exact WKB periods?. Moreover, the TBA 
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equations together with the exact quantization condition lead to the exact 
energy spectrum of the one-dimensional Quantum Mechanics +4. 

The Schrödinger-type equation with polynomial potential also appears 
as the quantized Seiberg-Witten curve of the Argyres-Douglas (AD) theo- 
ries in the Nekrasov-Shatashvili (NS) limit of the Q background® 7, where 
the non-zero parameter of the Q-background plays the role of Planck con- 
stant. The TBA equations satisfied by the WKB period (quantum period) 
thus becomes very useful to compute the exact BPS spectrum, whose wall- 
crossing was observed when the moduli space parameters are varied. Due 
to the wall-crossing, the form of the TBA equations also changed, which 
we will denote by the wall-crossing of the TBA equations. 

To study the exact BPS spectrum of the higher rank AD theories, it 
is important to explore the relations between the exact WKB periods of 
the higher order ODE and the TBA equations in the quantum integrable 
model, whose wall-crossing will be useful to understand the wall-crossing 
phenomena of more general non-perturbative supersymmetric gauge theo- 
ries. In®, we have focused on the higher order ODE with quadratic poten- 
tial, which can be regarded as the quantized SW curve of the (A,, A1)-type 
AD theory. The corresponding TBA equations are found to be (A,, A1) 
type TBA equations. Moreover, the wall-crossing of TBA equations was 
observed for the third order ODE with cubic potential. The wall-crossing 
of the TBA equations for more general third order ODE has been studied 
in®. For the monomial type potential, the TBA equations obtained from 
the (A2, A2)-type and (Ag, A3)-type ODE become the D4-type and E¢-type 
TBA equations, respectively. 

This proceeding is organized as follows. In section 2, we study the WKB 
analysis of the third order ODE and introduce the Borel resummation of 
the WKB period. In section 3, we derived the TBA equations satisfied 
by the Y-functions, and identify them with the exact WKB periods in 
the minimal chamber. In section 4, we study the wall-crossing of the TBA 
equations. The TBA equation in the maximal chamber will be constructed. 
The section 5 is is devoted to conclusions and discussion. 


2. WKB analysis of the third order ODE 
We consider the third order ODE on the complex plane: 
d8 
(Cm + p(a)) p(z) =0, p(x) = woe? +u" +- tuner (1) 


where e and u; (i = 0,--- , N +1) are complex parameters. This ODE can 
be regarded as the quantized SW curve of the (Az, Ay )-type AD theory in 
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the NS limit of the Q-background. We will denote this ODE by (A2, An)- 
type ODE in this paper. 

The WKB analysis begins by imposing the following ansatz of the wave 
function 


1 x 
v(e) =e (4 [ Pear), (2) 
where P(x) satisfies the Ricatti equation 
p(x) + P? + 3ePP’ + 2 P®?) =0. (3) 
The solution can be expanded in power series of e€ 
z 5 "Dy (x), (4) 
n=0 
which can be determined recursively by 
Po = a 
=i ni 
Pn 32 ~~ [0 2 iPi + 5 NE Da i—jPiPj 
-3p a=1 j=0 (5) 


+3 ree + p-a] » nei. 


It is natural to regard Pdz as the meromorphic differential on the WKB 
curve % 


y? + p(x) =0, (6) 
which is nothing but the SW curve of (A2, Aw) AD theory. We introduce 
the period of Pdz along the one-cycle y: 


" J Pla)de = F eT, (7) 
Y n=0 


which will be denoted by the WKB period. Since the basis of the mero- 
morphic differentials on the WKB curve can be generated by 
N+1-i 
x 


Ou, y dx = “3 ee a= 1,2, (8) 


we can obtain the quantum correction of the WKB periods by acting the 
differential operators of u; on the periods of y*: 


2 
m =y Ol Tig, (9) 
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where 
(Ha), = f yde. (10) 
ah 


The operator is called Picard-Fuchs operators, whose details can be found 
in®9, 

Since the WKB periods is an asymptotic formal series in €, we need to 
perform the Borel transform 


BII = >> 


n>0 


1 
yen 

which Laplace transform along the direction y is called the Borel resum- 
mation 


1 ooet? i 

sosi BE, (12) 
The WKB period is Borel summable when the Borel transform converges. 
There arises a discontinuity for the resummed WKB period 


disc I (€) = 85+ (Iy) (€) — sy- (IL,) (€)) 


HA (s(IL,)(e'?¥e) — s(IL,)(e"?-e)), 


(13) 


when B[IL,] is singular along the direction y. By using the Borel-Pade 
technique, we are able to compute the singularity structure of the Borel 
transform numerically. In Fig.1, we plot the singularity structure of the 
Borel transform of the WKB periods for the potential p(x) = —x? + 7x +6, 
which will be compared to the TBA results later. 


3. TBA/WKEB correspondence of the third order ODE 


In this section, we derive the TBA equation derived from the solutions of 
the ODE, and identify the Y-function with the WKB periods. 
The ODE (1) is invariant under the rotation: 


cow tt, u > wtu, i=0,1,...,N +1 (14) 


with w = e?"i/(N+4), Since the solution to the ODE is changed, the rotation 
is very useful to generate the solutions of the ODE 


bw te, {wui}; 6) = O(a, {ui}; ee), (15) 
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1 
N 
2 ooo 


Fig. 1. The singularity structure of the Borel transformed WKB periods I}, , (blue) 
and Iy; ı (yellow). The Borel-Padé technique is applied to order e160 terms of the 
formal power series computed by using the Picard-Fuchs operators. Here the potential 
is p(x) = =x? + 7a +6. The figure was first presented in®. 


where we have rewrite the rotation in terms of e. At the infinity of the real 
and positive axis, the subdominant solution of the ODE behaves as 


go (a as 


=) |] 00, larga] < (16) 


T 
N+4’ 
where we have set ug = 1 without loss of generality. The rotated solution 
Qp defined by 


Ti 


r(x, {ui}; €) = polz, {ui}; e” 


are the subdominant solution in the sector 


Re) (17) 


Qk T 
Sr = fx € C; larg z| ana (18) 
Using these solutions, we introduce the T-functions 
Tok = W[d-1, 60, ¢1]-* = 1, Tik = W[d-1, 0, 9k41)", 
Ton = W[oo, k41 Orpa] 07, Tse = Wide, desis Orpa]! = 1, 
(19) 
where fl*l({u;};€) = f({u;j};e3 *e) and WH] is the Wronskian 
fra fra frs 
W [fies fka fka] = det | 3x frkı 3x fro Ox frka |- (20) 


OF fri O? fies O? fks 


48 


We introduce the Y-function by the cross ratios of the T-functions 


Ta—1.kTa 
Yap = EELE og 1, k=1,...,N. (21) 
Ta,k-1Ta,k+1 
By using the identities of the Wronskian, it is to find these Y-function 
satisfy the (Az, Ay )-type Y-system 9: 
1+ a—1,k)(1 4 


yi ee ee (22) 
i (1+ Yak- (1 Ga Yak+1) 


Using the WKB approximation of the solutions ¢, and the Stokes 
graphs, we find the asymptotics of the Y-function are determined by the 
WKB period. Together with the discontinuity structure, we find 


[a—k] 
log Yak = EM ; Va,k = Y2-k,k Sere Yatl—k,k (23) 


a=1,2, k=1,...,N, 


where 7, is the one-cycle encircling the branch points x,_ 1 anticlockwise 
and x; clockwise, respectively, on l-th and (J + 1)-th sheets of the WKB 
curve. In the following of this section, we will test this relation numerically. 

Based on the asymptotics of the Y-functions, it is easy to convert the 
Y-system into TBA equations 


log ¥i,4(0 — idp) = [mi ple? + K * Li, — Kp k1 * L1p-1 — Kk kpi * L141; 
(24) 


2 and ¢; = arg(™m1,x). 


where x denotes the convolution, Ma, k = ga (ka) 
The kernel is defined by 


a, 


1 4/3 coshé 
K(0) = ————__.,_ K 0)= K(@-1 — i 2 
(0) = a Keel) = KO — ilhe — rs). (25) 
At small €e = e7? , we can expand the TBA equations by 
log Yı (0) = mi pe? + 5 ie (26) 
n=1 


where the coefficient m?) can be computed by using the solutions of the 


TBA equations. In table 3, we compare the e-expansion of the two hand 
sides of (23) for the potential p(x) = —x? + 7x + 6. 


The e-expansion of the two hand sides of (23) for p(x) = ~x? + 7a + 69. 
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n n-1 
n% min 


n 
1,1 

2 0.21721574367 | 0.21721574362 

6 || —1.5195679452 | —1.5195679452 

8 || —20.486617772 | —20.486617762 

12 20065.209702 20065.206052 

14 1160395.6762 1160393.4227 


Moreover, we can express the left hand side of (23) by using the TBA 
equations: 


log Y11(0) = mı 1e? +f dð K 0 — 6g’ ia a ib1)L11(6') 


=00 


-f do’ K (0 — 0" io a0, 
log Yı 2(0 — a = m2 +f do’ K (0 — 0 + ip — 


3 )L12(6') 


-J do’ K(0— 6 + ids — Zap : 


3 

(27) 
When arg(m 1,1) = arg(™m1,2) = §, ie. the case studied in Fig.1 and table 
3, we find the poles of the kernel of the TBA equations are located along 
the directions 0 = —*,-£, Bi 5z for log Y1,ı(0) and 0 = —3, — 32, 5-3 for 
log Y1,2(0 — 5). Comparing with the singularity structure of Fig.1, we find 
these poles are nothing but the discontinuity of the Borel resummed WKB 
period of the right hand side of (23). 


4. Wall-crossing of the TBA equations 


In this section, we study the wall-crossing of the TBA equations for the 
case of (Az, Az) ODE. We parameterize the zeros, xo(t), xı(t) and 22(t), 
of the potential in the following way 


ao(t)=3-—t, a(t) =—14+V3it, z(t) =-—2+t-vV3it, O<t<1 
(28) 
and denote the potential by p(x, t) = (x — xo(t))(x — zı (t))(a@ — x2(t)). For 
a given value t, we can compute the classical periods and mass of the the 
corresponding Y-functions. We find that the pole of the kernel in TBA 
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equations (24) are located at 


1) 
t = 0.162117..., eg ee 3( 2 ) = 0, (29) 


on 1) 
t = 0.397459..., gyn 3( 722 ) — 0, (30) 
3 I 


which are nothing but the marginalstability walls of the (A2, A2) AD theory. 
Crossing the first wall at ¢2—¢1 = 3, one needs to pick up the contribution 
of the pole and modify the TBA equations: 


i = = mi 
log Yı, (0 = idı) =|m1 ıle? + Kx Lia — Ki * Liz = L12(0 = 3 En idı), 


: pen ps Tmi 
log Y12(0 =n ida) =|m,2\e° — Kay x Li + Kx Li = I11(0 + 3 = ida) 
(31) 
It is more convenient to introduce the new Y-function 
1 


PEE E 1 
Yi 2(0 — zi) : 


vy) = Yı (0) (1 T zs Yı 1(0 + my) 
’ 3 


¥{) (0) = Yı 2(0)(1 


14 mm ty 
Yi,2(0-—% Y1,1(0) 
1 ? 
Yı,1(0)Y1,2(0— 3) 


YOW = 


(32) 
such that the modified TBA equations can be rewritten as a three-TBA 
equations system 


A —(1 —(1 a —(1 
log Y{) (0 — igi) =|mi |e? + Ran. -Kız REG? + Ky 12 KLO: 


% —(1 —(1 = —(1 
log Ye — iġ2) =|my,2|e? + K x T. — K21* be — Kyj2* ine (33) 


$ —(1 —(1 —(1 
log YO (6 — idia) =|rmigle® + Kx LY) + Ki T - KhaN). 


These new Y-functions are identified to the WKB periods in the following 
way 


log Y (9) = eTl, > log Y} (0) = eels, , (0 aS 


(1) 0 (34) 
log Yio (9) =e hy, 3441.2 (9), 


which have been tested numerically. Moreover, we can cross the second wall 
at ġ2 — ¢, = 27/3 and arrive at the maximal chamber, where four-TBA 
equations system of vo : eu Ye ) and y ) were found. At the monomial 


potential, i.e. p(x,t = 1) = x? — 8, we can compute the classical periods 
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and the masses of the Y-function, which leads to a symmetry between the 
Y-functions 


} . f i Ti 
log ¥{9 (0 — idı) = log Y{ (0 — ids — mi) = log YP (0 — id — z) (35) 


The TBA equations at the monomial potential thus become a two-TBA 
equations system, which has the same form as the (symmetric-) D4-type 
TBA equations !1}, 


5. Conclusions and discussion 


In this paper, we have found the correspondence between the WKB periods 
of the third order ODE and the Y-functions satisfying the TBA equation 
in the quantum integrable model. We have studied the wall-crossing of the 
TBA equations on the moduli space. For the (A2, A2)-type ODE, we found 
the TBA equations at the monomial point become the D-type TBA equa- 
tions. Our observation can be regarded as the duality between (Ag, A2)- 
type and D4-type AD theory at the quantum level. The TBA equations 
of different integrable models are connected through the wall-crossing. It 
would be interesting to see how these integrable models are unified at the 
IM side. 
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In this note, we discuss the infrared physics of 5d M = 1 Yang-Mills theories 
compactified on S!, with a view toward 4d and 5d limits. The Coulomb phase 
boundaries in the decompactification limit are given particular attention and 
related to how the wall-crossings by 5d BPS particles turn off. On the other 
hand, the elliptic genera of magnetic BPS strings do wall-cross and retain the 
memory of 4d wall-crossings, which we review with the example of dP2 theory. 


Keywords: Gauge theory; Wall-crossing. 


1. Introduction 


When one realizes 5d M = 1 theories by geometric engineering as M-theory 
on a local Calabi-Yau? ®, BPS objects are realized by M2 and M5 branes 
wrapping 2-cycles and 4-cycles respectively. The former gives electrically 
charged particles, including dyonic instantons, while the latter gives mag- 
netic strings. Compared to their 4d counterpart ™S, these 5d theories look 
very simple; the 5d prepotential is at most a piecewise cubic function of the 
real Coulombic vev’s, while on S', one must deal with the special Kahler 
geometry of complex vacuum expectation values. 

The simplicity of 5d M = 1 theories is gratifying but at times appears 
too simple in that the 5d theory and the same theory compactified on a 
circle St seem superficially very disparate. The latter acquires a much richer 
character. This is partly because the compactification produces particle- 
like monopoles from magnetic strings wrapping S! and allows wall-crossing 
in the 4d sense. In the compactified theory, the wall-crossing involve BPS 
particles carrying several kinds of charges and could be very complicated. 
But on the other hand, in the 5d limit where Rs — ov, the ubiquitous 
wall-crossing phenomena of 4d‘ turns off in 5d, as far as particle-like BPS 
states are concerned. 

In this note, we wish to study 5d rank one theories compactified on a 
circle to study how the decompactification limit emerges in the large radius 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
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limit and how the complicated wall-crossing turns off and characterizes the 
co-dimension-one boundary of the Coulomb phase. We will see shortly the 
two are closely related to each other, and one recovers the known fact that 
5d Coulomb phases end where a BPS magnetic string becomes tension- 
less. For some theories, this tensionless limit coincides with the symmetry 
restoration of the non-Abelian gauge symmetry, but not always. 

On the other hand, the absence of wall-crossing really addresses the 
stability of 5d BPS particle states, whose central charges are all “real.” 4d 
monopoles uplift to monopole strings, and the central charge density thereof 
is by and large imaginary. In fact, the Kaluza-Klein particles also come with 
an imaginary central charge. This explains why, upon S! compactification, 
the wall-crossing turns on immediately. For a finite size of this circle St, 
the magnetic strings remain extended, so the relevant counting should be 
given by the elliptic genus!®!9. In the second part of this note, we will take 
a close look at the wall-crossing of this elliptic genus across a point where a 
charged matter becomes massless, or equivalently, across a flop transition. 

This note is a short version of Ref! and is organized as follows. In 
Section 2, we will revisit the question of the wall-crossing starting from 
the circle-compactified theory and addressing what to expect from the 
decompactification limit. Exactly what we mean by the absence of 5d 
wall-crossing is mulled over, from which we draw the anticipation that the 
boundary of the 5d Coulomb phase should be universally characterized by 
tensionless BPS strings. Section 3 will take the simplest types of rank one 
non-Abelian theories on a circle and investigate the Coulomb phase detail. 
Here we will see how, in the decompactification limit, the 5d Coulomb phase 
boundary is realized to make sure that certain marginal walls collapse, as 
suggested by the general discussion on wall-crossings. The wall-crossing of 
magnetic BPS string across such a flop transition is the topic for Section 4, 
where we illustrate using the example of dP2 near a massless quark point. 


2. Wall-Crossing or Not 


When it comes to BPS spectra, a very distinctive feature of the 5d the- 
ory is the absence of wall-crossing for BPS particles. Let us see how this 
happens and what this absence means. From the low energy dynamics of 
BPS objects, the usual 4d wall-crossing occurs due to runaway Coulombic 


directions emerging at special values of FI constants??? 


, which is in turn 
related to how the central charge phases are complex and phases of a pair 


can align at a co-dimension-one wall. With mutual intersection number, i.e., 
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with nonzero Schwinger product of charges, this leads to a wall-crossing. 

On the other hand, the 5d BPS particles would be all represented by 
M2 wrapped on 2-cycles so that a pair of 2-cycles in a Calabi-Yau 3-fold 
cannot have a mutual intersection number. While D4 on a 4-cycle can have 
an intersection number against D2 on a 2-cycle, the former is an extended 
object in the form of the BPS string since it is really M5 wrapping a 4-cycle. 
In fact, most BPS objects that would have entered the wall-crossing quiver 
gauge quantum mechanics in the 4d limit are made up of magnetic strings. 

Once compactified on a circle S! of sufficiently small radius Rs, however, 
the low energy dynamics are again well captured by BPS quiver quantum 
mechanics for D4-D2-D0 bound states, for which the wall-crossing are nu- 
merous. We will ask exactly how the wall-crossing turns itself off back in the 
limit R5 —> co. One quick answer to this is that since all magnetic objects 
become a string, as pointed out already, the question becomes moot if we 
ask for wall-crossing among 5d BPS particles. However, we can ask a little 
more by keeping track of BPS strings wrapped on S! considered as BPS 
particles. As we will see below, approaching the decompactification limit 
from the compactified theory, the disappearance of wall-crossing occurs for 
multiple reasons, different for different BPS objects. 

The usual 4d wall-crossing, where the complex central charge Z of the 
bound state aligns with Z; and Z; and decays as 


Z > lit 4; , (1) 


is possible only if D-branes that constitute particles į and j share a net 
number of open strings attached to them. This is in turn counted by the 
intersection number between the two cycles wrapped by the D-branes. In 
the R+! gauge theory viewpoint, on the other hand, the intersection num- 
ber translates to the Schwinger product of electromagnetic charges, 


ligj — qji FO, (2) 


so at least one of the two constituents must carry a magnetic charge. As 
such, one can imagine three logical possibilities for the 4d wall-crossing: 


e Z is non-magnetic, while both constituent Z; and Z; are magnetic, 
e Z and one of the two constituents, say Z;, are magnetic, 
e all three are magnetic. 


Also, the marginal stability wall of such decay, found by equating phases 
of Z; to that of Zj, would extend between Z; = 0 and Z; = 0 locus in the 
Seiberg-Witten moduli space. 
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For the first case where Z; and Z; are magnetic, since the magnetic 
central charge can be approximated as: 


Zmagnetic x tR5Tmonopole (3) 


which is purely imaginary and Tmonopole is the tension of the monopole 
string. In the 5d limit it should remain finite otherwise Z; and Z; will be- 
come too massive for a non-magnetic particle Z to decay to. That indicates 
the tension of the monopole string Tmagnetic Should go to zero at the same 
time when R5 — oo. If we take the absence of such wall-crossings as given, 
this in turn implies that no magnetically charged BPS strings should be- 
come tensionless, perhaps except at the boundary of the Coulombic phase. 
This is quite natural since from the (p,q) web realization of IIB theory, 
where the tensionless limit of magnetic strings translates to a collapse of 
a face, so the above claim that the Coulomb phase ends where a magnetic 
BPS becomes tensionless is quite natural already. 

For the second and the third, Z would represent a magnetically charged 
BPS string wrapped on St. The second type of marginal stability wall 
would emanate from a point in the moduli space where some components 
of charged BPS particles become massless. For instance, imagine a quark 
hypermultiplet with mass jr in the defining representation of the gauge 
group. Since for large Rs, the central charge of Z; is dominated by Zmagnetic 
which is purely imaginary, the wall emanating from the massless quark point 
following 


Arg(Z;Z*) = 0, (4) 


would initially extend along the circular Wilson-line direction of the small 
period 1/R5. Taking the Rs — oo limit, this means that there is a sense 
of discontinuity for magnetically charged strings across @ = uf, where ¢ is 
the 5d moduli on the Coulomb phase. 

Note that ¢ > py is precisely where an Sp(1) doublet fermion would 
contribute a Jackiw-Rebbi zero-mode”® to the monopole. With ¢ < pf, this 
zero-mode is lifted, so the d = 1 + 1 low energy dynamics of the monopole 
(dyon) strings change qualitatively across this point. This discontinuity 
translates to the conventional wall-crossing once the BPS string wraps a 
circle, S', which may be captured via the elliptic genus of such a magnetic 
string; later, we will review this with a concrete case of a magnetic string 
in the dP, theory, or Sp(1) theory with a single massive flavor. 

Marginal stability walls of the third type where magnetic strings decay 
into a pair of magnetic strings deserve different considerations. In 4d, such 
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decays involving three types of dyons are known, and these walls actually 
extend into the weak coupling region ?’. 
crossings, two independent adjoint scalar fields spanning the Coulombic 
moduli space are essential; The decaying states may be visualized by an 
analog of planar (p,q) string web, which cannot be drawn when the adjoint 
scalar is real?+. Such walls would not survive the decompactification limit 
since half of the Coulombic directions, corresponding to the Wilson line vev, 
collapses due to Rs —> oo. Thus, such 3rd type of marginal stability would 
also turn off in the decompactification limit, leaving behind the simplest 
possible chamber. 


However, for this type of wall- 


3. Sp(1) Theories on a Large St 


In this section we will analyse the pure 5d Sp(1) theory on a large S! and 
illustrate how does the second type of marginal stability walls behave in 
the strict 5d limit. In particular, we will restrict to the Sp(1) theory with 
zero theta angle, namely the F0 theory, for a general discussion involving 
Sp(1) theory please refer to Ref?. 

The Coulomb phase of 5d M = 1 supersymmetric gauge theory is fa- 
mously described by the Intriligator-Morrison-Seiberg (IMS) prepotential®, 
which is one-loop exact. In particular, the prepotential for Sp(1) theory is: 


Fims(¢) = pod? + fys ; (5) 


where ¢ is the real Coulmob phase moduli and po = 877/g? is the 5d 
inverse coupling-squared and also the instanton mass. The first derivative 
gives the monopole string tension Tmono as, 


2r Ob 27 , 


tT mono = 


where we used 27inono on the left hand side as a reminder that the 4d 
monopole central charge upon a compactification is imaginary in the large 
radius limit. The second derivative produces the pure imaginary 5d cou- 
pling as 


_ 1 0? Fims _ i( uo + 4¢) 
a a) a a K 


with Tsa = 1877/9? off- 
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After compactification on St, the prepotential will receive instanton 


corrections, which becomes®!! 14: 


0° Fsa,s1 (a, uo) 
0a? 


=8+ Ninn 


e 4nmaks e72nn(2a+pH0) Rs 
1 
m,n>0 


——EE 3 
_ ee (—2m == 2n) ’ (8) 


where a is the complexified version of the 5d moduli ¢ combining with the 
Wilson line around the circle, and in the asymptotic region of the Coulombic 
moduli space one has Re(a) ~ $. Nm,n is an invariant which counts the 
number of the holomorphic 2-cycles in Calabi- Yau picture, where the exact 
numbers can be found in Ref.'+. We are mostly interested in the interface 
between 4d and 5d, so we take Rs to be sufficiently large in the following 
discussion whereby |g?/27Rs| < 1. 

Let’s assume the bare coupling square is positive in the following dis- 
cussion. Since we are working with large Rs (uo Rs >> 1) , one can set n = 0 
in (8) to obtain: 


OF sq gi e 47 aRs 
“aa = 8+ 16 Soe — ar} 3 (9) 
where we use the fact that the only non-zero Nm,o is Ni o = —2 for FO TE 


The effective 4d coupling for the compactified theory is obtained by inte- 
gration: 
2 4 sinh? (2ra Rs) 
= 27k, x i— log | ——— 10 
Tad TLt5T5d t 0g ( (QRsAqcp)? )’ (10) 


OF. 4 : P 
where Tsa = +e is defined parallel to (7) and the integration con- 


stant is determined by comparing (7). Here Agop is the effective QCD 
scale for the compactified theory defined as 27R5AQcp = e~THoRs/2 such 
that if Re(a) < 1/27Rs5 and the 4d approximation can be trusted, (10) 


indicates: 
1 1 2a)? 
z aS log Co , (11) 
J94 eff An QCD 


which is the correct 4d coupling implied by 8-function. 
The moduli space for FO theory with positive uo is depicted as figure 
1*. There is a strongly coupled region ¢ ~ Agcp near the tip of the cigar, 


a Actually for a generic Sp(1) theory with matters, the moduli space is a double-copy 
of figure.1 and there are two distinct holonomy saddles. We refer Ref. 15 for a detailed 
discussion of the moduli space. 
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and the theory is effectively the strongly coupled pure Sp(1) Seiberg-Witten 
theory. Around Agcp there are two singularities where the mass of (0, 1) 
monopole and (2,—1) dyon will separately become zero, and there is a 
marginal stability wall connecting these two singularities as shown in the 
figure 1. 

Now let’s see what happens in the decompactification limit R5 — oo. 
Since the instanton contribution vanishes in this limit, we can identify 
Re(a) = ¢. The period of a isa ~ a+ i/2Rs , therefore as Rs — oo 
the imaginary part of a is zero such that a = ¢ in this limit, and the cylin- 
der will become the half-line in figure 1. Also Agcp becomes much smaller 
than 1/27Rs5 in this limit therefore the two singularities and the marginal 
stability wall both shrink to the endpoint ¢ = 0 in the 5d moduli space. At 
that point both W-bosons and monopoles become massless and the Sp(1) 
gauge symmetry is restored. 


S! x Ris 


? 


Agep a 1/Rs 


Symmetry restored 
o=0 


Fig. 1. Moduli space for F0-theory with positive coupling-squared. In the large Rs 
limit, entire strong coupling region is pushed to the far left of 1/Rs, which itself is 
pushed toward ¢@ = 0. No trace of the marginal stability wall is left behind. 


Ria 


’ 


We have discussed the first kind of marginal stability wall where the W- 
boson decays into a pair of (0,1) monopole and (2, —1) dyon and seen how 
the wall collapses to the endpoint of the 5d moduli space. More generally, 
for any pair of magnetic particles Z; and Z;, their central charge are of 
order ~ Agcp in the core region or dominated by the magnetic part: 


Zmag ~ iR52a(2a + uo), (12) 


in the asymptotic region Re(a) > 1/27 Rs. 
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On the other hand, since Z is non-magnetic, its central charge is a 
combination of the electric, instanton and KK parts: 
Zelec = 2a, Zinst = Ho, ZKK = LA (13) 
Rs 
The wall, if any, must remain in the region where Zelec, Zmst; ZKK are 
comparable with Zmag, end where the magnetic particles Z; or Zj become 
massless. Zelec ~“ Zmag requires a ~ Agcp while Zinst ~ Zmag can happen 
only if Re(a) ~ 1/27R5. ZK scales down as 1/Rs also. Therefore, no walls 
of the first kind can extend beyond a ~ 1/27Rs5, which in turn collapse to 
ġ = 0 in the decompactification limit. 


4. Wall-Crossing of a Magnetic BPS String 


In this section we will explore the second type of marginal stability wall, 
where the bound state Z and one of its component Z; are magnetic. As 
discussed in the previous section, one may consider adding a flavour with 
positive mass pf to the 5d theory, and on the Coulomb branch it becomes 
massless at @ = uf. From the M-theory picture, the quark is represented by 
a M2 brane wrapping on a compact 2-cycle, and at the massless point in the 
Coulomb branch the volume of the 2-cycle vanishes, which is responsible 
for a flop transition of the local Calabi-Yau. 

The correct description of these magnetic objects is, in fact, well known 
to be governed by (0, 4) nonlinear sigma models which has been well-studied 
in various literatures!®1’. The elliptic genus that counts the BPS states 
for such a theory is defined as1819; 


Z(T,7, y) = Tens F2(—1)F en Dasta qto- glo Gt ePrives”, (14) 


where the monopole string is represented by the M5-brane wrapping on the 
compact divisor S in the local Calabi-Yau X, q* are membrane charges, 
Day is the intersection matrix between 2/4-cycles and s“/2 is a half-integer 
charge shift due to the Freed-Witten anomaly ?5. 

The evaluation of the elliptic genus is based on a fact that the elliptic 
genus Z(7,7,y) is subject to a 6-function decomposition (see!?°for a re- 
view). Eventually the elliptic genus is totally determined by the geometry 
of the compact divisor S and the way it is embedded in the ambient Calabi- 
Yau X. The change of geometry during flop transition will induce a jump 
of the elliptic genus, and for dP2 geometry that is: 


O14 (T, ya) 


Zap, = -Zz 
? n(T) 


, (15) 
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where the 6-part is exactly the partition function of a chiral fermion gen- 
erated by the Jackiw-Rebbi zero-mode as discussed in section 2. 
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de Sitter Duality and Holographic 
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We perform the resummation of the infrared logarithms in the inflationary uni- 
verse. Applying the renormalization group, we derive the stochastic equations 
as the effective theory at the horizon. We focus on the conformal zero mode 
to respect local Lorentz symmetry. Under Gaussian approximation, we derive 
the fundamental equation for the Universe (EqU). We also derive the identical 
equation from the first law of thermodynamics in a dual geometric picture. 
We believe it is a convincing evidence for de Sitter duality between quantum 
stochastic physics on the boundary and classical thermodynamics in the bulk. 
The equation for the Universe (EqU) possesses the solution with the ultraviolet 
fixed point. It also contains the inflationary universe with the power potentials. 
We discuss possible scenarios for the very early universe with decreasing e. We 
argue inflationary universe subsequently dominates to maximize the entropy 
and e problem is naturally solved. 


1. Introduction 


dS space may be decomposed into the bulk and the boundary, i.e., the 
sub-horizon and horizon. From a holographic perspective, we consider the 
conformal zero mode dependence of the Einstein-Hilbert action: 


1 T 
-—— d‘z(Re’ — 6H’e*”) ~ — (1 — 4w? 1 
acy | VIAR -6R ~ (1-404), A 
where the gauge fixing sector is suppressed because it does not produce 
the IR logarithms in the background gauge. Our gauge fixing procedure 
and the propagators are explained in+. The semi-classical dS entropy was 
obtained by rotating the background spacetime d$* to 94 in (1)°. 


*E-mail address: kitazawa@post.kek.jp 
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The quadratic part of w constitutes a Gaussian distribution function for 
the conformal zero mode, 


pla) = [= exp (= $u). (2) 


It may represent an initial state of the Universe when the dS expansion 
begins. In order to describe the time evolution of the conformal factor of 
the Universe, we introduce a new parameter &(t). 


4€(t) ( _ AO), (3) 


p(t), w) = a 


€ is the only parameter in the Gaussian approximation. We work 
within the Gaussian approximation since it is an excellent approxima- 
tion for gravity with the small coupling g. The von Neumann entropy 
S = —tr(plogp) ~ 1/2 log(g/&) becomes larger as € becomes smaller. Thus 
the diffusion triggers an instability in de Sitter space. 

In terms of the distribution function, the n-point functions are defined 
as follows 


(w" (t)) boundary = J tootoo". (4) 
In particular, the two-point function of the conformal mode is given by 
2 g 
t))boundary = z7x- 5 
(w* (t) boundary O (5) 


The negative norm of the bulk conformal mode indicates that the p 
diffuses toward the future. In fact the perturbative quantum expectation 
of (1) gives g(t) ~ g(1 — 2yHt) ~ g(1 — 3gHt). The duality between 
the inflation and quantum gravity is also based on this one loop effect %”. 
However such an estimate is reliable only locally Ht < 1. In cosmology it 
is essential to resum all powers of IR logarithms (Ht)” to understand the 
global picture. For such a purpose, we find our holographic approach is up 
to the task. 

We investigate the dynamics of conformal mode after integrating the 
bulk mode. The two point function gives rise to IR logarithms. 


3g [^ dk 
<W“ >bulk = p 
Ha(t) 
39 ry SI Ne) (6) 
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We assume there is no time dependent UV contributions. We focus 
on the Hubble scale physics where a(t) = 1/ — TH = exp(Ht). We recall 
that —Tu ~ 1 at the Horizon. While the wave functions of the bulk mode 
oscillates, those of the boundary mode do not. That is why we call it the 
zero mode. Our strategy is to integrate out oscillating mode first. 

The finite bare distribution function is given by subtracting the bulk 
mode contribution above the Hubble scale. We thus construct low energy 
effective theory around the Hubble scale. Such a theory is holographic and 
finite after the subtraction. 


psg =exp(— ===). (7) 


The renormalization scale of the low energy effective action is the Hubble 
scale. As pg is independent of the renormalization scale, the renormalized 
distribution function obeys the following renormalization group equation. 


ae ae dae (8) 


where O denotes a derivative of O with respect to the cosmic time t. The 
factor 3g/4 in the diffusion term is the projection factor to the IR re- 
gion. The conformal mode w consists of a minimally coupled field X. 
w = V3X/4+ Y/4%8. We neglect Y because it has the effective mass 
msp = 2H?. There is no drift term in the reduced space consisting only 
of X. In fact X is nothing other than the curvature perturbation Ç. 

The gravitational FP equation (8) is obtained by integrating the quan- 
tum bulk modes inside the horizon. It turns out to be a diffusion equation 
due to the lack of the drift term. The solution is the Brownian motion as 
it is jolted by the horizon exiting modes. The FP equation is a dynam- 
ical renormalization group equation. We can sum up the IR logarithms 
log” a = (Ht)” by this equation to find a running coupling g(t). 

The FP (diffusion) equation shows that the solution is the Gaussian 
distribution with the standard deviation increasing linearly with the e- 
folding number N(t)°. 

“gilt = INg) = a (9) 

The standard deviation is related to N as 1/ = 6N in (9). It is con- 

sistent with a standard Brownian motion prediction. The von Neumann 
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entropy thus increases logarithmically, 


iL 1 1 
6S = log- = —1 N(t). 1 
S 3 8 e 5 log 6 (t) (10) 


Identifying the von Neumann entropy of conformal zero mode with the 
quantum correction to dS entropy, we obtain the bare action with the 
counter term 

1 1 1 
By requiring the bare action is independent of the renormalization scale: 
namely N , we obtain the one loop £ function. 


B= eal) = -3N (12) 


We find the running gravitational coupling as 
2 
N) = ——. 13 
g(N) loa) (13) 
The holographic investigation at the boundary shows that g is asymp- 
totically free toward the future!?!!. The renormalization group trajectory 
must reach Einstein gravity in the weak coupling limit for the consistency 
with general covariance!*®. We find that it approaches a flat spacetime in 
agreement with this requirement. 


2. Equation for the Universe and de Sitter Duality 


The Gaussian distribution of the conformal zero mode is characterized by 
the standard deviation 1/€. Although there is no inflaton in Einstein grav- 
ity, we propose to identify the inflaton f? as f? « 1/€. In our interpretation, 
the inflaton is not a fundamental field but a stochastic variable. The two 
point function at an equal time grows due to the Brownian motion: IR 
logarithmic fluctuations 1/§ ~ N. While the inflation theory is specified 
by the inflaton potential, the dynamics of quantum gravity is determined 
by the FP equation which describes the stochastic process at the horizon. 
We thus argue the classical solution of the inflation theory satisfies the FP 
equation as well. 

It is likely that there are multiple elements in the universality class of 
quantum gravity/inflation theory. The inflation era of the early Universe 
may be one of them. As we discuss shortly we find a pre-inflation era which 
is indispensable to launch inflation era which in turn necessary to trigger 
the big bang. 
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We evaluated the time evolution of entropy to the leading log order in 
(10). In order to take account of the higher loop corrections in g, the FP 
equation should be generalized. It turns out to be just necessary to intro- 
duce the dimensionless gravitational coupling g(t) and e-folding number: 
N = f dt'H(t’). 

2 
ow a sa ta a) 
Since the equation is local, we consider the following local solutions. 


E(t) Ae(t) a 
moO ga 


We next put the ansatz into the FP equation and find the condition for the 
background to satisfy. 


(15) 


— 2’ = Lg. 1 
ƏN 8€ 1 (e 
We obtain a remarkably simple equation 
ð |, ght) 
log = 6€. 1 
aay on 2 = 6 (17) 


(17) determines the evolution of von Neumann entropy S = 5 log # with 
respect to N. This formula shows the validity of our postulate that von 
Neumann entropy of conformal zero mode constitutes the quantum correc- 
tion to de Sitter entropy. We call it the equation for the Universe (EqU). 
The increase of the entropy S = 1/g can be evaluated by the first law 


TAS = AE where AE is the incoming energy flux of the inflaton??. In 
this way, the one of the Einstein’s equation is obtained: 

H(t) = —4nGy f?. (18) 
From this formula, we obtain 

2e = — ? log g(t) (19) 
We add the same quantity 1/ N to the both sides of the equation. 

1 o ~ 
2e + = = -l N). 2 
c+ = -gy los(gW) (20) 


For power potential inflationary mere (20) is rewritten as 


6€ = + tog! E 
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It is precisely our EqU (17). We have derived it in a dual geometric picture 
here in place of the original quantum stochastic picture. 
The equation (17) has inflationary solutions with power potentials. 
g=cN?, aa (21) 
12N 

Here c is an integration constant. m denotes the power of the potential: 

f™. It is convenient to replace N by Ñ where N = N.—N. Ne denotes the 

e-foldings at the end of inflation. (17) becomes as follows after the change 

of the variable, 


—~— log (Ñ) = 6€(N). (22) 


Although the dS entropy can be explained by quantum effects alone for 
the weakly coupled inflaton solution, the strongly coupled inflaton solution 
is a dual object in the sense that geometrical description is reliable. 

It proves the consistency of the EqU and the confirmation of de Sitter 
duality. Since (17) and (19) are equivalent, the solutions of the former 
satisfy the latter. 


3. Inflation and UV completion 


In the literature, dN formalism is widely used to investigate the curvature 
perturbation. It underscores the validity of the stochastic picture of the 
inflation. '? 16 Let us consider the fluctuation of the curvature perturbation 


¢. 


H 
= ô = — ó > 
Ç= ôN 7 f 
saro E=: (23) 


4r? 


We obtain in the super-horizon regime: 


e < (t) >= 2 < (Zep > 


ON ON f 
1 ð 2 

= MZ ON ~ OV te 
H? g 


= GME aH (24) 


We recall the following identity holds at the horizon exit t = t,. 
D)e = k. (25) 
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It is nothing but choosing our renormalization scale as logk = Ht. Let 
P~k":~!, The scaling dimension of P can be estimate as 


d ‘ 
k— (4 log p(t) — 2 log o(ts)) 


a (4 log p(tx) — 2 log J(ts)) 


dN, 
=2(n — 36), (26) 
where 
p 6 
C= 20) n=- +e. (27) 
2H? po 


In power inflation potential models V ~ f™, € and e scale as 1/N due to 
the universality of the random walk. They belong to the same universality 
class as the following scalar two point functions scale as 


I m 

Tia A (20l 

The scaling exponents agree with our FP equation and the ô N formalism. 
o 1 m+2 

l-—ns = l N(t)) = 2e + —~ = ——.. 29 

ns = gp 8UNO) = 2+ ey = Sa (29) 


There is a UV fixed point in our renormalization group. We study it 
next. FP equation (17) enables us to evaluate higher order corrections to 
the 8 function. The expansion parameter is 1/ log N. We can confirm that 
the following gf and ¿p satisfies (17), 

2 1 E 1 
fog N | TaN xl- meN 
Thus, the 6 function, e and the semi-classical entropy generation rate are 
given by 


Es = (30) 


gf = 


2 4 1 1 
= ~~ gf = — YH oes = —2(1 - 2) — —. 31 
j O log NI log’ N = log(NV)3 ( Tog N Tog? N ) 
10 1 
= z pf. 2 
€f = -3N 2807) = -Nh 32) 
1 1 
E EE EN 33) 


aN = Ny NA 


A remarkable feature is that the coupling has the maximum value g = 
1/2 at the beginning. It steadily decreases toward the future as the 6 
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function is negative in the whole region of time flow. It has two fixed 
points at the beginning and at the future of the Universe. The existence of 
the UV fixed point may indicate the consistency of quantum gravity. The 
single stone solves the € problem!” as well since it vanishes at the fixed 
point. The 8 function describes a scenario that our Universe started the 
dS expansion with a minimal entropy S = 2 while it has S = 10!?° now. It 
corresponds to N = e? ~ 7.4,a = e ~ 1.6 x 10°. The just born Universe 
is rather large which reflects the critical coupling g = 1/2 is rather small. 
In terms of the reduced Planck mass, H? /4r? M2 = 1. 

Since we work with the Gaussian approximation, our results on the UV 
fixed points are not water tight as the coupling is not weak. Nevertheless we 
find it remarkable that they support the idea that quantum gravity has a 
UV fixed point with a finite coupling. In fact 4 dimensional de Sitter space 
is constructed in the target space at the UV fixed point of 2+ € dimensional 
quantum gravity!8. 4 dimensional de Sitter space also appears at the UV 
fixed point of the exact renormalization group !??°. 

Such a theory might be a strongly interacting conformal field theory. 
However, it is not an ordinary field theory as the Hubble scale is Planck 
scale. Our dynamical ( function is closely related to the cosmological hori- 
zon and physics around it. The existence of the UV fixed point could solve 
the trans-Planckian physics problem. A consistent quantum gravity theory 
can be constructed under the assumption that there are no degrees of free- 
dom at trans-Planckian physics?!. In this sense, it is consistent with string 
theory and matrix models. The Universe might be governed by (30) in the 
beginning as it might be indispensable to construct the UV finite solutions 
of the FP equation. 


4. The Inflationary Universe and Pre-History 


The inflaton may be identified with the stochastic variable f whose corre- 
lators show characteristic features of Brownian motion. < f? >= N and 
gx N™/2 =< f™ >. 6€ = 1 — n, measures the extra tilt of the scalar two 
point function k"*—! with respect to k. 

The left-hand side of (17) can be identified with 1 — ns where ns is the 
scalar spectral index. Let us recall that 1 — ns is expressed by the slow-roll 
parameters € and 7, 


1— ns = be — 27. (34) 
The equation (17) is equivalent to (34) for the f™ inflaton potential where 


n = (m—1)/(2N). 
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We thus conclude: 


i 10 m 
g=cN?, € 5 an 099) rr 
m+2 
1—n, = 6€ = ——. 35 
s 2N (33) 


In this case, the concave power solutions can be obtained by formally re- 
placing m by 1/n. We admit it is not entirely clear why concave potentials 
are relevant. It is possible that the convex potentials are already excluded 
by observations. 

In contrast, there is more room for concave potentials to accommodate 
the observational information with the judicious choice of n. 


ô | (a) 1 

=—--=—lo =— 

E 2 ƏN g\g InN’ 
1. 
=+2 

ee ns (36) 
2N 


In our view the inflaton behaves as an particle in the Brownian motion. Its 
trajectory is made of infra-red fluctuations. Nevertheless the convex and the 
concave potentials in the inflation models seems to form two distinct groups. 


The concave potentials are promising avenue to explore right now ?’. 


Table 1. 


4N 8Ñ 12Ñ 
1 3 2) 
M 
a 2Ñ 4Ñ 6Ñ 


In Table 1, we list expected € and 1 — n, in the power potential model 
with f”, n = 1,2,3. We note 1 — n, is bounded from below by 1/N 
while r = 16e is not. It is consistent with the current observations. It is 
important to establish the bound on n. We argue the order of magnitude 
estimate can be trusted for the concave potentials as they are in the weak 
coupling regime. The semi-classical de Sitter entropy favors smaller n as 
the corresponding entropy grows faster: S ~ (1/ Ñ). 

Although it is still highly speculative, we argue that our Universe is 
likely to be located close to the point m = n = 1. Since the convex 
potentials are excluded by observation, it must be a concave potential with 
a small n. The convex potentials correspond to strongly coupled systems. 
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Here the naturalness and the anthropic principle may come in. It is hard 
to believe that our Universe comes out of a strongly coupled system. It 
is much easier to accommodate hierarchies in the weakly coupled system. 
Our Universe might to be located near the boundary between the stable 
and unstable universes like the standard model of particle physics. Since 
the Universe can stay at the fixed point forever, it is likely to be not far 
from it. 

Our FP equation is expected to go beyond the resummation of leading 
IR logarithms which results in the logarithmic decay of g. It is still a great 
surprise to find inflation theory as its solutions. Fortunately EqU can be 
derived in a dual geometrical picture. This fact gives a non-perturbative 
evidence for the dS duality. 

The solution (30) is UV complete. g is attracted to the fixed point as 
we roll back the history of the Universe. However, it cannot terminate the 
eternal inflation as € = —(1/2)0logg/ON ~ 1/(2N log N) decreases with 
time. On the other hand, the solution (21) is not UV complete but it can 
end the inflation as € ~ m/4Ñ increases with time. These solutions generate 
the entropy in different ways: We consider the leading de Sitter entropy 
1/g(t). The t dependence of g(t) depends on the quantum corrections 
through the EqU equation. 1/g ~ log N for the former and 1/g ~ 1/Ñ = 
for the latter. From the perspective of the dominant entropy principle, (30) 
is chosen initially and (21) is chosen after log N ~ 1/N2*. That is to say, 
gi in (30) describes the newly born Universe and gə in (21) describes the 
inflation era. 
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In recent years, the extension of the notion of symmetry using the picture of 
topological defects and its applications have been actively studied. One di- 
rection is the so-called non-invertible symmetry, in which non-invertible topo- 
logical defects are treated as symmetry. We have constructed non-invertible 


topological defects from duality in 4D lattice pure Zz gauge theory. This work 


is in collaboration with Y.Nagoya and S. Yamaguchi. 1 


Keywords: Non-invertible; symmetry; topological defect. 


1. Introduction 


Symmetry is one of the important tools in the non-perturbative analysis 
of quantum field theory. In recent years, the extension of the notion of 
symmetry and its applications have been actively studied”. In extended 
symmetry, topological operators are considered as symmetries. There are 
also several types of extended symmetries, depending on the properties of 
their topological operators. The one direction is called the non-invertible 
symmetry. The unitary topological operator corresponding to the ordinary 
symmetry has an inverse transformation due to the group structure of the 
symmetry. There are general topological operators that do not have such 
inversions, and these types of topological operators are treated as non- 
invertible symmetries. 

The most famous concrete example of non-invertible symmetry is the 
duality defect associated with the Kramers-Wannie duality in the 2- 
dimensional Ising model. Since this duality defect vanishes when acting 
on a single spin operator, it is a non-invertible for which there is no inverse 
transformation. It is also known that there is a non-invertible symmetry 
defect called Verlinde line in the 2-dimensional RCFTs. Thus, while exam- 
ples of non-invertible symmetry in 2-dimensions were known, examples of 
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non-invertible symmetry in 4-dimensions were less understood than those 
in 2-dimensions. 

One approach to constructing topological defects in lattices is the work 
done by Aasen, Mong, and Fendley(AMF).*:+ In this paper, We explain 
how to use this AMF method to construct a high dimensional topological 
defect. + 


2. 4-dimensional Zz lattice gauge theory 


In this paper, we examine the 4-dimensional pure Zs lattice gauge theory. 
We prepare the 4-dimensional cubic lattice and assign a link variable Um = 
+1 to each link m. The partition function of this theory is given by, 


=e | KX [[ Un]. (1) 


{U} iEP mE 


Here P is the set of all plaquettes, K is a positive constant parameter, and 
i is the set of four links included in the plaquette i. It was shown by 
Wegner that there is a duality in this theory, called the Kramers-Wannier- 
Wegner duality(KWW duality). According to the KWW duality, these 
theories with the parameter K and K are equivalent when they satisfy the 
relation 


sinh 2K sinh2Kk = 1. (2) 
The situation K = K = K. is called self-dual point. 
1 
Ke = —5log(-1+ V2). (3) 


Next, I explain how to extend the AMF method to work in the 4- 
dimensional pure Zə lattice gauge theory. To consider the duality de- 
fect, we prepare two kinds of lattices as shown in Fig.1. To be more 
specific, we introduce coordinates (x1, 22,73,%4) of R4. On this Rt, we 
define two lattices A := {(x£1, £2, £3, £4)|£1, £2, £3, £4 E€ 2Z} and Â := 
{(£1, £2, £3, £4)|£1, £2, £3, £4 E 2Z + 1} that are dual to each other. In 
each lattice, the line segment consisting of the two nearest points is called 
a link, and the smallest square consisting of the links is called a placket. A 
line segment connecting points on different lattices is not called a link. 

We assign link variables Um = +1 to the links on A, but constant 
weights rather than variables to the links on the dual lattice A For this 
reason, lattice A given a variable is called an active lattice, and lattice A 
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Fig. 1. Schematic figure of a lattice. In this figure, the lattice is depicted as 2- 
dimensional, but the actual lattice treated in this paper is a 4-dimensional one. The 
black lattice represents the lattice A, and the blue lattice represents the lattice A. They 
are duals of each other. 


given a constant weight is called an inactive lattice. This inactive lattice is 
an auxiliary lattice prepared to construct the duality defect. 

Consider a 16-cell as the basic unit on these lattices as shown in Fig.2. 
A 16-cell consists of 16 tetrahedra such that each has one active link and 


one inactive link as edges. The surface of the 16-cell is homeomorphic to 
SS, 


Fig. 2. 3-dimensional stereoscopic projection of a 16-cell. The black plaquette represent 
an active plaquette and the blue plaquette represent an inactive plaquette. 


Since there is a one-to-one correspondence between active plaquettes 
and 16-cells, we assign Boltzmann weights to these 16-cells. If the four link 
variables in a 16-cell are a; = 0,1, (i = 1,2,3,4), we define the Boltzmann 
weight assigned to this 16-cell as 


W (a1, a2, a3, a4) = exp(K(—1)(1ta2tsta4)), (4) 


78 


This definition is equivalent to the definition of Boltzmann weights in 
Eq. (1) ,if the link variables for the link m are identified as Um = (—1)*™. 
We also assign a constant weight to each site or link, in addition to the link 
variable. We assign the weights for active links, active sites, inactive links, 
and inactive sites as s, l, 5, and l, respectively. These weights are deter- 
mined later to make the duality defect topological. The partition function 
is 


z=] TE«l/ Tell Dell We 


{a} active active inactive inactive 
sites links sites links 
[Wano ana aaa) (8) 
iEC 


where C is the set of all 16-cells, and jı(i), j2(i), j3(i), and j4(i) are the 
four active links in the 16-cell 7. aj is the link variable assigned to the active 
link j. 


3. Duality defect 


We discuss topological KKW duality defects in the Ze lattice gauge theory 
following AMF method?*. This duality defects are 3-dimensional opera- 
tors located at the boundary between two regions. The active lattice and 
the inactive lattice are swapped across the duality defect as shown in Fig.3. 
Because of this property, the building block of the duality defect is a tetra- 
hedral prism with a doubled tetrahedron (see Fig.4). Each building block of 
duality defect contains two active links. Let a, à = 0,1 be the link variables 
assigned to these two active links. We assign a weight D(a, à) (a, @ = 0, 1)to 
this building block. 

The entire duality defect is constructed by connecting these building 
blocks together. This defect is a 4-dimensional object in the lattice, but 
when the lattice spacing is set to zero, it becomes a 3-dimensional operator. 

Since the duality transformation just changes the description of the 
theory and does not change the observables, we expect that such topological 
duality defects exist. In order to make the duality defect topological, we will 
impose a commutation relation on the defect and find its solution. From 
now on, we will consider the commutation relation on a single 16-cell. On 
the surface of a 16-cell, there are 16 tetrahedra. Consider an arrangement 
A of a certain duality defect. We assume that part of A occupies part of 
the surface of the 16-cell we are focusing on. We also consider a deformed 
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Fig. 3. Schematic figure of the duality defect. A duality defect is placed at the boundary 
between the two regions. The active lattice (black dots) and the inactive lattice (blue 
dots) are swapped across the duality defect. The building block of the duality defect is 
a tetrahedral prism, depicted as a green parallelogram in this figure. 


Fig. 4. The building block of duality defects. the building block is a tetrahedral prism 
with a doubled tetrahedron. The 3-dimensional surface is composed of tetrahedrons each 
of which includes an active link and an inactive link. 


configuration B in which the tetrahedrons which are not filled by A are filled 
and vice versa. Outside of a focused 16-cell, the configurations of A and B 
are assumed to be same. We assume that if the topologies of configuration 
A and B are the same, then the weights of those two configurations are 
equal. We call this relation commutation relation(Fig.5). 

The commutation relation does not hold for some configurations, such 
as when the topology changes due to deformation. Therefore, we impose 
the following conditions on the commutation relation. 


e The filled tetrahedrons on the 16-cell of A and B are both non- 
empty sets. 

e The configurations A and B restricted on the surface of the focused 
16-cell are both simply-connected. 

e There is no connection such that the duality defects are connected 
only by sites or links. 


In the following, we prepare some definitions to express the commutation 
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Fig. 5. Schematic figure of a commutation relation. The circle represents a 16-cell. A 
green line represents a duality defect. The commutation relation means that the value of 
the duality defect does not change when the topology is transformed without changing 
the topology. 


relation as a mathematical expression. First, consider a 16-cell where the 
duality defect is not placed. Let m = 1, 2,3,4 be the four active links on the 
focused 16-cell, and ñ = 1,2,3,4 be the four inactive links. Let the set of 
links be M = {1,2,3,4} and N= {1, 2,3, 4}, respectively. The active link 
and inactive link pair (m, ñ) corresponds one-to-one to the tetrahedron. 
The set of tetrahedra on the surface of a 16-cell U can be represented as 
follows 


U = {(m,fi)|\m = 1,2,3,4, ñ = 1, 2,3, 4}. (6) 


Next, consider a 16-cell with a duality defect in configuration A. Suppose 
that I C U is occupied by a duality defect. The duality defect is defined on 
a tetrahedral prism doubling the tetrahedron of the 16-cell surface, so an 
additional active link appears. Since this additional active link corresponds 
one-to-one to the original inactive link, we use ñ as the label. And let 
E = {ñ|(m, ñ) € I} be the set of such additional active links. Thus, the 
weight of the duality defect on the focused 16-cell in coordination A can be 
expressed as 


W(a1,42,03,04) || D(am, an). (7) 
(m,ñ)EI 


In configuration B, the defect is placed at J = U \ I. Therefore, the weight 
of the defect in configuration B becomes 


W(aj,45, 43,45) [[ D(am,aa). (8) 
(m,ñ)EI 
Note that the Boltzmann weight variables are different for configurations 


A and B because the theories on the 16-cells are dual. We also consider 
E = {m| (m, ñ) € I}. 
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In this case, the commutation relation is written as follows 


5 W (a1, a2, a3, a4) 8% 191 3% 7 II D(am, Gn) 
M\E (mayer 
= X W(G;, ay, Gg, àa) 13e || Dlam,aa)- (9) 
N\E (m,njet 
Here the number of active sites, active links, inactive sites, and inactive 
links in the configuration A are denoted by a1, 61, &1, and 81, respectively. 
The number of active sites, active links, inactive sites, and inactive links 
on the configuration B are denoted by ag, 82, &2, and 2, respectively. The 
sums )/y7\p) 2o g g are defined as follows 


S-I oo YY. w 


M\E m€M\E am=0,1 N\E nEN\E 44 =0,1 


In a physically sensible solution, these values satisfy 


D(a,a) 40, 1s,l,3>0, KER, KO. (11) 
By solving for the commutation relations, we obtained the following unique 
solution 

D(a, a) =(-1)", (12) 

1 1 
Se BS 13 
5 fa (13) 
l=1, 3=1, (14) 

1 

K = K. = -3 loøg(—-1 + v2), (15) 
W (a1, a2, a3, a4) =exp(K(—1)(te2ta8444)), (16) 


Thus, a topological duality defect was obtained. Note that the value of K 
is determined to be the critical value Ke. 

Finally, we show that this topological defect is non-invertible. To do 
this, consider the situation as J = U as the placement of the duality defect 
in configuration A (see Fig.6). In this case, the topology of the duality 
defect changes with the deformation, so the commutation relation is not 
always satisfied. In fact, the calculation results in the following 


DD W (a1, a2, a3, a4) 8°18 38ST II D(am,ãñ) 


aı1,a2,a3,a4=0,1 (m, ñ)EU 


W (àz, àz, 43,4;)s*I*5*1*. (17) 


Sl- 
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For asymmetry defect placed on a closed manifold with no operator inserted 
inside, the weights are identical to the empty configuration”. However, for 
this duality defect, their ratio is 1/ V2. So, this defect is non-invertible. 


Fig. 6. Schematic figure of Eq.(17). The circles represent the 16-cell. The green line 
represents the duality defect. 


4. summary 


We construct topological KWW duality defects in the Zə lattice gauge 
theory following AMF method**. And this duality defect is non-invertible. 
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We show that the split /non-split transition is, with some exceptions, a conifold 
transition at a Do, (k > 2) or an Ey codimension-two singularity, and these 
non-split models have local deformed conifolds. These clarify, in these cases, 
that the origin of non-local matter generation and that the deformations of 
conifold singularities correspond to diagram automorphisms of the expected 
simply-laced Dynkin diagrams in the split sides. 


Keywords: F-theory, matter, singularity, resolution and deformation 


1. Introduction 


F-theory'? is a geometrical framework of nonperturbative compactifica- 
tions* of type IIB theory with general 7-branes. In six- or lower-dimensional 
F-theories, if a fiber type involves the condition that an exceptional curve 
splits into two irreducible ones, then the fiber type has two types, “split” or 
“non-split” , depending on whether the exceptional curve can split globally 
or not®. In the split models, each intersection diagram of exceptional curves 
that arises after the resolution corresponds to the expected ADE Dynkin 
diagram implied by Kodaira’s classification’(Table 2). In the non-split 
models, since the two split exceptional curves are identified by monodromy, 
each ADE gauge symmetry in the split models is reduced to the correspond- 
ing non-simply-laced gauge symmetry by being subject to a projection by a 
diagram automorphism. The fiber type In (n = 3,4,...), Iž (n =0,1,...), 
IV or IV* can involve such identification of exceptional curves ®®. 
Non-split models have some puzzles which require a new understanding 
of the non-local generation of charged matter. In split models, if the mat- 
ter fields locally exist at all codimension-two singularities’, the number of 
the matter fields matches the anomaly-free condition®. This is one of the 
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aIn F-theory, strong- and weak-coupling regions coexist on the base of the total space 4&5. 
bIn F-theory, there are 7-branes on codimension-one loci (in the base: codimension-two 
in the total space) where the elliptic fibers become singular. These codimension-one loci 
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reasons why the massless matter fields are localized at all intersections of 
7-branes in F-theory. However, in non-split models, there is a puzzle that 
the anomaly-free spectrum and the naive counting are not consistent °!? 1°, 

Moreover, at a Do, (k = 2,3,...) or an Ey codimension-two singularity, 
some conifold singularities® remain after all codimension-one singularities 
are blown up in the split models, but not in the non-split models. In the 
split models, since we can yield new two-cycles by small resolutions of the 


9.17 we can obtain an intersection diagram of excep- 


conifold singularities 
tional curves which is different from one on a codimension-one singularity 
and explains the enhancement of the symmetry. In the M-theory dual, an 
M2-brane wrapping a new two-cycle generates local matter multiplets 18. 
On the other hand, in the non-split models, no additional blow-up is re- 
quired at these singularities, so the intersection diagrams remain the same 
and there is no new two-cycle where an M2-brane can wrap around!?. 


12,19: we discuss these puzzles 1? 


This paper is based on our recent papers 
and show the split /non-split transition is a conifold transition at a D2, or an 
E; codimension-two singularity with some exceptions (I3,,_,)'°. Since this 
implies these non-split models have local deformed conifolds 1°, we can yield 
new three-cycles which are non-local in terms of an elliptic fibration. These 
clarify, in these cases, the “non-local” matter fields may satisfy the anomaly- 
free condition and the deformations of conifold singularities correspond to 


diagram automorphisms of the ADE Dynkin diagrams in the split sides. 
2. Puzzles on matter generation in non-split models 


In this section, we discuss non-split models using a non-split Jg model as 
an example!?. Specifically, we examine their matter spectra and the result 
of the resolution at a Dg singularity of a non-split Jẹ local equation. Thus, 
we show the non-split models have some puzzles regarding the generation 
of charged matter fields. 


2.1. Puzzle on matter spectra near Dg points in Ig model 


In this section, we consider the six-dimensional F-theory on an elliptically 
fibered Calabi-Yau three-fold over a Hirzebruch surface F,, !-° in which the 
unbroken gauge symmetry is As or C3. These setups of F-theory compacti- 


are called codimension-one singularities®®!!. In general, there are intersections of the 
codimension-one loci in the base. These codimension-two loci, where codimension-one 
singularities are enhanced, correspond to intersections of 7-branes. So, The codimension- 
two loci are involved in matter generation and are called codimension-two singularities. 
This conifold singularity exists not only in the base but also in the elliptic fiber direction 


and appears naturally where matter generates without any special tuning of parameters. 
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fications are dual to the Es x Es heterotic theories on K3 with (12+n, 12—n) 
1-3. In the heterotic dual, these massless spectra can be calcu- 
lated using the index theorem (Table 1). They satisfy the anomaly-free 
condition for the Eg with 12 + n instantons ny — ny = 30n + 112°. 


instantons 


Table 1. The anomaly-free massless matter spectrum of heterotic theory on 
K3 in which the unbroken gauge symmetry is A5 or C3. Specifically, we con- 
sider the case where 12 +n instantons are distributed among (8 +n —r,4+r) 
(0<r<n+2) in (Ag, A1) or (G2, A1). 


Gauge group | Number of hypermultiplets: nz 
A5 320 + (n +2-— r)15 + (2n +16 +r)6 + (3n —r + 21)1 
C3 $14’ + (n+1-—r)l4+ (2n +164 3r)6 + (4n — 2r + 23)1 


In F-theory, an elliptic fibration can be described in Weierstrass form 
y’ =+ forty, (1) 
where x and y are coordinates describing an elliptic curve. And f and g 
are holomorphic functions on a Hirzebruch surface (a P fibration over P’) 
I . n 
fw) = £ 2" fa4n(4—i)(w), g(z,w) = 2 9124n(6-J) (w), 
i=0 j=0 
where w and z are the affine coordinates along the base and P' fiber, and I 
and J are the largest integers that satisfy I < 8, 8+n(4—I) > Oand J < 12, 
12 + n(6 — J) > 0, respectively. In this subsection, all subscripts denote 
the degree of the polynomial in w. We only consider i < 4, j < 6, and near 
z = 0 so that we focus on the Eg with 12 + n instantons in the heterotic 
dual! °. There are singular fibers over the zero loci of the discriminant 


A := 4f? + 279? =0. (2) 


These singularities can be classified by the vanishing orders z of f, g, and 
A on the codimension-one singularities (Kodaira’s classification”). 


We obtain a split Te equation (as Table 1) from the equation (1)°1?: 
fig(z,w) = —3téhá pao + 62th ppo Hn-r44 + 327t, Qteeahe r42 
tye n44) ol z” (tr fn—r+8 J 3Ur44Hn—r+4) “le tees 3) 
GI (z, w) = hd eee ai 620 Ne esta ea > OPE tee (gale r42 
=t, H? n44) T 2, (9ur+4hz pig Hn-r+4 = ey Rene eer 


2t, He r44) Ez (Bur, shine T tr fn-r+8Hn-r+4 
= fet hn-r+2 N 3trurpa H2 npa) 


+25 (fgtrHn—r44 =F Urtafn—r+8) + 9122°, (4) 
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Table 2. Kodaira’s classification of singularities of an elliptic surface”. 


ord(f) ord(g) ord(A) Fiber type Singularity type G 
>0 >0 0 smooth none 
0 0 n In An-1 
>i 1 2 II none 
1 >2 3 III A 
>2 2 4 IV A2 
2 >3 n+6 hs Dn+4 
>2 3 n+6 ry Dn+4 
>3 4 8 Iv* Es 
3 >5 9 III* Ex 
>4 5 10 II* Eg 
24 >6 > 12 non-minimal a 


Table 3. Massless matter content in the I model of F-theory, if the matter fields 
locally exist at all codimension-two singularities. In Is, each representation of [¢ 
simply decomposes into an irreducible representation of C3 at each zero locus. 


| Tg (As) Te (Ca) 
| Enhancement in I | Matter rep. | Multiplicity | Matter rep. Multiplicity 
IV* (Ee) 520 r Fs +6) r 
Iž (De) 15 n+2-r 14 2n+4-—2r 
I7 (A6) 6 2n+16+r 6 2n+16+r 
1 3n +21-r 1 An + 23 — 2r 


where tr, hn-r+2, Hn-r+4, Ur+a, fn—r+s, fg, and gig are the sections of 
appropriate line bundles over the base P’. And the discriminant (2) is 


Ar = 2°he-42Pantrsis +2 tha 42Qan420 + 2° Rans2a + O(2°), (5) 


where Pən+r+16; Q3n+20, and R4n+24 are some non-factorizable polynomi- 
als. At the zero loci of tr, hn-r+2, and Pən+r+16, the codimension-one 
singularity is enhanced to IV* (E6), I} (De), and Iy (Ag), respectively. 
These zero loci are the codimension-two singularities. In general, any two 
of tr, hn-r+2, and Pon+4,416 will not be zero at the common locus, so we will 
get the matter content (Table 3) from (3), (4), and (5). It can be verified 
from Tables 1 and 3 that if the matter locally exists at all codimension-two 
singularities, then the anomaly-free condition is satisfied. 

A non-split Ig equation, which reduces the gauge group from As to Cs, is 
obtained by the replacement of the section Re 42 > hon-2r+4 (not chang- 
ing tr and Pən+r+16) in the split Je equation since the two split exceptional 
curves can no longer be factorized®. Thus, we consider each matter repre- 
sentation of A; simply decomposes into an irreducible representation of C3, 
respectively (Table 3). In this case, the number of the zero loci where the 
codimension-one singularity is enhanced to Iž is doubled, but we cannot 
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assign it to two separate zero loci since 14 of C3 is not a pseudo-real but a 
real representation. Therefore, it is difficult to realize the local generation 
of matter at all codimension-two singularities in the non-split models. 


2.2. Intersections of exceptional curves in non-split Ie 


As we will see explicitly in the next section, we examine, in the non-split 
Ig model with a Dg codimension-two singularity, the intersection diagrams 
of exceptional curves near hon—-2r+4 = w = 0 (in Section 3, bzo = w = 0). 
At a point of w Æ 0, in the split Jẹ model, we have five exceptional 
pir Cpa» and Cy, (the top of Figure 1). In the non-split Js model, 
since the two split exceptional curves C% cannot be factorized in terms of 


Pl 


the polynomial ring of w and replace each other at w = 0, C7, are identified; 


the same applies to C,- This corresponds to the diagram automorphism of 
the Dynkin diagram: Cp, := į (C} +C;,) (i = 1,2) (the middle of Figure 1). 
At w = 0, no conifold singularity remains after the codimension-one 
%11 by lifting 
up the exceptional curves C’s and 0’s (as (x,y,z) > (#12, yiZ, 2) in Section 
3.2), we obtain Cp, > dp,, Cpo — Spy, Cpa + dp, (the bottom of Figure 1). 
Therefore, in the non-split models, even at the codimension-two singu- 


curves C 


singularity is blown up!?. Thus, as in the previous works 


larities, the intersection diagrams of exceptional curves remain the same. 
Thus there is no new two-cycle where an M2-brane can wrap around. 


eo—_e_e—e_® 


Ch, Cp, Cp; C, Cp, 
eee 
Cp, Cp, Cp; 
Cpi > Sp; 
o—_o——s 
Sp, Sp, Ops 


Fig. 1. Intersection diagrams of the exceptional curves (based on Ref. 12): (Top) w 4 0 
before the identification; (Middle) w Æ 0 after the identification; (Bottom) w = 0. 


3. Split/non-split transition as a conifold transition 


In this section, to show the split/non-split transition is a conifold transi- 
tion in the I2 model, we will perform a concrete blowing-up process of a 
codimension-one singularity, near a Dj; codimension-two singularity. After 
that, we will show the non-split models have nontrivial three-cycles; we will 
discuss the origin of non-local matter in the non-split models!?. Finally, 
we will discuss all other fiber types in which the non-split models exist 1°. 
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3.1. The local equation 


A local elliptic fibration with a codimension-one singularity along (x,y,z) = 


(0,0,0) can be described in “Degline form” 7°: 
b b b 
= 2) Big A a 6 
(x,y,z, w) = Y re 4 4” T 3” j 4 0, 
1 
bs := 7 (bobs — bi): (6) 


1 
A = = (abs — 9bzbabe + 8b + 2788) , 


where b; is a section. b; is expanded as b; = bjo + bj12 + + bijz (j = 
2,4,6), where b; x is a ((j — k)n + 2j)-th degree polynomial in w. 

There is a relationship between “Degline form” (6) and Weierstrass 
form (1) in the previous section as follows: f = -4 (b3 — 24b4), g = 
aa (b3 36b2b4 4 216b¢). f and g are sections of the same line bundle as 
b4 and bg, respectively. And then, in the dP» fibration they are expanded 
as f = fao + farz+---+ faaz*, 9 = 960+ 96,02 +--+ +96,62°, where fax, 
ge,k are written as fg4n(4—k), 9124n(6—k) in the previous section. 

The conditions for singularities of the split and non-split fiber types are 
summarized in Table 4'°. Thus, substituting bz := 4[(w—e€)(w+e)—z], ba := 
2z*, and be := 4z?* into equation (6), we obtain the local equation with an 
Izk codimension-one singularity along (x,y,z) = (0,0,0) for arbitrary w: 


(x,y,z, w) E +z? + [(w — e)(w + €) — z| £? + 2*o + z™ =0, 
bg = z7* [4(w — e)(w + €) — 1 — 4z], 
A = 2™ [(w — 6} (w + €)?(4w? — de? — 1) (7) 
2(w — e)(w + €)(6w? — 6e — 1)z 
+(12w? — 12 — 1)2? + O(2°)], 


l | 


where € € R is the parameter for a split/non-split transition since b2, can 
be factorized when e€ = 0 and not when € Æ 0. In this section, for simplicity, 
we consider k > 3, but we can make the same discussion for k = 2. 


3.2. Resolutions of Izk codimension-one singularities 


We consider the resolution of a codimension-one singularity? +! of the local 
equation (7) along (x,y,z) = (0,0,0) for arbitrary w. This equation (7) 
also has one D2, codimension-two singularity at w = +e and another one 
at w = —e in the non-split side, while it has one at w = 0 in the split side. 
So, we obtain the exceptional curves C’s at w Æ +e and ĝ’s at w = te. 
We replace (x,y,z) = (0,0,0) for some fixed w with a P?, by replacing 
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Table 4. Singularities of the split and non-split fiber types. I>), , denotes the “over-split” 
type which is explained in Ref. 19. 


erez ord(b2) ord(b4) ord(be) ord(bg) ord(A) ee ee 
To (k > 2) 0 k 2k 2k 2k b2,0 = Cio Tee 
b2,0 generic ge 
b2,0 = C10 
Ton4i(k > 1) 0 k 2k 2k+1 2k+1 ba, = C1,0C3,k Uspta 
b6,2k = ER 
b2,0 generic 
[ b4,k = b2,0C2,k Wea 
b6,2k = b2,002 k 
b2,0 = ĉio 
[ bak = Ch 0C2,k I3k+1 
been = ci oc k 
b21 = 4(p2,1 + q2,1 + 12,1) 
+ b4,2 = 2(p2,1q2,1 + 92,172,1 wg 
I 1 2 3 4 6 R Th 
be,3 = 4p2,192,172,1 
b2,1 = 4(p2,1 + 92,1) 
ba,2 = 2(p2,192,1 + 14,2) I5 
b6,3 = 4p2,174,2 
62,1, 64,2; b6,3 pens 
generic 0 
T3,-3(k>2) 1 k+1 2k 2k+1 2k+3 b6.2h = C3.b TE 
be,2% generic xR 3 
T3,_9(kK>2) 1 k+1 2k+1 2k+2 2k+4 bs,2h4+2 = CL pq TAPS 
bg,2n+2 generic en 
IV 1 2 2 3 4 b6,2 = C34 Iv* 
b6,2 generic Lys 
Iv* 2 3 4 6 8 bea = c33 Iv*s 
be,4 generic Iv*"* 


C? with Ĉ3 = {((#,y,2),(€:0: ©) € C3 x Pla: y : 2) = (Ein: O}. 
We will be blowing up the codimension-one singularity in inhomogeneous 
coordinates defined in the three different affine patches of P?, for example, 
(x: y:z)= Eine =r ror? 1) (1, z #0). Then, to replace C? 
with Ĉ3, we simply replace (x,y,z) with (x12, yz, z) in the equation (7). 
To not change the canonical class, the equation after the blow-up is defined 
as follows: 2~?®(x1,2, y12z, z,w) =: ®, (1,41, 2,w) = 0. The other patches 
are also similar. In this (7) case, the codimension-one singularity remains 
along (xj, yj,z) = (0,0,0) (j = 1,---,#—1). Thus, we perform similar 
processes k times for the resolution of the codimension-one singularity. 

Performing these operations, we obtain the j times (j = 1,--- , k) blown- 
up equation in patch 7,..., is 

Qe 


j 
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Delan us mw) = -47 +a] + [lw ew e) — ap pany + E 
a 
=0. (8) 
These equations have codimension-one singularities along (£j, yj,z) = 
(0,0,0) in j = 1,---,k— 1, and not in j = k. In the non-split side (e€ Æ 0), 
we also obtain the exceptional curves at w # +e and w = +e as follows: 


Cp, :2=0, y2? = (w— + 
gage, O 


bp, :Z=0, yy =0, w= e 
Co. 2=0, yg = (w — e)(w + e)z? + zk +1 
x 5 (Gj=k). (10) 
bp, 12 =90, Yk =k +1, w= te 


3.3. Local deformed conifolds in non-split models and 
split/non-split transition as a conifold transition 


To see conifold singularities, we consider the j times (j = 1,--- , k) blown- 
up equation in the patch of jz... zz 
SA 
j-1 
2 Mw? j—1 j 
C= Be -zalEj-1 Yj Zj w) = =y} — Žj£j-1 ee ey 
7-1 
k—j+1,(k-J) | ,2(k-J+1) ,.2(k—J) 
+z” any +z 7 T31 2-1) 
where (£o, Yo, zo) means (x,y,z). In the split side (e = 0), we can confirm 


that there are k conifold singularities at a D2, codimension-two singularity: 
Var : (x, Y1, 41, w) = (0, 0, 1, 0), 
Ug; * (£j—1, Yj, Zj, w) = (0,0, 0, 0) (j =2,---,k—1), 
Ury_1 | (Lk—2,Yk—1; Zk—1, W) = (0,0, 1,0). (12) 


Thus, the intersection diagrams of exceptional curves are obtained by blow- 
ing up all conifold singularities (small resolutions), which explain the en- 
hancements of the symmetries. The smooth split I2 model corresponds to 
the resolved side of a conifold transition. 

In the non-split side (e 4 0), by appropriate variable transformations of 
equations (11) and considering only their lowest order, we obtain 


yt aajitw-?=0 (f=1,---,k-1). (13) 


These are local deformed conifolds!®. Thus, it is natural that there is no 
conifold singularity in the non-split models; there are three-cycles instead 
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of two-cycles. The smooth non-split [2, model corresponds to the deformed 
side of a conifold transition. 

The split/non-split transition parameter € also plays the role of a coni- 
fold transition parameter, so the split /non-split transition is a conifold tran- 
sition at a Do, codimension-two singularity in the I2 model. This clarifies 
the deformation of conifold singularity is associated with a diagram auto- 
morphism of the expected simply-laced Dynkin diagram in the split side. 

Moreover, in the non-split models, it turns out that new S°’s, which are 
non-local in terms of the elliptic fibration, are generated. Thus, this can 
be regarded as the origin of non-local matter which cannot be generated 
from two-cycles. Since it is important to consider two Do; codimension-two 
singularities as a pair to find these three-cycles, the matter generated from 
three-cycles may satisfy the anomaly-free condition. However, this requires 
a new way to generate matter fields from three-cycles with finite size. We 
will leave this issue that should be clarified in the future. 


3.4. All other fiber types with non-split models 


Similar results are obtained in IV and I3,_; (k > 3) at a Dox point, and 
in IV* at an Ey point!®. In Inx_1 (k > 2), a split model has no Do, point 
in general. But we get a Do, point if we adjust a Do,_1 and an Aək—1 
to overlap at the same point. We call such a specially tuned fiber type an 
“over-split” 19. At this D2, point, we obtain a similar discussion through an 
“over-split” model. In I3,_, (k > 1), conifold singularities do not appear 
after the resolution of the codimension-one singularity. Therefore, this is 
the only case where similar arguments cannot be applied 1°. 


4. Conclusion 


In this paper, we have considered the non-split models have two puzzles 
on matter generation, using Ię specifically: first, if the matter fields locally 
exist and matter representations in the split models simply decompose into 
irreducible representations of the corresponding non-simply-laced gauge 
groups at all codimension-two singularities, the number of matter fields 
mismatches the anomaly-free condition®!* !®; second, there is no conifold 
singularity at all Do, or Ey points after all codimension-one singularities 
are blown up, so the intersection diagrams of exceptional curves remain 
the same!?. Next, using Ig, specifically, we have shown the split /non-split 
transition is a conifold transition at a D2, codimension-two singularity, and 
nontrivial three-cycles S? are generated in the non-split models. These have 
been clarified, in these cases, the deformations of conifold singularities are 
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associated with diagram automorphisms of the corresponding simply-laced 
Dynkin diagrams in the split sides. Moreover, we have discussed how these 
solve the above puzzles '®. However, in the 3, models, we cannot discuss 
the above. Also, we need a new understanding of matter generation due to 
three-cycles. These are issues that need to be clarified in the future. The 


analyses based on the deformations?!:?? could help in solving these issues. 
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Entanglement entropy in Schwarzschild spacetime 
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In this work, we discuss the entanglement entropy in the Schwarzschild space- 
time, and its relation to the vacuum state of matter fields. Recently, it was pro- 
posed that there is either large violation of the additivity conjecture in quantum 
information theory or disentangled states of black holes in the AdS/CFT cor- 
respondence. Here, we consider the additivity conjecture in the Schwarzschild 
spacetime. Usually, the entanglement entropy is calculated assuming that the 
vacuum state is in the Hartle-Hawking vacuum. We discuss the entanglement 
entropy in the Schwarzschild spacetime, and its relation to the vacuum state 
of matter fields. We show that the other vacua than the Hartle-Hawking vac- 
uum should be taken into consideration in order to consider the additivity 
conjecture. 


1. Introduction and summary 


Recently, it was proposed that there is either large violation of the additivity 
conjecture in quantum information theory or a set of disentangled states 
in the black hole spacetime [1]. A simplest statement of the additivity 
conjecture|[2—8] says that the minimum output entropy 


Smin N) = min S(N(p)) (1) 


of two quantum channels Mı and Mz which map states in Hilbert spaces A; 
or Az to those in another Hilbert spaces B, or Bə satisfies the additivity 
condition; 


Smin (Mi @ No) = Smin (M1) + Smin( N2) , (2) 


where S(p) is the con Neumann entropy of the density matrix p. 

Hayden and Penington studied the additivity conjecture in the 
AdS/CFT correspondence [1]. Two Hilbert spaces A; and A are identified 
to those of two conformal field theories which correspond to two boundaries 
of the AdS black hole spacetime. They considered quantum channels Mı 
and Mə which take the partial trace in each CFT. The output entropy 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
use, distribution and reproduction in any medium, even commercially, provided that the 
original work is properly cited. 
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of the channels are the entanglement entropy, which can be calculated by 
using the Ryu-Takayanagi formula [9, 10]. Then, they argued that the en- 
tanglement entropy does not satisfy the additivity condition (2). Hayden 
and Penington proposed that there is either large violation of the additiv- 
ity conjecture or a set of disentangled quantum states which correspond to 
disconnected geometries in the gravity side [1]. 

Here, we consider the additivity conjecture in the Schwarzschild space- 
time. We consider the quantum states of matters in two exteriors of the 
horizon in the Schwarzschild spacetime (See Fig. 1(left)). In a similar fash- 
ion to the case of the AdS/CFT correspondence, we take the partial trace 
of the states to obtain the entanglement entropy of the Hawking radiation. 
The entanglement entropy is usually calculated by assuming that the state 
is given by the Hartle-Hawking vacuum. The total entanglement entropy of 
the Hawking radiation is smaller than the sum of the entanglement entropy 
in each exterior, and hence, the additivity conjecture is not satisfied. 

This is because we considered only the Hartle-Hawking vacuum, and the 
other vacua should be taken into account. In the other static vacua than the 
Hartle-Hawking vacuum, the quantum energy-momentum tensor becomes 
very large near the horizon. By solving the semi-classical Einstein equation, 
we find that two exteriors of the horizon are disconnected (Fig. 1(right)), 
as is proposed by Hayden and Penington. Thus, it is important to consider 
general static vacua in order to reproduce the additivity conjecture in the 


Schwarzschild spacetime. 


Fig. 1. The Penrose diagram of the Schwarzschild spacetime in the Hartle-Hawking 
vacuum (left) and the semi-classical Schwarzschild spacetime in the other static vacua 
(right). 


1 
‘Island 7 
on 


aes 
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This paper is organized as follows. In Sec. 2, we briefly review the 
additivity conjecture in the AdS/CFT correspondence. In Sec. 3, we cal- 
culate the entanglement entropy in the Schwarzschild spacetime. In Sec. 4, 
we solve the semi-classical Einstein equation to obtain the semi-classical 
Schwarzschild solution. In Sec. 5, we consider the additivity conjecture in 
the Schwarzschild spacetime. This contribution is based on [11]. 
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2. Additivity conjecture in AdS/CFT correspondence 


The AdS/CFT correspondence relates the AdS spacetime with conformal 
field theory. In the cases of the eternal black holes in asymptotically AdS 
spacetime, the geometry have two boundaries, which correspond to two 
conformal field theories. The Einstein-Rosen bridge which connects two 
exteriors of the event horizon implies that states in two conformal field 
theories are entangled with each other. This entangled state is interpreted 
as the thermofield double state at the Hawking temperature, 


Ib) = $ PF In, )|n_) . (3) 


We consider the additivity conjecture for states in two CFTs A; and 
Az. Two quantum channels M; take partial trace in each CFT and maps 
states in total system of each CFT to those in subregions Bı and B2. The 
output entropy of the channels are nothing but the entanglement entropy 
of region B;, 


S(Ni) = S(Bi) , S(N, @ M2) = S(B, U B2) . (4) 


By using the Ryu-Takayanagi formula, the entanglement entropy of the 
region B in CFT is given by the area of the minimal area surface yg whose 
boundaries are anchored at the boundaries of B in the AdS boundary, 


_ Area(ys) 
zE (5) 


EE 6.0 
OC CE 


Fig. 2. The geometry of a time slice in the AdS black hole spacetime and the Ryu- 
Takayanagi surface for Bı or Bo (left), and that for Bı U Be (right). 


S(B) 


For the entanglement entropy of Bı, the minimal surface lies in one 
exterior of the event horizon (Fig. 2 (left)). For the entanglement entropy 
of Bı U B2, the minimal surface extends between two boundaries through 
the Einstein-Rosen bridge (Fig. 2 (right)) [12]. Thus, we have 


S(B, U B2) < S(B,) + S(B2) . (6) 
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Therefore, the additivity conjecture is violated by the typical states of black 
holes in the asymptotically AdS spacetime. 


3. Entanglement entropy in Hartle-Hawking vacuum 


Now, we consider the additivity conjecture for states of matters in the 
Schwarzschild spacetime. The Schwarzschild metric is given by 


ds? = — (1 = =) dt? + (1 Š -D dr? + r2dQ2 . (7) 


The quantum channels M; map states in total system into those of the 
Hawking radiation, which are defined as states in the region R, or equiv- 
alently in r > b. The entanglement entropy of the Hawking radiation is 
given by [13-17] 


Area non logal. 
S a 5 + XE Cees (8) 
aR or CN ap, or 
where 
(non-local) __ i—j41 © 
Š matter ra (= 1) j 5 2 logl( (U; — Uj \(Vi — Vj) t+ i2 3 D8 lguv (xi)| 
8A 
(9) 


and the summation is over the endpoints of region R and the island I 
[18-23] (if there is the island), which are labeled so that i — j = +1 for 
neighboring labels. For the Hartle-Hawking vacuum, the coordinates in (9) 
should be the Kruskal coordinates. 

We first consider the entanglement entropy of the Hawking radiation in 
two exteriors, Rı U R2. Before the Page time [24, 25], the configuration 
without the island dominates, and the entanglement entropy is calculated 
as [26] 


2mb? c 16r? (b — rp) to 
S= +=] h h? . 10 
Ge 6° | bD O a oy) 
After the Page time, the configuration with an island dominates. The 
position of the island is determined so that the entanglement entropy is 


extremized. The entanglement entropy is evaluated as [26] 
Qrr2 2mb? c 16r} (b — rn)? b—Tp, 
Sa + = |] h | 11 
GN $ GN 6 hoe ( b ) Th ( ) 


Next, we calculate the entanglement entropy of the Hawking radiation 
in one of two exteriors, Rı. There is no configuration with islands. The 
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entanglement entropy is calculated by introducing the IR cut-off A and is 
given by [27] 
mb? 


c 
= — + 4log A 
S Get ee : (12) 


which implies the entanglement entropy is infinitely large after taking A > 
oo. Therefore, the entanglement entropy of the Hawking radiation in the 
Hartle-Hawking vacuum does not satisfy the additivity condition, 


S(Rı U R2) & S(R1) + S(Ro) . (13) 


This apparent violation of the additivity conjecture is because we consider 
only the Hartle-Hawking vacuum. The Hartle-Hawking vacuum is given in 
the form of the thermofield double state in terms of the states in each of 
two exteriors. By using these states in each exterior, we can construct other 
static vacuum states. Since we need to consider the state with the minimum 
output entropy, the other vacua than the Hartle-Hawking vacuum must be 
taken into account. 


4. Semi-classical Schwarzschild spacetime 


In the other vacuum states than the Hartle-Hawking vacuum, the energy- 
momentum tensor on the classical Schwarzschild spacetime diverges at the 
event horizon. This implies that quantum effects becomes very important 
near the horizon, and we should solve the semi-classical Einstein equation, 


1 
zI R = 8TG y (Tuv) - (14) 


In the s-wave approximation, the expectation value of the quantum energy- 


Ruw — 


momentum tensor is completely fixed by the conservation law and Weyl 
anomaly. For the most general spherically symmetric spacetime, 


ds? = —C (u, v)du du + r° (u, v)dQ? , (15) 
the energy-momentum tensor is calculated as [28, 29] 


ck? 


_ 1/292 -1/2 |, 
(Tuu) 4872 aris Oye "192722 ’ mo) 
(Tos) = 01/2820 21/2: 4: cK (17) 
“= aa ie 1927? r2 ’ 
(Taw) = -ss (Cu C — O,C0,C), (18) 


sie aie 
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where « is the integration constant, which is related to the asymptotic 
condition. Then, it is straightforward to solve the semi-classical Einstein 
equation to obtain the expression near the Schwarzschild radius [30-33], 


C(r.) = e™/™ + O(a?) , (20) 
1 


rameter -a(g e 
Ar, 


where a = SGX and r, is the tortoise coordinate rą = }(v—u). The Hartle- 
1 


Hawking vacuum is given by k = x— where the radius approaches to the 
Schwarzschild radius rp in rą — —oo. In the other static vacua, K Æ 7o 
the radius diverges for k < = or goes to zero for K > a before C (rx) 
goes to zero. In either case, the geometry has no event horizon but naked 
singularity. Thus, two exteriors of the event horizon are disconnected in 
the other static vacua than the Hartle-Hawking vacuum. 

Away from the Schwarzschild radius, the metric is approximated by the 


classical Schwarzschild solution (7), where 


Th 


r(ra) 


rear = ra +n log (=) f (23) 
Th 


Clr.) =1- (22) 


5. Additivity conjecture in Schwarzschild spacetime 


Now, we calculate the entanglement entropy in static vacua with «x Æ ae 
As we have seen in the previous section, two exteriors are disconnected 
by taking quantum effects in the energy-momentum tensor into account. 
Then, the entanglement entropy of Ri U Rə is simply given by 


S(Rı U Rə) = S(R1) + S(R2) k (24) 


Thus, the additivity condition is satisfied. 

Next, we calculate the entanglement entropy to see that the entangle- 
ment entropy in a vacuum with Kk Æ Se gives the minimum output entropy, 
or equivalently, is smaller than that in the Hartle-Hawking vacuum. In these 
vacua, the coordinates in the formula (9) should be chosen as 


Uk te, Von te, (25) 


where (u,v) coordinates are related to t and r, as v =t+r, and u = 
t—r,. In the configurations without islands, the entanglement entropy is 
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calculated as [11] 


b? b— 
Sao log ( rt) ; (26) 
N 


This entanglement entropy becomes very large for very small «, and then, 
a configuration with an island gives the dominant saddle point. The entan- 
glement entropy for small « is approximately given by [11] 


Tb? c b—rph c 48rr? 
S= l l R), 27 
TE + to ( b ) + $ og ( Ear (27) 
The entanglement entropy (26) is smaller than that of the Hartle- 
Hawking vacuum for k > E and hence, gives the minimum output en- 
tropy. Therefore, the additivity conjecture is satisfied in the Schwarzschild 


spacetime if the other static vacua than the Hartle-Hawking vacuum are 
taken into consideration. 
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Holographic index calculation for Argyres-Douglas and 
Minahan-Nemeschansky theories 
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We calculate the superconformal indices of the Argyres-Douglas theories and 
Minahan-Nemeschansky theories realized on N coincident D3-branes in 7-brane 
backgrounds via the AdS/CFT correspondence. We include the finite N cor- 
rection as “giant gravitons”, which are D3-branes wrapping around 3-cycles. 
We take account of a single giant graviton for simplicity, and our method nicely 
reproduces known results and gives predictions for theories whose indices are 
unknown. 


Keywords: Argyres-Douglas theories; Superconformal index; AdS/CFT corre- 
spondence. 


1. Introduction 


In this proceeding, we calculate the superconformal index of the Argyres- 
Douglas theories and Minahan-Nemeschansky theories using a brane real- 
ization in the type IIB superstring theory. We consider a brane system with 
N D3-branes and a 7-brane. Let x°, x!, x”, and z? be the coordinates along 
D3-branes, and X, Y, and Z be the 3 complex coordinates of C? transverse 
to the D3-branes. The D3-branes are located at X = Y = Z = 0. We also 
introduce a 7-brane with the worldvolume Z = 0 (Table 1). 


Table 1. The brane setup. 

0 1 x Y Z 
7-brane V y vv 
D3-branes Vo Vv 


<<} 
< ssa 


There are 7 types of 7-branes with constant axiodilaton. Correspond- 
ingly, there are 7 types of 4-dimensional superconformal field theories, which 
we denote by G|N]. (G = Ho, Hi, H2, D4, Eg, E7, Es). This brane system 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
use, distribution and reproduction in any medium, even commercially, provided that the 
original work is properly cited. 
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is quater BPS and M = 2 supersymmetric theories are realized on the D3- 
brane worldvolume. H,,[N] (n = 0,1,2) are Argyres-Douglas theories)? 
and E,,[N] (n = 6,7,8) are Minahan-Nemeschansky theories?*. D4[N] is 
an SQCD with Sp(N) gauge group. H,,[N] and E,,[N] are strongly coupled, 
and it is worth analyzing these SCFTs via AdS/CFT correspondence. 
The 7-brane induces a deficit angle tag on Z-plane® ®. Table 2 shows 
ag and the gauge symmetry on 7-brane for each G. From the viewpoint 


Table 2. Theories corresponding to each dificit angle ag. 


G Ho Ay Hə D4 Es Ez Eg 
oe 13 1/2 2/3 I 4/3 3/2 5/3 
gauge sym. on 7-brane || None SU(2) SU(3) | SO(8) | Es E7 Esg 


of 4-dimensional theories, the gauge symmetry on 7-brane becomes fla- 
vor symmetry. The rotational symmetry along transverse directions is 
SO(4) x SO(2) pe SU(2)r x SU(2)r x U(1)r;, where SU(2)r x U(1)r, is 
R-symmetry, and SU (2) p is a flavor symmetry. Let Rx, Ry, and Rz be the 
generators rotating X, Y, and Z-plane, respectively. U(1)r, is generated 
by Rz, and the Cartan generators of SU(2)r and SU(2)pr are (Rx + Ry) 
and 3(R x — Ry), respectively. 

The near-horizon geometry of the brane system is AdS5 x S3 g, where 
See is a 5-sphere in XYZ space with deficit angle tag on Z-plane. The 
7-brane wraps around the singular locus Z = 0, and the gauge symmetry 
G lives on the locus. 


<> 4Y i 
X 
Fig. 1. 7-branes in S2 of G[N]. 


The superconformal index is defined by® 


rankG 
T=Ar epe JI d . (Uxruyuz = 1) (1) 
g= 
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where H is the Hamiltonian, J and J are the angular momenta, and p; 
are Cartan generators of G. It expresses the BPS spectrum concisely as a 
function of fugacities q, y, Uz, Uy, Uz and xi. The factor e2™/+/) is +1 for 
bosonic states and —1 for fermionic states. 

Table 3 shows the theories whose superconformal indices have already 
been calculated!®!4. The method explained below can be used to calculate 


Table 3. Theories whose superconformal in- 
dex have been already calculated. 


the leading finite N correction to all G|N]. 

In the following sections, we will briefly explain how to calculate the 
index on the AdS side and show a few results. For detailed explanations and 
more results, refer to the original paper !° in collaboration with Y. Imamura. 


2. Large N limit 


Let us first consider the large N limit, which has already been analyzed 
in16-17, There are two contributions to the index (Figure 2). 


(a) Kaluza-Klein modes of the gravity multiplet in the bulk 
Closed strings give the gravity multiplet in the bulk. They are ex- 
panded into spherical harmonics in 9° o: Due to the deficit angle, 
the boundary condition associated with the angular coordinate on 
Z-plane should be appropriately modified 1617. The superconfor- 
mal index can be obtained by summing up the contributions from 
all modes. 

(b) Kaluza-Klein modes of vector multiplets on the 7-brane 
Vector multiplets of gauge group G live on the 7-brane. We can 
expand them into spherical harmonics on the 7-brane worldvolume 
S? C $2. in a similar way to (a). The superconformal index can 
be obtained by summing up the contributions from all modes. In 
the D4, case, the 7-brane is a stack of an O7-plane and 4 D7-branes, 
and the vector multiplet comes from open strings on the D7-branes. 
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Fig. 2. The contributions in the large N limit. (a) The closed string in the bulk. (b) 
vector multiplets on the 7-brane. Open string description of (b) is justified only in the 
D4 case. 


3. Finite N corrections 


For finite N, we include extra contributions from “giant gravitons” 18, which 


are the D3-branes extended in $3 o- We follow the prescription proposed 
in!®. Namely, we consider D3-branes wrapped around three particular 3- 
cycles, X = 0, Y = 0, and Z = 0. It is expected the complete finite N 
index Zain] would be obtained if we sum up over all wrapping numbers 
(Nz, Ny, Nz) as 


(nanya) 
Ten] = TLE {oo} 5 Lown} , (2) 


Ng Ny Nz 


where T@o0] is the large N index, and a is the index from gi- 


ant gravitons with specific wrapping numbers. It is difficult to calculate 

Tenue) for general i b dh tak t of onl 
GIN] or general wrapping numbers, and here we take account of only 

leading corrections, ier and ie 


— (1,0,0) (0,1,0) 
Tan) = Zajo (l + Zain) + Zein) t (3) 


(Due to the deficit angle, the volumes of 3-cycles X = 0 and Y = 0 are 


smaller than that of Z = 0, and A becomes subleading.) T and 
(0,1,0) 
T, 


Gin] are related to each other by the Weyl reflection uz + uy. 
We can split T into 3 factors (Figure 3). 
(c) The classical factor from the maximal giant graviton 
The volume of the 3-cycle X = 0 is proportional to N, and a giant 
graviton wrapped around the cycle carries H = Rx = N. This 
corresponds to the factor qu. 
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(d) The vector multiplet on the giant graviton 
U(1) vector multiplet live on a giant graviton. Its fluctuations can 
be expanded into spherical harmonics. We sum up the contribu- 
tions from them. 

(e) The degrees of freedom along the intersection of the 


giant graviton and the 7-brane 
In the D; case, open strings between the giant graviton and the 


D7-brane give chiral fermions living along the intersection. Their 
contribution to the index is the character of the basic representation 
of the D4 current algebra. For G # D4, we cannot directly derive 
the contribution, and we assume that it is the character of the basic 
representation of the G current algebra. 


Fig. 3. Three ccontributins to the finite N correction. (c) The classical contribution. 
(d) The vector multiplet. (e) The intersection modes. Open string description of (e) is 
justified only in the D4 case. 


4. Results 


Now, we can calculate the index by combining the five contributions (a)-(e). 
We will first show the known result for Ho[1]1°. 

£ 6 Ł 17 -4 u 
Troji] =1 + Uz q5 m ue xłq™ + Uz °q 5 


12 ı2 & 2 T 
where xy/ and yf are the characters of the n-dimensional SU(2); and 
SU(2)r representations, respectively. 

Let us first compare this with the large N index without the giant 
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graviton contribution included: 


É o Ł 17 = =a 
Tmj] =1 + usg® — už xi qT +u xi + (uF xi + uz” )q? 
12 6 3 z 
+ Qu q” + (uf -uz xf xia” — Qukx{g +. (5) 
We find the discrepancy between (5) and (4) at the order q? as shown by 
the underline. 
Next, let us look at the result including the giant graviton contributions: 


É 6 + 17 =4 11 
Taa =1 + ub g’ — uk x{q™ +u: q? 
12 6 29 
+u? q? +ubxfq@+us q+. (6) 
The discrepancy becomes of order qe, which is higher than the previous re- 
sult. This implies the contributions from giant gravitons nicely reproduces 


the finite N correction. 
Finally, we show the result of E[1] 


Ads -1.E. E = -1.B 
Lap =1 + uz “xaasd + (—1-— X248)q? + (uz! + uz X248)X1 4 


9 
+ (uz + u2xe $o00)0' + (-2 — dis) X19? + Ol), (7) 


where 28 is the character of the n-dimensional irreducible Eg represen- 
tation. The expected error of our calculation is of order q°, and all terms 
shown in (7) are expected to be correct. This has not been obtained by 
other methods. 


E 
1 2 
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Topological defect junctions in 4-dimensional pure Zz lattice 
gauge theory 
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We explore topological defects and their properties in 4-dimensional pure Z2 
lattice gauge theory. This theory has the Kramers-Wannier-Wegner(KWW) 
duality. Duality defects associated with the KWW duality are constructed 
and shown to be non-invertible topological defects. In this paper, we explore 
the crossing relations including the duality defects. We construct 1-form Z2 
center symmetry defects and defect junctions. Crossing relations are derived 
from these defects and defect junctions. We also calculate some expectation 
values of topological defects by crossing relations. 


Keywords: Generalized global symmetries, Non-invertible symmetries 


1. Introduction 


The concept of symmetry in quantum field theories is very important. 
There are many applications to the non-perturbative analysis of quantum 
fields theories. In recent years, there has been progress in the generaliza- 
tion of the notion of symmetries. The important notion of generalizations 
is topological defects. One of the generalizations is so-called non-invertible 
symmetries. There are a lot of studies of non-invertible symmetries in 2- 


14 


dimensions! !*. Non-invertible symmetries in 2-dimensions are relatively 


understood than in higher dimensions. There are several studies of non- 
invertible symmetries in higher dimensions !> 2. 

In this paper, we study 4-dimensional Zə pure lattice gauge theory. 
There is a duality discovered by Wagner?’ in this model. This duality 
is similar to so-called Kramers-Wannier duality?4?°, therefore we call this 
duality “KWW duality”. Another important property is a 1-form Zo center 
symmetry 7°. The charged objects are Wilson loops. 

In this paper, we aim to investigate the crossing relations among topo- 
logical defects in 4-dimensional Zg lattice gauge theory based on this work?” 
which is in collaboration with Masataka Koide and Satoshi Yamaguchi. 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
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Crossing relations are local relations between topological defects. These re- 
lations determine the algebra of symmetry structure. The crossing relations 
are not closed within only duality defects. For the crossing relations, we 
need Z symmetry defects and junctions. Defect junctions occur when two 
types of defects meet. We determine the weight of junctions by junction 
commutation relations and find the crossing relations by these defects and 
defect junctions. Related examples are discovered by the works2*:?9. 

This paper is organized as follows. In Sec. 2, we explain our setup and 
KWW duality defects. In Sec. 3, we construct Zə symmetry defects and 
defect junctions. In Sec. 4, we explain crossing relations and calculate some 
expectation values of KWW duality defects. 


2. 4-dimensional Zə pure lattice gauge theory 


In this section, we explain our setup of 4-dimensional Zə lattice gauge the- 
ory. We prepare two kinds of 4-dimensional cubic lattices. One lattice is 
an active lattice A := {(£1, £2, £3, £4)|£1, £2, £3, £4 E 2Z}. Another lat- 
tice is an inactive lattice A := { (£1, £2, £3, £4)|£1, £2, £3, £4 € 2Z +1}. A 
schematic picture of these lattices is shown in Fig. 1 


Fig. 1. A schematic picture of our lattices. Even though this figure is drawn in 2- 
dimensional, the lattices are 4-dimensional. Black dots represent a lattice A. Blue dots 
represent a lattice A. These lattices are dual to each other. 


We assign the link variables Um = (—1)° = +1 to each link on A. 
We also call am(= 0,1) a link variable. We don’t assign the link variables 
to each link on A. There is a one-to-one correspondence between 16-cells 
and plaquettes on A. We regard the 16-cells as fundamental components of 
the total lattice. We assign the Boltzmann weight to the 16-cell with link 
variables a, = 0, 1(k = 1,2,3,4) and real parameter K as 


W (a1, a2, a3, a4) = exp(K (—1) (1 te2tastaa)), (1) 
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We define the partition function as 


Z=S\( [J sc UP oC I aC TL OT] Wenw, ene, ai, ai), 


{a} active active inactive inactive i€C 
site inks sites links 


(2) 
where C is the set of all 16-cells and j1(2), j2(i), js(2), j4(i) are links of a 
plaquette 7 in A. The partition function sums over all configurations of link 
variables. s and / are the weights of the site and the link on A, respectively. 
3 and / are the weights of the site and the link on A, respectively. 

We construct KWW duality defects on the codimension 1 surfaces. 
Weights of links, sites and KWW duality defects, and K are determined by 
commutation relations of KWW duality defects. Commutation relations 
are requirements that KWW duality defects can be smoothly deformed. 
Because we expect the KWW duality defects are topological defects, we 
require these commutation relations. The nontrivial solution is 


o 1 l= 1 
TENE YR 
The solution of KWW defects has the following properties. Wilson loops 
have nontrivial action of KWW duality defects. When the duality defect 
acts to a Wilson loop, a ’t Hooft loop appears where the Wilson loop 
was. The most important property is that KWW duality defects are non- 


. 1 
s=1, [=1, K 5 log(1 + V2). (3) 


invertible. 


3. Z2 1-form symmetry defects and defect junctions 


In this section, we construct Zo 1-form symmetry defects and defect junc- 
tions. There are 1-form Zo center symmetries”° in this model. We need 
to construct 1-form Zə center symmetry defects as codimension 2 objects. 
We consider that Zə symmetry defects are supported on triangles which 
are formed by a link in A and the midpoint of the adjusted link in A. We 
deformed these triangles to double the link in A and the midpoint of the 
adjust link in A on these triangles. On the other hand, the sites at both 
ends in A of these triangles are not doubled. These deformations are shown 
in Fig. 2. 

We assign weights to Z2(b, c) = of. = (1— ôb,c) to each building block of 
Zə symmetry defects. b,c = 0,1 are link variables in the building block. We 
also assign the weight z = V2 to each pair of each link on A. Zə symmetry 
defects assigned those weights are topological and invertible. Topological 
means that Zə symmetry defects can be smoothly deformed on octagons 
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|g 


Fig. 2. A schematic picture of the Zə symmetry defect. Blue dots and links represent 
sites and links in A. Black dots are the midpoints of links in A. We assign weights a 
triangular prism and a red line. We deform triangles to double the links and sites on 
triangles. 


surrounded by triangles. We can show these defects are invertible from the 
fact that the S? expectation value of a Zz symmetry defect is 1. When the 
Zə symmetry defects act on a Wilson loop, the sign flips. Therefore, we can 
consider the constructed Z2 symmetry defects are actually the symmetry 
defects associated with a 1-form Zə center symmetry. 

We consider configurations of defects where KWW duality defects and 
Zz symmetry defects meet. Junctions occur in such configurations. There 
are two types of junctions depending on whether a sharing link is in A or 
A. Fig. 3 and Fig. 4 show schematic pictures of two types of junctions. 
We denote each weight of junctions J(a) and J(b,c). a,b,c = 0,1 are link 
variables in these junctions. 


Fig. 3. A schematic picture of a junc- Fig. 4. A schematic picture of a junc- 
tion. We call its weight J(a). aisalink tion. We call its weight J(b,c). b,c are 
variable in this junction. link variables in this junction. 


These weights are determined by junction commutation relations. Junc- 
tion commutation relations are a requirement that Zə symmetry defects can 
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be smoothly deformed along with KWW duality defects as shown in Fig. 5. 
We expect there are topological junctions and for this reason, we require 
these junction commutation relations. The nontrivial solution is 


J(a) = (-1)", J(b,c) = oF... (4) 


Fig. 5. A schematic picture of the junction commutation relations. Green surfaces 
represent KWW duality defects. Red surfaces represent a Zə symmetry defect. This 
relation is a requirement that the junctions can be smoothly deformed. 


4. Crossing relations 


In this section, we explain crossing relations between KWW duality defects 
and Zə symmetry defects. Crossing relations can be explicitly calculated by 
weights of defects and junctions. We find three types of crossing relations. 

The first crossing relations are as shown in Fig. 6. Zə symmetry de- 
fects that have the end to KWW duality defects can be removed when the 
boundary of a Zə symmetry defect on KWW duality defects is homologi- 
cally trivial. 

The Second crossing relations are as shown in Fig. 7. We call these 
relations solid torus equations. We consider a decomposition of S? into two 
3-dimensional solid tori. We call these solid tori Vı and V3. We place a 
KWW duality defect on V; and this is the left-hand side of Fig. 7 . The 
first term of the right-hand side of Fig. 7 is the configuration of KWW 
duality defect on V2. The second term of the right-hand side of Fig. 7 is 
the configuration of KWW duality defect on V2 with a Z symmetry defect 
on Də whose boundary is a non-trivial cycle on V2. 


116 


A 


Fig. 6. A schematic picture of the crossing relations. Green surfaces represent KWW 
duality defects. A red surface represents a Z2 symmetry defect. A boundary of the Z2 
symmetry defect is on the KWW duality defect. When the boundary of Zz symmetry 
defects on the KWW duality defect is homologically trivial, the Z2 symmetry defect can 
be removed. 


Fig. 7. A schematic picture of solid torus equations. Green surfaces represent KWW 
duality defects. A red line represents a Z2 symmetry defect. The left-hand side is a 
configuration of a KWW duality defect on V;. The first tern of the left-hand side is a 
configuration of a KWW duality defect on V2 and the second term of it is a configuration 
of a KWW duality defect on V2 with a Z2 symmetry defect on D2 whose boundary is a 
non-trivial cycle on V2. 


The third crossing relations are as shown in Fig. 8. We consider the 
decomposition of $3 into S? x D? and D! x S?. We consider the configura- 
tion of KWW duality defects on S° x D? and D! x S?. Expectation values 
of these configurations are 1/\/2 times different. 

We can calculate some expectation values of a configuration of a KWW 
duality defect by crossing relations. 

One example is an expectation value of a KWW duality defect on S°. 
Let us consider two 9°? expectation values of KWW duality defects. Each 
S? has D3. Therefore we can use a crossing relation of Fig. 8. After using a 
crossing relation of Fig. 8, we can find the expectation value of S° is 1//2. 
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Fig. 8. A schematic picture of the crossing relation including S° x D3. 


This expectation value means KWW duality defects are non-invertible. 

Another example is an expectation value of a KWW duality defect on 
S? x $1. We can calculate this expectation value by using crossing relations 
of Fig. 6 and Fig. 7. A schematic picture of a calculation is as shown in 
Fig. 9. S? x St contains D? x St. This is a solid torus. Therefore, we 
can use a solid torus equation of Fig. 7 and we use this in the first equality 
of Fig. 9. A Zə symmetry defect on D? and its boundary is on a KWW 
duality defect on 9. We use crossing relations of Fig. 6 to this term. We 
can see that the S? x S1 expectation value is 1 by using the result of an S° 
expectation value of a KWW duality defect. 


S2 x St S8 


Fig. 9. A schematic picture of a calculation of an expectation value of a KWW duality 
defect on S? x S1. In the first equality, we use the solid torus equation of Fig. 7. In the 
second equality, we use a crossing relation of Fig. 6. In the third equality, we use the 
result that the S3 expectation value of a KWW duality defect is 1/\/2. 


5. Conclusions 


In this paper, we construct topological defects in 4-dimensional pure Zz lat- 
tice gauge theory. There are KWW duality defects and 1-form Z symmetry 
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defects. KWW duality defects are codimension1 non-invertible topological 
defects. We construct 1-form Z symmetry defects and defect junctions. 
Finally, we explain three types of crossing relations. These relations can be 
calculated by weights of junctions and defects. We calculate some expec- 
tation values of KWW duality defects by crossing relations. 
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We construct half-geodesic Witten diagrams with four external scalar fields 
in the Rindler-AdS black hole and show that their late-time behaviors agree 
with the Regge limit of conformal blocks. We also argue connection with 
pole-skipping phenomena by demonstrating that the near-horizon analysis in 
the Rindler-AdS black hole can determine the Regge behaviors of conformal 
blocks. This proceeding is based on a collaboration with Keun- Young Kim and 
Kyung-Sun Lee. 


Keywords: AdS/CFT correspondence; Black hole; Conformal block 


1. Introduction and summary 


The Regge limit of conformal four-point functions is a well-studied subject 
for consistency of quantum field theories. Out-of-time-order correlation 
functions (OTOCs) in conformal field theories (CFTs) are related to the 
Regge limit, and it has been proposed that exponential behaviors in four- 
point OTOCs are bounded by the consistency!. It has also been proposed 
that the exponential behaviors in theories with Einstein gravity duals can 
be calculated from pole-skipping points in retarded Green’s functions of the 
energy density. Here, the pole-skipping points are defined by intersec- 
tions between poles and zeros of momentum Green’s functions. 

In our paper*, from the viewpoint of holography, we study a relation 
between the Regge limit of conformal blocks and the pole-skipping points 
in CFTs on Rindler spacetime. First, we construct half-geodesic Witten 
diagrams in the two-sided Rindler-AdS black hole and demonstrate that 
their late-time behaviors agree with the Regge limit of conformal blocks. 
Second, we show that the near-horizon analysis in the Rindler-AdS black 
hole for the pole-skipping points can capture the Regge behaviors of confor- 
mal blocks. Our result is a generalization of the original proposal regarding 
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the pole-skipping points. 


2. Half-geodesic Witten diagrams and the Regge limit of 
conformal blocks 


We define scattering amplitudes of spin- exchange half-geodesic Witten 
diagrams in the Rindler-AdS black hole as 


WR e = | dX dX Goa (Y(A), W.L; Aw) Ga (Y(A), Wr; Aw) 
Y 


wer (rra e A.0) () 


x Goa (Y (0), Vr; Av) Goa (YQ’), Vr; Av) , 


where Gg is a scalar bulk-to-boundary propagator, and G, is a spin-f 
bulk-to-bulk propagator. Instead of the entire black hole spacetime, we 
integrate over two half-geodesics y8, and y between boundary points and 
centers of Penrose diagrams. After some approximations and calculations, 
we evaluate (1) in a late-time limit: 


2 
wE, = pauls) Can iy (2) 


Awi H(A — 1)” 


xelf-Itn—(A—-1)d oF (a —1, —-1A4+1- sie?) ; 


(2) 


which agrees with the Regge limit of conformal blocks®®. Here, tg and 
d are differences of time and space between two operators in four-point 
OTOCs. 


3. Near-horizon analysis in the Rindler-AdS black hole 


Let us consider an ansatz of symmetric traceless spin-f fields 


hy..op (v, T, x) = eu eG _ 1RD os (x) ’ (3) 
j=0 


where v is a coordinate in the incoming Eddington-Finkelstein coordinates 
of the Rindler-AdS black hole. By substituting this expansion into the 
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classical equations of free spin-£ fields, we obtain 
0=(V,V" — A(A — d) + Dhe... (v, r, x) 
=e” Dua — (l+ iw)(d — 1) — A(A — d)] AY, (x) 
-27 fiw + (€ — 1)] AG, x) — e fiw + (€ — 1] AM, (x) +e. 
(4) 


The near-horizon analysis is a holographic method to search pole-skipping 
points by imposing that coefficients in the classical equations are zero’. 
From the near-horizon analysis of (4), we obtain conditions for the pole- 


skipping points: 
-iw = 2-1, [Oe — (A -1)(A-d+1)]h®, =0. (5) 


These conditions are consistent with the Regge behaviors of conformal 
blocks (2) with respect to tr and d6, and therefore the pole-skipping 
points in CFTs on Rindler spacetime are related to the Regge limit of con- 
formal blocks. 
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Entanglement entropy (EE) in field theory is a measure for quantum entangle- 
ment between spatially separated regions. While there are a lot of studies on 
EE in CFTs and free theories, EE in general interacting field theories requires 
further investigation. It is of very interest in order to relate the effect of the 
entanglement with low-energy physics. In this talk, we introduce our study on 
EE in interacting field theories with a subregion of a half-space. There, spe- 
cific contributions to EE can be expressed in terms of renormalized correlation 
functions of operators. The contributions are expected to be dominant when 
we discuss low-energy effective theories. 


Keywords: Entanglement entropy; Field theory; Correlation functions 


1. Introduction 


In the last few decades, there has been much investigation on the entan- 
glement entropy (EE) , which is one of standard measures for the bipartite 
quantum entanglement between subsystems. In particular in field theory, 
while two causally separated regions are irrelevant at the classical level, the 
quantum entanglement between them can bring some nontrivial influence 
on observations made in the subsystem. EE can capture such a nontrivial 
quantum effect, and has been studied in many contexts: quantum phase 
transition, information paradox of the blackhole, holography, and so on. 
So far, EE has been widely discussed in CFTs 1? or perturbations from 
them*®, and in free theories’ 13. On the other hand, EE in interacting 
field theories which are far from CFTs is rather less understood. There are 
conceptional and technical difficulties. EE measures the entanglement with 
a physical cutoff scale, and we are often interested in not the entire EE 
but its “universal” part. In an interacting field theory, we are to deal with 
divergent radiative corrections by renormalization, and the physical quan- 
tities are described without explicit appearance of the cutoff scale. Here, 
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a natural question is how the universal part of EE is related to the renor- 
malized physical quantities. It is of course crucial to directly investigate 
the physical implication of EE to realistic observations. Despite its impor- 
tance, not so many studies have been made on it because it is very difficult 
due to the lack of any useful symmetries or simplifications. On the other 
hand, it was reported in Ref. 14 that a part of EE in a scalar field theory 
is expressed with the renormalized mass parameter in a specific setup. 

In this talk, we present our series of works on EE in interacting field 
theories!’ +’, where we apply a technique called the orbifold technique 18419. 
It was originally proposed to evaluate EE in free field theory in the case 
of A being the flat half-space, and we apply it to interacting field theories. 
There, by investigating the general structure of Feynman bubble diagrams, 
we check that the area-law of EE holds at all order. Then we extract an 
essential part of EE which comes from two-point correlation functions of 
various operators. The two-point functions naturally contains the radiative 
corrections and is described by renormalized parameters. As a result, we 
associate a dominant part of EE with renormalized correlation functions 
of fundamental and composite operators. We show our analysis with an 
example of an interacting massive scalar field theory, while the basic idea 
is applicable to more general theories, as long as we consider the half-space 
subregion. 

In Section 2, we explain our application of the orbifold technique to 
interacting field theories, and derive the area-law of EE to all order. We 
describe in Section 3 our main result, where EE is associated with the 
renormalized correlation functions. Section 4 is devoted to the summary 
and outlook. 


2. Orbifold technique in interacting field theories 


In field theory, a state is defined on a Cauchy surface, namely a spacelike 
hypersurface. Consider separating the surface into two subregions: A and 
A. Those who can make an observation only in A find him/herself in the 
state pa = Trg ptot, where prot is the density matrix for the total system. 
The corresponding EE is defined as the von Neumann entropy of pa: 


Sa = —Tra (pa log pa). (1) 


A standard way to calculate it is the replica trick’, where we consider an 
n-fold covering. It consists of n replicas of the Euclidean system sewed up 
together at the subregion A. EE is then obtained by calculating the free 
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Fig. 1. Our setup to investigate EE. The subregion A is a half-space, whose boundary 
is expressed with a red line. The coordinates perpendicular and parallel to the boundary 
are x and zj, respectively. 


energy of the theory on the n-fold, taking the derivative of it with respect 
to n, and taking the limit n — 1. 

In particular, when we focus on the case where A is a half-space in d- 
dimensional flat space (see Fig. 1), we can employ a variation of the replica 
trick called the orbifold technique. In this technique, we take an analytical 
continuation on n to 1/M with a integer M. Then, EE is computed as 

ð (MFO) p 

~ OM M-1 l ) 

where F™) is the free energy of the theory on an orbifold R? /Zm. In the 

following, we represent the two-dimensional coordinates perpendicular to 

the boundary as Œ = (x1, %2), and the others parallel to the boundary as 
TIJ ` 

The advantage of introducing the orbifold is that we can formulate the 
theory on it almost in the same manner as on the ordinary flat space. 
The only difference is that we introduce a projection operator, that is a 
symmetric summation of a rotation operator: 


M-1 
5 1 . Ge — sin 22 zı 
P = — me H Ei z M Mi s 3 
mî © sin +4 cos 4 £2 (3) 
As is mentioned above, we discuss the example of a massive scalar field 
theory. The propagator on the orbifold is represented as: 


© 


, Ma. 

Go" (2,9) = = Yo Golz, yim), (4) 
m=0 

- dip ei((G"P)-e—p-y)+ipy (x —y)) 

Go(z, y; m) -|i d+ ptm (5) 
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It is easy to see that Go with nonzero m does not conserve the perpendicular 
components of the momentum. Rather, it carries a twisted momentum: 
(p,p) —> (g'™p,p). On the other hand, interaction vertices consist of 
coupling constants and the (d+ 1)-dimensional delta functions representing 
the momentum conservation. Note that the factor of the delta function for 
p’s is invariant under a simultaneous rotation by ĝ: 


Ô’ (pı + p2 +--+) = 6 (9(pi + p2 +--+). (6) 


By using this invariance and Eqs. (4) and (5), we can prove that the general 
formula for a L-loop bubble diagram B, which contributes to the free energy, 
takes the following form: 


V M-1 

d—1 = 

a M 2 aie mz), (7) 
mi; mL= 


L 


L 
B(m,-+- ,mz) = [Tepe pnts? (So = ve) x (8) 
l=1 l=1 
Here, V4—ı denotes the area of the boundary 0A, Va—ı = 57~1(0), whose 
divergence is not important in the following discussion. B is a contribution 
from a fixed configuration of twists (m,---,mz). I(p) is an ordinary inte- 
grand for the diagram. It is essential that the configuration is characterized 
by the twists of the loop momenta. 
Among the configurations of twists, the trivial one (mı = 0,- , mg; = 
0) yields the delta function with vanishing argument. Therefore, the con- 
tribution to the free energy is proportional to (Vy_1 x 8? (0)) = Va41, the 
volume of the bulk. However, it does not contribute to EE, because its 
M-dependence comes only from the prefactor in Eq. (7), and it is canceled 
in Eq. (2). On the other hand, all the other configurations have nontrivial 
arguments in their delta functions, and they are proportional to the area 
of the boundary. By combining this fact with the dimensional analysis, 
we conclude that the leading contribution from B to EE is proportional to 
Va-1/ e171, Since this is a statement for the general bubble diagram, we 
have seen that the area-law for EE holds at all order. 


3. Two specific contributions to the entanglement entropy 


In principle, we can obtain the free energy on the orbifold by computing 
B’s and summing them over the twists configurations. This is, however, a 
technically difficult task even for diagrams with a few loops. Instead, we can 
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Fig. 2. Examples of the configurations of twists we call the propagator contribution. 
The blue cycles and letters denote the twists on the corresponding loop momenta. The 
twists in the upper diagrams are realized by the twists on a propagator as shown in the 
lower diagrams. 


extract some specific contributions of physical importance. In the following, 
we discuss them and later their dominance in the whole contribution. 
While the twists are acting on the loop momenta in each B, in some 
configurations we can attribute them to one twist acting on a momentum 
on a single propagator, as is shown in Fig. 2. Let us call such class of con- 
tributions the propagator contribution. The point is that we can sum them 
up in a systematic way, as is explained below, and can straightforwardly 
compute their contributions to EE. In those diagrams, most of the parts 
are identical to those in the flat space; only one propagator gets replaced 
with Eq. (5). As a naive attempt, let us sum them up by three steps: (i) 
considering the bubble diagrams in the flat space, (ii) taking the derivative 
of them with respect to the propagator, and (iii) reconnecting them with 
Eq. (5). Note that in Step (ii), we obtain the interaction correction part 
in the exact two-point function G. Then, their contributions to FM) and 
EE would be summed up to take the following form: 
M-1 
FO) ~= | attnatty (G5(G—G)GS") Y Golm), 0) 


m=0 


S 1/e d?@-1 
Sprop,int ~ naf ort [UGo + ([Go)? + + (BGo)° qp | i (10) 


In Eq. (9), we have omitted the arguments (x, y) in the both factors. In Eq. 
(10), we have represented the contribution in the momentum representation 
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and expand it with respect to the 1PI part ©. The arguments of all the 
terms in the square bracket are (k = 0, ky). 

In fact, this is not correct and there is a subtlety in the attribution of 
the configurations to a twist on a single propagator. In some diagrams, the 
twists can be attributed more then one propagator. The method of taking 
the variation, as we have considered above, counts each of such attributions 
as distinct contributions. However, since what we have to sum up actually 
is the twist configuration, it leads to over-counting. Thus, Eq. (10) requires 
to be modified. Fortunately, we can show that the over-counting factors for 
diagrams are exactly classified according to the 1PI expansion. It means 
that we obtain the correct propagator contribution by dividing each term 
in Eq. (10) by the corresponding over-counting factor: 


Vii f° dtk 


1 1 
Sprop,int = ETa [co + zC)’ + 3(2Go)” +e 


(27)! 
ee 1/e diik 
= | ayer all — 5G). (11) 


Furthermore, we can deal with the free part on the same foot. It is 
because its contribution to the free energy, (1/2) log Go 1 is expanded by 
Schwinger parametrization to the series of ring-shaped diagrams, each of 
which consists of the propagators simply connected at their endpoints. As 
a result, their contribution to EE is identical to the one derived in the 
literatures. By combining it with Eq. (11), the propagator contribution is 
summed up to take the form of 


Danm i atk 
ep r)a! 


Note that the radiative corrections entering G have been calculated in the 
diagram in the flat space, and the final expression does not depend on the 
fact that we have introduced the orbifold in the middle of the analysis. 
After all, a potential danger of the additional counterterm on the singu- 
larity of the orbifold decouples, and we can consider the counterterms and 
renormalization just as in the flat space. Therefore, Eq. (12) is indeed the 
contribution written with the renormalized two-point function. 

We also extract contributions to EE other than the propagator con- 
tribution. Consider a class of twist configurations as is shown in Fig. 3. 
In each of them, the twists are attributed to a twist of a loop momen- 
tum flowing at some channel in a single vertex. We call them the vertex 


log G(k = 0, ky). (12) 


contribution. It can be explicitly derived almost in the same way as the 


131 


= 3) (3 


sÍ 


Fig. 3. Examples of the configurations of twists we call the vertex contribution. The 
twists in the upper diagrams are realized by the twists on a channel of a vertex as shown 
in the lower diagrams. The dot lines denotes the delta function to decompose the vertex. 


propagator contribution. 

We take the variation of bubble diagrams with respect to the vertex, 
namely the bare correlation function. It gives the exact correlation func- 
tion. Here is a difference; we organize the endpoints into two, and reconnect 
them with the twisted version of the delta function. In general, there are 
choices of how to organize the endpoints. It corresponds to the channel 
we consider, at which the momentum gets twisted. Accordingly, the orga- 
nization turns the correlation function into the two-point function of the 
composite operators. The connection by the twisted delta function comes 
from the fact that the momentum twist act on the channel of the vertex, 
which is expressed by decomposing the vertex with the delta function. In 
this case, we have the subtlety and over-counting, and have the correspond- 
ing 1PI-like expansion, and obtain the correct contribution by dividing the 
terms by the over-counting factor. 

This computation results in the formula for the vertex contribution, 
which is written in terms of the renormalized two-point correlation function 
of the composite operators. For instance, in ¢4-theory, we have :¢7: as the 
only composite operator, and the vertex contribution to EE is given by 


Va m d’-tky 


Set =f y 


log Gg2y2(k = 0, kj), (13) 


where G'g242 is the two-point correlation function for :¢?:. In general in- 
teracting theories, we have several vertices and each of them has multiple 
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channel to be considered. These give rise to many composite operators, 
which are mixed with one another. Even in such cases, we can go through 
the same discussion and get a similar expression of the vertex contribution. 
The operator mixing is taken into account by representing the correlation 
functions in the form of a matrix, whose elements are correlation functions 
of operators. As for operators with nonzero spin, we have to take trace 
over the space of spins. It leads to the introduction of some coefficient ma- 
trix. Together with the propagator contribution, which is nothing but the 
contribution from the fundamental operators, the contribution we discuss 
is written as 


Vai [5 dki, n A 
Sprop+vert = naf o" log G, (14) 
Gab = Go,o,(k = 0, ky). (15) 


Here, Go, o, denotes the renormalized correlation function of operators Oa 
and O+», that is normalized to be dimensionless by the cutoff e. ¢ is the 
numerical coefficient matrix, which is determined by traces over the spaces 
of spins. In the case where the operators are all scalar, ĉis a unit matrix. 

The result shows that the contributions of EE we have discussed is 
understood in terms of the renormalized correlation functions of various 
operators. This is an all-order analysis, and implies that a part of divergence 
in naive calculation of EE comes from the ordinary radiative corrections, 
and in a sense irrelevant once we study EE in terms of renormalized physical 
quantities. On the other hand, the explicit appearance of the cutoff scale 
in Eq. (15) simply suggests that EE should be measured with some scale to 
be finite. Note that even in the free field theory, we have to introduce some 
cutoff scale to define EE. Therefore, it is rather a matter of the definition 
of EE. 

There are many contributions to EE other than those we have discussed. 
They are expected, however, to be sub-leading, since all of the correspond- 
ing diagrams are of relatively higher loops. Their contributions to the free 
energy should be less dominant. In particular, when we consider a theory 
described by an Wilsonian effective action, where the quantum fluctuation 
has already been integrated out, such contributions are absent in the first 
place.* It is thus certain that Eq. (15) captures the dominant of EE. 


*In that case, we have infinitely many operators from vertices in the effective actions. 
However, since the correlation between higher dimensional operators tends to have a 
shorter correlation length, we can approximate EE by a finite number of the low- 
dimensional operaotors. 
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4. Summary and outlook 


In this talk, we have discussed EE in interacting field theories with the 
subregion of a half-space. We have reported in this case that a part of 
EE can be expressed in terms of renormalized two-point functions with the 
fundamental and composite operators. The composite operators arises in 
the decomposition of interaction vertices with all possible channels. This 
contribution to EE is certain to be dominant, and if we treat the Wilsonian 
effective action, it captures the whole EE. We expect that this result should 
be a significant step to associate EE to observations in a realistic physics. 

The most crucial question is whether and how we can generalize the 
result to the cases of the general subregions. Of course, our analysis with 
the orbifold technique highly depends on the fact that the subregion is the 
flat half-space. On the other hand, the final result of Eq. (15) it self does 
not. We further expect that, the qualitative implication that the dominant 
part of EE is described with the two-point functions of operators, should 
hold in more general cases. The reason is that once we express Eq. (15) 
in the position space representation, we do not see the restriction of the 
form of the subregion. A possible way to investigate it is to reformulate 
things directly in the position space representation. It is a challenging but 
interesting future work. As another direction, It is important to study 
the structure of the renormalization in detail. While our analysis holds to 
all-order, it is within a perturbative picture. It is worthwhile to pursue a 
nonperturbative proof for our statement. 
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In this paper, we explain how to derive the AdS5 xS supercoset sigma model 
in the Green-Schwarz formalism. The derivation is a generalization of the 
procedure for the PCM case developed by Delduc et al [arXiv:1909.13824]. 
This paper is based on the original paper [arXiv:2005.04950]. 
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1. Introduction 


Integrable systems give a fascinating playground for investigating nonlinear 
dynamical systems. Determining whether a given theory is classical inte- 
grable requires to construct a Lax pair which ensures the existence of an 
infinite many conserved charges, but finding a Lax pair is a difficult task 
due to the lack of a guiding principle. Recently, a nice way to handle this 
issue has been proposed by Costello and Yamazaki! which is based on a 4D 
Chern-Simons (CS) theory?:? with a meromorphic 1-form w. According to 
their proposal, by specifying the 1-form w and the boundary conditions of 
the gauge field of the 4D CS theory, one can systematically construct a 2D 
classical integrable field theory. This Costello- Yamazaki proposal has been 
explored in subsequent works, and it has been shown that 4D CS theory can 
describe various 2D integrable field theories and integrable deformations in- 
cluding the Yang-Baxter (YB) deformation* ® and the \-deformation!?™. 
For a recent review on this subject, see?. 


In this paper, we will explain how the AdS; xS? supercoset sigma model 
can be derived from a 4D CS theory by following the original paper!® 
(see also!*). As mentioned in the above, it is necessary to determine a 
memorphic 1-form w. Remarkably, this 1-form w is identified with a twist 
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published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
use, distribution and reproduction in any medium, even commercially, provided that the 
original work is properly cited. 
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function characterizing the Poisson structure of the integrable system by 
Vicedo!®. The Poisson structure of the AdS5 xS supercoset sigma model 
has been studied in!®!”, and we can read off the twist function from their 
results. Then by considering a 4D CS theory with the twist function and 
specifying the boundary conditions we give the AdS5x S supercoset sigma 
model in the Green-Schwarz formalism. The derivation we use here is a 
generalization of the procedure developed by Delduc, Laxcroix, Magro and 
Vicedo!®. 


This paper is organized as follows. Section 2 explains how to derive 2D 
NLSMs from 4D CS theory. In section 3, we derive the AdS; x S? supercoset 
sigma model by generalizing the preceding result on the PCM. Section 4 is 
devoted to conclusion and discussion. 


2. 2D integrable field theory from 4D CS theory 


In this section, we will give a brief review about a derivation of 2D integrable 
field theories from a 4D CS theory by following 1*8. 


2.1. 4D CS theory 


We start with introducing the action of 4D CS theory. Let G? be a com- 
plexified semisimple Lie group, and g? be the Lie algebra for G? with a 
non-degenerate symmetric bilinear form (-,-) : g7 x g7 + C. Then the 
action of the 4D CS theory is given by 

S[A] = -> wACS(A), (1) 

4T JMxcPi 

where a 2D Minkowski space M has a metric q; = diag (—1, 1) with coor- 
dinates (t,o) and a global holomorphic coordinate of C C CP! = CU {00} 
is denoted by z. The gauge field A on M x CP! is taking a value of g©, 
and C'S(A) is the CS 3-form 


CS(A) = (a dA + ZA ^ A) (2) 


Here, w is a meromorphic 1-form on CP! 

w = p(2)dz, (3) 
where (z) is a meromorphic function on CP!. It is noted that the mero- 
morphic function (z) can be identified with the twist function character- 


ising the Poisson bracket associated with the 2D integrabe field theory 15. 
For later use, we denote the sets of zeros and poles by 3 and p, respectively. 
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The fact that w is proportional to dz indicates that the action (1) is 
invariant under the transformation 


A= A+xdz, (4) 


where y is any g©-valued function defined on M x CP!. Hence, by using 
the gauge symmetry (4), the gauge field A can be expanded as 


A=Aidot+A_do + Az dZ. (5) 


where we introduced the light-cone coordinates defined as o® = 4 (rt +0). 
In the following discussion, we will ignore the z-component of A. 


Equations of motion 
Next, let us give equations of motion for the action (1). 


By taking a variation of the action (1) with respect to the gauge field, 
we obtain the bulk equations of motion 


Fy-=0, wF;=0, F(A)=dA+AQAA, (6) 


and the boundary equation of motion 
dw \ (A, 6A) =0. (7) 


Here, the factor w in (6) is kept because there is a possibility that 0: A4 are 
distributions on CP! supported by 3. Note that the boundary equation of 
motion (7) has the support only on M x p C M x CPt. This follows from 
the fact that since dw can be expanded as 


dw = Ozp(z) dz A dz, (8) 


only the poles of p can contribute to the distributions according to the 
relations 


t-o=-eS(4), rep. (9) 


za 


2.2. Lax form 


Classical integrablity of a given 2D classical field theory requires the ex- 
istence of a Lax pair £ on the 2D spacetime which is valued on g© , and 
satisfies the following properties: a) it is on-shell flat, b) it meromorphically 
depends on an auxiliary Riemann surface. In this subsection, we will ex- 
plain how the Lax pair underlying a 2D integrable field theory is introduced 
in the context of the 4D CS theory. 
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By following t18, the Lax form £ can be introduced by performing the 
formal gauge transformation 


A= dg +L, (10) 
where a smooth function g : M x CP! > G®, and we take a gauge choice 


In order to see whether the g©-valued 1-form £ introduced in (10) satisfies 
properties a) and b), we rewrite the bulk equations of motion (6) in terms 
of £ as 

O4,£_ —0_£L44+[£L4,£_] =0, (12) 


The first constraint (12) is regarded as the on-shell flatness condition for 
the Lax form on M. The second constraint (13) indicates that £ is a 
meromorphic 1-form with poles at the zeros of w. In this way, the 1-form 
L satisfies the properties a) and b), and can be regarded as the Lax form 
on M by identifying z on CP! with a spectral parameter. 


2.3. From 4D to 2D via the archipelago conditions 


Here, let us give a general formula for the 2D action of the underlying 
integrable field theory by performing a dimensional reduction of the 4D 
action (1) along CP?. 


For this purpose, by substituting (10) to (1), we rewrite the 4D action 
(1) in terms of £ as 
a 


S{A] = wh Igetl= zl, a NMI AGE). (14) 
xCPl 


4T J MxcP1 
In order to reduce the 4D action (14) to a 2D one on M, we need to impose 
the archipelago conditions 18 
defined as follows: There exist open disks Vz, U, for each x € p such that 
x € Vz C Uz and 


on ĝ. The archipelago conditions for ĝ are 


i) U, NU, =o if x fy for all x,y E€ p, 
ii) ĝ = 1 outside M x UzepUz, 
iii) ĝļ|mxu, depends only on o* and the radial coordinate |s| where €, is 
the local holomorphic coordinate defined as €, = z — 2, 
iv) ĝļlmxv, depends only on o*, that is, ĝs = ĝmxv, = ĝ|mMx {z} - 
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Here, we assume that g satisfies the archipelago condition (see, for exam- 
ple,'® for the existence of such g). Then, the 4D action can be reduced to 
the 2D action by performing an integral over CP!1® as follows: 


S LA = D (resz (9 L), gz dga) 
rep 


+5 5 (resz w) f Iwz [gx] : (15) 


aes Mx[0, Ra] 
The action (15) is invariant under the transformation 

gz œ Izh, L= hLh+h tdh, (16) 
where h : M —> G®. The symmetry can be regarded as the residual gauge 
symmetry of the Lax form. 


As discussed in!%, in order to determine a 2D integrable field theory, 
it is necessary to specify not only the boundary conditions, but also the 
behavior of ĝ. This point will be explained in the next section. 


3. The GS action of the AdS5 x SŽ supercoset sigma model 
from the 4D CS theory 


In this section, we will give a derivation of the Green-Schwarz (GS) ac- 
tion of the AdS; x SË supercoset sigma model from the 4D CS theory by 
following 13. 


3.1. The AdSs x S° supercoset sigma model 
Here, let us review a supercoset construction of the AdS; x S? superstring. 

The classical action of the AdSs x S° superstring in the GS formalism 
has been constructed based on the following supercoset !9 

PSU (2, 2|4) 
SO(1,4) x SO(5) ` 

As is well known, the super Lie algebra su(2,2|4) can be decomposed to 
vector subspaces with respect to the Z4-grading structure: 


g=9 7 GgY og Gg,  g® = so(1,4) x s0(5), (18) 


(17) 


where g® @ g®) and g™ @ g® are the bosonic and fermionic parts of 
su(2, 2/4) , respectively, and the commutation relations of g‘”) satisfy 


(gh, g] E g™ (m+n=k mod 4). (19) 
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Note that the Z4-graded property ensures the classical integrablity of the 
AdS; x S? superstring”?. 


The GS action (32) of the AdS; x S° supercoset sigma model is given 
by 
Stal =f Strides) dot dor, (20) 
M 


where j = g~ ‘dg € su(2, 2|4) is the left-invariant current for g € SU (2, 2|4). 
Here, d+ are the linear combinations of the projection operators Pix) : g > 
g™ (k = 0,1, 2,3) like 


d+ = Pa) + 2P) F Po) - (21) 


The associated Lax pair is?° 


+ (a GO 4271 EE ag”) do~ . (22) 
We can see that the on-shell flatness condition for (22) is equivalent to the 
equations of motion of the action (20). 


The Poisson structure of the AdS5 x S? superstring has been considered 
int®17, By computing the Poisson brackets of the Lax pair (22), we can 
obtain the twist function of the AdS; x S° supercoset sigma model* 


23 
val?) = aye (23) 


The twist function (23) has the following poles and zeros: 


p={+1,-1,+i,-2}, 3 = {0, co}, (24) 


where the poles are double poles and the zeros are triple zeros, respectively. 


3.2. A derivation from the 4D CS theory 


Now, let us reproduce the GS action (20) of the AdS; x SË supercoset sigma 
model from the 4D CS theory (14) with the meromorphic 1-form 


W = Ystr(z) dz. (25) 


* ystr(z) is slightly different from ¢string(z) in (2.10) ofë. These are related via ystr(z) = 
+ dstring (2) . 
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In this case, the gauge field A in (14) takes a value in g = su(2,2|4), and 
the bracket (-,-) in the 4D action (1) is replaced by the supertrace Str. 


As mentioned in the previous section, we need to choose a solution to 
the boundary equation of motion to specify a 2D integrable model. The 
associated boundary equations of motion for the meromorphic 1-form (25) 
are 


e ( (Ai, O¢, Ai), (Aj, 8,43)» =0, DEP, (26) 
where the double bracket is defined as 


((2,y),(0",¥)))p = (aes ws) Ste(a - 2’) + (resp Epos) (Stra y") + Stefa" + y)) 
a 2 (Str(x-y’) + Str(a’-y)) . (27) 


A possible solution to the boundary equationa of motion (26) is given by 
Aļz=p = 0 (p Ep). (28) 


This boundary condition obviously solves the boundary equation of motion 
(26) and leads to the GS action of the AdS; x S? supercoset sigma model 
as we will see in the following discussion. 


Next, let us construct the associated Lax pair by solving the constraint 
(13) with the boundary condition (28). For our purpose, we consider the 
following ansatz for the Lax pair as 


L= T yi! + yP! +z vi +27 vi") dat 


+ (a2 vE + tv 4 yM 4 eV) don, a 


where VI” (n = —1,0,1), V” : M — su(2,2|4) are undetermined func- 
tions. It is easy to see that the ansatz (29) satisfies the bulk equation 
of motion (6) with the twist function (23). By solving the relation (10) 
between the Lax pair and the gauge field at each pole under boundary 
conditions (28), vE are determined as follows 


y a Jit a a E E + Ja, pital — Jit — J24 + 93,4 — Ja, 
E 4 il ae 4 , 
vil — JLz — t j2, — J3,+ + ija, yea Jz t ij2,+ — J3,+— ija, 
to = 4 , £ _ 4 
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where we denoted the left-invariant curret g~'dg at each pole of the twist 
function (23) by 

ji =O “dglers jo = G9 1dG\ =i, (31) 
J3 ae | | ee JA = dolz 


Note that the above ansatz (29) is not the only possible solution that sat- 
isfies the boundary conditions (28). In fact, we can also consider Lax pairs 
for the pure spinor formt. 


Next, we substitute the Lax pair (29) into the master formula (15), and 
then obtain the 2D action 


Slax] = =/ Str] 5 (o.,+ = (1+ t)5o(2),4 + 5o(3),4 — (1 - ‘)Je(4),+)Jo(),- 
16 Jm = 


—(-deq),- + (1 = 4)Jjo(2),— — Jo(3),- + (1 + Dieta- jeans dot Ado , 
(32) 
where o € S* is a cyclic permutation of the set {1, 2,3,4}. The action (32) 
is clearly invariant under the cyclic permutations of jẹ. 


As discussed int’, it is necessary to impose a relation between jp (k = 
1,...,4), so that the 2D action (32) can reproduce the GS action of the 
AdS; x S? supercoset sigma model. This can be achieved by requiring the 
following relation, which respects the cyclic symmetry of the 2D action (32): 


jr = 70) (k=1,...,4), (33) 


where j € ṣu(2,2|4) is the left-invariant current for g € SU (2,24), and 
the map fs : su(2,2|4)© — su(2,2|4)€ is an automorphism of su(2,2|4) 
satisfying the Z4-grading property f+ = Id. For our purpose, we will take 
the Z4-grading automorphism f, such that each subspace g) (k = 0,1, 2,3) 
is the eigenspace of f, satisfying 


fog) = ig”. (34) 


The explicit expression of fs can be written down after taking a supermatrix 
realization of su(2,2|4) (For the details, see?!). 


By imposing the constraints (33) satisfying (34), we can easily show that 
the 2D action (32) and the Lax pair (29) reduce to the ones (20), (22) of 
the AdS; x SŽ supercoset sigma model in the GS formalism, respectively. In 
this way, the 4D CS theory with the moromorphic 1-form (25) can describe 
the AdS; x S? supercoset sigma model. 
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4. Conclusion and Discussion 


In this paper, we have reviewed how a 2D classical integrable field theory 
can be systematically derived from a 4D CS theory by following!®. Then, 
by generalizing the previous procedure on the PCM case, we have repro- 
duced the AdS;xS° supercoset sigma model from a 4D CS theory with a 
meromorphic one-form (25). In addition, we should note a recent work ‘4 
which has succeeded in introducing the world-sheet metric into an underly- 
ing 2D integrable field theory by extending a 4D CS theory. Thus, we can 
treat AdS;xS° superstring theory beyond the AdS; xS° supercoset sigma 
model in the context of the 4D CS theory. 


It has been investigated how 4D CS theory describe various 2D classical 
integrable theories. As a next step, it is natural to examine ways to quantize 
2D classical integrable field theories and their quantum integrability in this 
framework. The problem of quantization has been extensively investigated 
for so-called ultralocal models, such as the Faddeev-Reshetikhin model?? 
and massless Thirring model. Note that these 2D theories can be derived 
from 4D CS theories with a meromorphic 1-form w without zeroes +2425, It 
is known that ultralocal integrable field theories can be described by taking 
the appropriate continuous limit of integrable lattice models???3. Since 
4D CS theory can also describe integrable lattice models in terms of the 
expectation value of the Wilson loops*, it would be an interesting problem 
to understand this picture in the framework of 4D CS theory. 


When w has zeros, e.g. PCM or coset sigma model, the quantization 
of the corresponding classical integrable theory becomes more subtle (For 
this issue, see™!5 in the context of the 4D CS theory). Quantizing non- 
ultralocal classical integrable theories in a first-principles way is still an 
open problem and an important task. 
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Quantum entanglement is closely related to the structure of spacetime in quan- 
tum gravity. For quantum field theories or statistical models, we usually con- 
sider base space entanglement. However, target space instead of base space 
sometimes directly connects to our spacetime. In these cases, it is natural to 
consider a concept of target space entanglement. To define the target space en- 
tanglement, we consider a generalized definition of entanglement entropy based 
on an algebraic approach. This approach is reviewed and is applied to the first 
quantized particles, in particular, fermions. This article is based on the paper 
JHEP 08 (2021) 0461. 


Keywords: Entanglement; Matrix models; Target space; Mutual information. 


1. Introduction 


It is widely believed that quantum entanglement is closely related to the 
structure of spacetime in quantum gravity. In the AdS/CFT correspon- 
dence, the Ryu-Takayanagi formula? states that entanglement about the 
base space in holographic CFTs is connected to the area of minimal surface 
in the bulk. As in this example, we often consider the base space entangle- 
ment in quantum field theories or statistical models. However, target space 
instead of base space sometimes directly connects to our spacetime, for ex- 
ample, perturbative string theories or matrix models. Thus, it is natural to 
investigate a notion of target space entanglement®°. See also recent Refs. 
1,6-9.* 

In Ref. 3, the target space entanglement is defined using an algebraic 
approach. We will review this approach in Sec. 2, and apply it to quantum 
mechanics of fermions in Sec. 3, Sec. 4 and Sec. 5. 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
use, distribution and reproduction in any medium, even commercially, provided that the 
original work is properly cited. 

aA concept of entanglement in string theories (matrix models) is investigated in 10 and 
revisited in 11. 
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2. Definition of entanglement entropy based on subalgebras 
of operators 


Let us recall the conventional definition of entanglement entropy (EE). 
Suppose that a total density matrix p is given for a Hilbert space H = 
Hp ®Hg. The EE for subsystem Hg is defined as the von Neumann 
entropy of the reduced density matrix pg = trg p as Sp = — tre pp log pB- 
This definition relies on the tensor product structure of the Hilbert space, 
H =Hp®Hg. However, total Hilbert spaces sometimes do not have such 
simple tensor-factorized forms. For example, the Hilbert space of a first- 
quantized (non-relativistic) particle in a space R@ is schematically given by 
a “direct sum” as H = span{|x)|xz € R¢}. Thus, even if we divide the space 
R into two subregions R? = BU B, it is difficult to take the “partial trace” 
on B. 

The algebraic approach enables us to define EE without relying on the 
tensor product structure (see, e.g., the references 112-14), The algebraic 
definition is based on the subalgebra of operators (observables). If a total 
density matrix p is given, and we have a restricted set of operators (sub- 
algebra A), an entropy S.4(p) associated with the subalgebra A is defined. 
This concept is natural, if we recall the meaning of entropy in information 
theory. The entropy is a measure of uncertainty about the whole informa- 
tion when we can only know partial information. If an observer can use 
only a subset of operators A, the whole information is not obtained. En- 
tropy S.4(p) quantifies the amount of uncertainty (or unknownness). In this 
sense, the usual EE, Sg = —trg pg log pg, for H = Hp ® Hg represents 
uncertainty for an observer who can probe only subsystem Hpg. That is, it 
is the entropy for the subalgebra L(Hg) ® 1y,- The choice of subalgebra 
A is arbitrary, and we do not need the tensor product structure. 

For general subalgebra A, the entropy S4(p) is computed as follows. 
First, the ‘reduced density matrix’ p4 is uniquely determined from p and 
A as an operator in A satisfying the following equation: 


tr(p40) = tr(pap) “OEA. (1) 


For example, if the total Hilbert space has a tensor product form as H = 
Hp ® Hg, and we take the subalgebra A as A = L(HB) 8 lup, then pa 
is given by pp ®1y,/dimHg. The point is that the definition Eq. (1) is 


bHere, L(V) denotes a set of linear operators on linear space V, and 1y does the identity 
operator on V. 
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applicable even when the Hilbert space does not have the tensor product 
structure. 

Furthermore, for a given subalgebra, we can decompose the Hilbert 
space into blocks of tensor products where the subalgebra acts nontrivially 
only on each tensor component as follows: 


H= Hs, 8Hp, st. A=PL(Hz,)® la. (2) 
k k 


This decomposition is uniquely fixed by the subalgebra A. We represents 
the projection onto each block by IIx. We define the density matrix pz on 
the projected space IIH as 


1 
prk := —Ilkpllk, (3) 
Pk 


where pk is a normalization factor defined as p, := tr(I,pIl,) and is a 
probability of being in the sector I],H for the given p. Since the projected 
space IIH has a simple tensor-factorized form as IIH = Hg, ® Hp, in 
the decomposition (2), we can consider the reduced density matrix of pp on 
Hp, as 


PB, ‘= tT B, Pk- (4) 
Then, the ‘reduced density matrix’ p4 satisfying Eq. (1) is given by 
15 
pa = D Pr PB (5) 
B i dim(Hp, ) ) 
We define the reduced density matrix pg on space Hg = @, HB, as 


pp := Q proB,- (6) 
k 


EE S4(p) is defined as the von Neumann entropy 
S.a(p) = -trs pplogpp = — X pr log pr +X. peS(pp,); (7) 


k k 
where S(pB,) := — tre, PB, l0g pB,- The first term in the r.h.s. of Eq. (7) 
is called the classical part, 
Sei(p, A) : ae log Pk, (8) 


and is the Shannon entropy of the probability distribution {px}. On the 
other hand, the second term in the r.h.s. of Eq. (7) is called the quantum 
part S,(p,A). The expression in Eq. (7) is similar to the symmetry resolved 


entanglement entropy 15-16, 
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2.1. Example: Entanglement in a single qubit 


As a concrete example of EE in the algebraic approach, we consider a single 
qubit. The Hilbert space is two-dimensional space, H = span{|0),|1)}. We 
usually consider entanglement between two qubits. The algebraic approach 
enables us to consider “EE” even for a single qubit. 

The full set of operators £(H) is L(H) = span{I, oz, oy,02}.° If we take 
the subalgebra as this full algebra, the decomposition Eq. (2) is trivial as 


H=HEC (9) 


with A = L(H) @1. In this case, pg in Eq. (6) is just the original p. Thus, 
the EE associated with the full algebra £(H) is just the von Neumann 
entropy of p, 


Secu (p) = — tr plog p. (10) 


In particular, if state p is pure, the entropy vanishes as Scu) (p) = 0. It 
means that the pure state is not ambiguous and is completely determined 
by quantum tomography if we can use any operators. 

Situation changes when we can use only a subset of operators. Let us 
suppose that we can probe only z-direction. This corresponds to taking 
subalgebra A = span{l,a,}. The decomposition (2) for this choice of the 
subalgebra is H = span{|0)} ® span{|1)} where A = span{1,0;} can be 


represented as A = span { iG a ® span { (l i \ The projection I; 
are II, = |k)(k| (k = 0,1). We then have py = (k|p|k) and pg, = Ik. 


Thus, the EE associated with the subalgebra A is 


Sa(p) = —po log po — pi log pr, (11) 


where the quantum part S,(p,A) always vanishes, and the entropy is just 
the classical Shannon entropy of the probability distribution that the qubit 
is measured in 0 or 1 for the given state p. Even pure states in general have 
non-vanishing entropy (except for the case where states are eigenstates 
of oz). The non-vanishing entropy reflects the fact that pure states are 
ambiguous for restricted observers who can probe only z-direction. In fact, 


Ie pata. 2 f : 0 
the observers cannot distinguish pure states with mixed states p = & ) : 
Pı 


1 0 


“We take a basis such that oz = e = 


) with oz|k) = (—1)*|k) (k = 0,1). 
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3. Entanglement of fermions with a fixed number 


We now consider target space entanglement of first-quantized N fermions 
by the algebraic approach. The Hilbert space of the single particle is rep- 
resented by H). It is given by H™ = span{|x)|a € M} where M is the 
target space of particles. The Hilbert space H? of N fermions is given by 
the N-th exterior power of H® as 


HM = A Ho, (12) 


which is spanned as H? = span{|z1) A- A |aw)|ai € M(i = 1,..., N)}. 

We take a subregion B in the target space M, and consider the EE of this 
subregion. Since the Hilbert space HA? does not have a tensor-factorized 
structure with respect to the target space coordinates, we adopt the alge- 
braic approach instead of the conventional definition. The subalgebra we 
take is the set of operators acting non-trivially only on particles in subregion 
B, which is represented by A(B). For example, when N = 1, the subalge- 
bra A(B) is given by A(B) = span {|y)(y'||y,y’ € B} ® span { fz dz|z) (z|} 
where B is the complement region of B. For general N, we can decompose 
HW) into a direct sum of the following subsectors as 


N 
HM = AHP. (13) 
k=0 


The subsector HW) consists of states where k particles in B and N — k 
ones in B as 


ny = span {|21) A---A|rn) |a1,.-.,0% € B, thy1,-.-,0n E€ B}. (14) 


To represent the subalgebra A(B), we introduce the following abbreviated 
notation: 


{}n) = lar) Ass Alan) E€ N HO. (15) 
The subalgebra A(B) is then given by 


N 
A(B) = QB Ar, (16) 
k=0 


where A; is a subalgebra on HW? and takes the form 


Reis { [dere dena Abe) (La Bel A (Chad) 


with Y1,---3 Yk: Yh Yk © B. (17) 
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Since the subalgebra A(B) is specified, we can compute the entropy 
S'4(B) associated with this subalgebra in the manner described in the pre- 
vious section.? We call this entropy the target space entanglement entropy 
Sp because the subalgebra A(B) is characterized by the subregion B in the 
target space of particles. In the second quantized picture, we can define 
the conventional entanglement entropy for subregion B. The target space 
EE Sz agrees with this base space EE 134, 


4. Fermions in the Slater determinant states 
To be more specific, we focus on pure states whose N-body wave functions 
are given by the Slater determinants as 


nee E Sa detil) ganai 18) 


where x;(x) are one-body wave functions normalized as 


| aexi@xs(@) = by. (19) 
M 


The target space EE for subregion B can be evaluated as Eq. (7) by 
computing p and S(pz,) for the pure states w. After some computations 
(see Ref. 1 for details), we can find that the entropy Sp follows the simple 
formula: 


Sp() = —tr[X log X + (1n — X)log(Lw — X)], (20) 


where X is a N x N matrix given by 


Xy= | dex. (21) 


We call X overlap matrix. It is easy to show that the eigenvalues A; of the 
overlap matrix are in the range 0 < àA; < 1. 

From the formula (20), we can find that the entropy has the upper 
bound® as 


Sp(w) < N log2. (22) 


The maximum entropy N log 2 is proportional to the number of particles 
N, and thus follows an extensive property like thermal entropy. However, 


4In this case, the projection H, in Eq. (3) is the projection to the subsector He in 
Eq. (13). 
eWe can also confirm that the classical part Se; is bounded as Sa (p; A) < O(log N).+ 
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this upper bound is too generic. We expect that EE for ground states is 
not extensive but sub-extensive in local models. In fact, we will see in the 
next section that the entropy of a ground state of N free fermions behaves 
as S ~ O(log N) in the large N limit. 


5. Entanglement for free fermions in a circle 


We now apply the formula (20) to N free fermions in a circle with length 
L, i.e., the target space M is a circle. The Hamiltonian is given by H = 
pee 2 and we consider its ground state. The one-body eigenfunctions 
are given by y,(v) = eF ”2 where n are integers. Supposing that the 
total number of particles N is odd (N = 2K +1), the N-body wave function 


for the ground state is given by the Slater determinant as 


1 
Plz tw) = = YS sgn(o)xX-K« (tea) Xk (tom) (23) 
VN! oESN 
Thus, the target space entanglement for a subregion B can be obtained by 
the formula (20) with the N x N overlap matrix 


Rie Ve de x2, (0)Xn (2), (24) 
K. 


where n,n’ runs in —K,..., 


5.1. Single interval 


In this subsection, we consider the case where the subregion B is a single 
interval J; in the circle. We parameterize the length of the interval as rL 
(0<r<1), i.e., r is the ratio of the interval to the circle. 

In the large N limit, the asymptotic behavior of the entropy can be 
obtained as 


1 
Si, ~ 3 log[2.N sin(ar)| + Yı (25) 
with 
TS if dw —" — log (> iw) 0.495018. (26) 
-œ cosh*(7w) T (4 — iw) 


We show the plot of the entropy with the large N result (25) in Fig. 1. 
It shows that the entropy is sub-extensive (not proportional to N). Fur- 
thermore, the large N behavior Eq. (25) agrees with the EE for the single 
interval in c = 1 CFTs on the circle!” if we regard N as a (dimensionless) 
cutoff. 
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2.55 


eee 
20 40 60 80 100 


Fig. 1. EE for the half region. The red dots are the EE S for N = 1,3,--- ,101. The 
blue curve represents the large N result (25) with r = 1/2. 


5.2. Entanglement entropy and mutual information for two 
intervals 


In this subsection, we consider two disjoint intervals J; and Jo in the circle. 


Suppose that the coordinates of the circle is x moving in —£ <a< Z 


with the periodic condition x ~ x + L. We take the two intervals as I 
(FL, -7 L) and b = (SFL, $F L). 
The EE for the subregion Jı U Ig can be analytically computed in the 


large N limit! as 


Stub ~ | 2log[2N sin(nr)] + log sin[r(d + r)] sin[r(d — r)] 


3 sin? (7d) aes 
(27) 
We can also evaluated the target space mutual information; 
(h; Ig) := S(hh) + SU2) — S(h U I). (28) 
The large N behavior is 
. 2 
Thi ta) ~ f log sin[r(d tar —r)] 2) 


The mutual information is finite even in the large N limit. In addition, 
Eq. (29) agrees with the result in a c = 1 CFT (free compact boson at the 
self-dual radius18), although the reason is not understood well. 

The plot of the target space mutual information (28) is Fig. 2. 


fSee Ref. 1. 
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0.004 + N = 101, r= 0.01 


0.002 + 


= 


Fig. 2. Mutual information for two intervals. We take N = 101 and set the parameter r 
as r = 0.01 (length of each interval is rL). The red dots represent the mutual information 
for some values of d. The blue curve represents the large N result (29). 


6. Brief conclusion 


The algebraic approach is a powerful method of characterizing entangle- 
ment. This approach might be useful beyond the target space entanglement. 
A similar idea to define entropy based on observables is also investigated 
as the observational entropy (see, e.g., Ref. 19). 

We have used the algebraic approach to define the target space entan- 
glement of particles. In particular, we consider non-interacting fermions, 
which can be regarded as the singlet sectors of one-matrix models. It is 
more interesting to consider entanglement in multi-matrix models, and its 
relation to holography. 
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Jacobi Forms on Even Positive Definite Unimodular Lattices 
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We study Jacobi forms associated to even positive definite unimodular lattices, 
in particular Eg lattice and the Leech lattice. We require the Jacobi forms to 
be invariant under the orthogonal group of the lattice, in particular Eg Weyl 
group and Conway group Cog. These objects are of interest in both number 
theory and string theory, for example elliptic genera of E-strings, Monster 
CFT and so on. We establish the explicit structure theorem and determine the 
generators for the Hg case with index t < 13 and the Leech case with index 
t <3. This is a short summary of two joint papers’? with Haowu Wang. 


Keywords: Jacobi forms; Even positive definite unimodular lattice; Eg lattice; 
Leech lattice; Conway group. 


1. Introduction 


Jacobi forms are some fundamental objects in number theory and also ba- 
sic tools in string theory whenever a torus T? and a global symmetry is 
involved. Originally, Jacobi forms are introduced by Eichler and Zagier in 
their monograph. These forms are holomorphic functions in two variables 
(7,2) € Hx C which are modular in 7 with respective to SL2(Z) and quasi- 
periodic in z. Later, Gritsenko*+ defined Jacobi forms of lattice index by 
replacing z with many variables associated with an integral positive definite 
lattice. The Jacobi form creates an elegant bridge between different types 
of modular forms. For example, Jacobi forms can be identified as vector 
valued modular forms through the theta decomposition. Jacobi forms also 
have many applications in mathematical physics, such as the elliptic gen- 
era of K3 surface, elliptic genera of 6d (1,0) SCFTs and the topological 
string partition functions on various Calabi-Yau threefolds. It is a natu- 
ral question to determine the structure of the space of Jacobi forms. This 
question was solved by Wirthmiiller® for Jacobi forms associated with root 
systems not of Eg type, where there is polynomial ring structure. For the 
Eg root system, the structure becomes much more complicated due to the 
even unimodular property. In particular, the space is not a polynomial 


© 2022 World Scientific Publishing Co. Pte. Ltd. This is an Open Access book chapter 
published by World Scientific Publishing Company. It is distributed under the terms of 
the Creative Commons Attribution 4.0 International (CC BY 4.0) License which permits 
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ring. This was studied before inf, but the structure is still not explicit. 
Besides, little is known about spaces of Jacobi forms associated with other 
irreducible lattices of large rank. We aim to tackle these problems, in par- 
ticular the Eg lattice and the Leech lattice. 


2. Jacobi Forms of Lattice Index 


We review the basic notions of Jacobi forms of lattice index. Let L denote 
an even positive definite unimodular lattice equipped with bilinear form 
(—,—). Here the unimodular means the lattice has determinant 1 and the 
even means all vectors have even norms. It is well-known such L can only 
exist in dimension 8k,k € Z. 


Definition 2.1. Let k € Z be an integer and t € N be a non-negative 
integer. If a holomorphic function y : H x (L ®C) — C satisfies the 
conditions 


(i) Quasi-periodicity: 
o(7,3 +27 +y) = exp (—tni(x,x)7 — 2tmi(x,3)) p(7,3), x,y EL, 
(ii) Modularity: for a,b,c,d € Z and ad — bc = 1, 


p (= a ) = (er + d)Fexp (mi2) p(T, 3), 


ctr+d’cr+d cr +d 


and the Fourier expansion of y takes the form 


Ley sae, eae" Cae ey, 
n=0 LEL 
then it is called a weak Jacobi form of weight k and index t. If p further 
satisfies that f(n,/) = 0 whenever 2nt — (¢,2) < 0, then it is called a 
holomorphic Jacobi form. If p is invariant under the orthogonal group of 
lattice O(L), then it is called a O(L) invariant Jacobi form. Here O(L) 
contains all automorphism of L that keeps the bilinear form. 


As L is unimodular, the theta decomposition (see? Corollary 2.6)) yields 
that every O(L) invariant weak Jacobi form of weight k and index 1 is a 
holomorphic Jacobi form and can be expressed as g(T)Oz(7,3), where g(T) 
is a modular form of weight k — $rk(L) on SL2(Z), and ©z(T,3) is the 
Jacobi theta function of L defined by 


Ox (7,3) = ‘> eT i(lse)r+2mi(e,3). 
LEL 
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For index t > 1, the structure of O(L) invariant Jacobi forms becomes 
nontrivial. The cases of interest are those with large O(L), i.e. L is highly 
symmetrical, then the space of O(L) invariant Jacobi forms will be relatively 
simple. 

We focus on two most interesting cases. The first is the Eg root lattice, 
where we require the Jacobi forms to be Eg Weyl invariant: 


p(t, 0(3)) = (7,3), o © W(Es). 
Es Weyl invariant Jacobi forms are frequently used in string theory, for 
example in E-string theory, which is the simplest six dimensional (1,0) 
superconformal field theory with manifest Eg symmetry®!9:!!. Due to this 
importance, it is highly interesting to determine the precise structure of 
space of Eg Weyl invariant Jacobi forms to high index. 

The second case is the Leech lattice, which is the unique even positive 
definite unimodular lattice in dimension 24 that has no roots. It was discov- 
ered by Leech!”, and its uniqueness was proved by Conway 1. This lattice 
has many remarkable properties. For example, it plays a role in construct- 
ing the fake monster Lie algebra’ and proving the monstrous moonshine 
conjecture !°, and it achieves the densest sphere packing in 24 dimension 1°. 
The group Cog is the automorphism group of the Leech lattice, whose struc- 
ture was first described by Conway !”. It is natural the require the Jacobi 
forms on the Leech lattice to be Conway invariant: 


Y(T, o(3)) = y(T,3), a € Coo. 

The quotient of Coo by its center gives a sporadic simple group of order 
4,157, 776, 806, 543,360,000. Therefore, the Leech lattice is highly sym- 
metrical, and we expect that the space of Conway invariant Jacobi forms 
will not be too large. Due to the importance of the Leech lattice and the 
Conway group, we also expect that Conway invariant Jacobi forms will have 
some applications in mathematics and physics. These motivate us to study 
such Jacobi forms. 


3. Eg Weyl Invariant Jacobi Forms 


Denote the vector spaces of W(£g)-invariant weak and holomorphic Jacobi 
forms of weight k and index t respectively by Jg i aS Bie len Let 

M,(SLə(Z)) = C[E4, Es] be the ring of modular forms on SL2(Z). For 
fixed index t, these are free M,.(SL2(Z))-modules 


J ,W(E ) w,W (Es) W(Es) ._ W (Es) 
J, ,Eg,t : =Q Eg,t i , Ji Eat — Jk Bst 


keZ keZ 
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We aim to describe the explicit structure of the above modules using Sakai’s 
forms A;, A2, B2, A3, B3, A4, B4, As, Bg. These were some algebraically 
independent W(£g)-invariant holomorphic Jacobi forms first constructed 
by Sakai in the study on E-strings!!. Due to its importance, we briefly 
explain how Sakai constructed these forms. One starts with the Eg Jacobi 
theta function Vp.(7,3) which is the unique W(£s)-invariant holomorphic 
Jacobi form of weight 4 and index 1. Acting the index raising Hecke opera- 
tors T_(t) on Vp,, one obtains W(£g)-invariant holomorphic Jacobi forms 
of weight 4 and arbitrary index t: X;(7,3) = 1+ O(q). Sakai’s forms A; are 
constructed as 


Ag (3) = X;(T,3), j =1,2,3,5; Aang) = Vp, (T, 23). 


To construct B+, one first takes an appropriate modular form gą of weight 
2 on the congruence subgroup To(t) of SL2(Z). Then the trace sum of 
gt(T)U ng (tT, t3) with respect to the cosets of I'9(t)\SL2(Z) defines a W (Ex)- 
invariant holomorphic Jacobi form of weight 6 and index t. That is the 
desired B;. 

Rather surprisingly, Del Zotto, Gu, Huang, Kashani-Poor, Klemm and 
Lockhart !® discovered an exceptional W(Eg)-invariant holomorphic Jacobi 
form of weight 16 and index 5 defined by the polynomial Py¢,5: 


864A? A> + 3825A1 BZ — 770A3B2E6 — 840A2B3 Es + 60A1 By Bg + 2145 E$. 


They checked numerically that Pie, vanishes at the zero points of E4 for 
general lattice variable 3 and then conjectured that the quotient Pig,5/E4 is 
holomorphic. They did not find other similar polynomials, so they further 
conjectured that any Jacobi form expressed as a polynomial in A;, B; and 
Es which vanishes at the zero points of E4 must be divisible by the above 
polynomial. In! we proved their conjectures. Utilizing this distinguished 
Jacobi form Pj¢,5/E4, we are able to give a full description of W (Es)- 
invariant Jacobi forms of arbitrary index in terms of Sakai’s forms. 


Theorem 3.1. 


(1) The quotient Pis, /E4 is a W (Es)-invariant holomorphic Jacobi form 
of weight 12 and index 5. 
(2) For any W(Eg)-invariant Jacobi form P € C[E6, Ai, Bj], if P/E. is 
holomorphic on H x (Eg @ C), then 
i? 


16,5 


€ C[E6, Ai, Bi]. 
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(3) Every W(Es)-invariant weak Jacobi form of index t can be expressed 
uniquely as 
Diy PEPYS 
AN: BE} 
where 
(i) tı is the integer part of t/5; 
(ii) Pa € C[E4, Ee, Ai, Bi]; 
(iii) P; € C[E¢, Ai, Bj] for0 <j < ti; 
(iv) N; is defined as follows 
5to, if t = 6to or 6to + 1, 
5to +1, t= Gig +2, 
5to +2, ift= 6to +3, 
Sto +3, ift= 6to +4 or 6to +5. 


N, 


To determine the explicit generators for each index, we first calculate the 
Fourier expansions of Sakai’s forms A; and Bj up to q?-terms. The Fourier 
expansions involve 268 Weyl orbits of vectors of norm $(v,v) < 54. We then 
express these Weyl orbits as polynomials in the eight fundamental Weyl 
orbits. Using these data and the above theorem, we successfully determine 
all generators of F A for 1 < t < 13 and following generating series 
for the space dimensions. 


Theorem 3.2. Let dy denote the number of generators of weight k of 
co. For 1 < t < 13 the Laurent polynomials 


PY =X dpa" 
keZ 


describing the weights of generators are determined as follows 


Pat, Perta +1, P¥ =a $ +e 4a 4407 +1 


3 


PY =z 164 > ag Zig 10 1 2g Sx Gig tygr 24i, 


PY =2g7 t6 + 2974 4 3e T 4 Oe Sta Sta aH, 


P% = 2075? + 1075 + 240-* + 322718 + 370-4 + 28171? + 290-1 
+ 282— 38 + 262-36 + 230-34 4 22x73? + 18273 + 160-28 + 140-76 


+ 1g Og + 8x7? Ge + 5x716 te Ye 


Qe + 2278 +e Hette. 


162 


The Laurent expansion of the following rational function at x = 0 gives the 
dimension of the space of weak Jacobi forms of arbitrary weight and given 
index t 


Fey = en = dim IPR at 


(1—xt)(1— x8) (1—24*)(1—2°) = 

We further compute the dimension of the space of W (Eg )-invariant holo- 
morphic Jacobi forms of weight 4 and small index and construct explicit 
generators. These so-called holomorphic Jacobi forms of singular (i.e. possi- 
ble minimal positive) weight are usually difficult to determine and construct 
in the theory of modular forms. 


Proposition 3.1. The dimension of the space Jae for t < 11 is for- 
mulated in Table 8.1. 


Table 1. The dimension of ie 


4. Conway Invariant Jacobi Forms on the Leech Lattice 


Let J? and JO, denote the spaces of Conway invariant weak and 
holomorphic Jacobi forms of integral weight and given index t respectively. 
The following is our main theorem?. 


Theorem 4.1. As free modules over M,(SL2(Z)), 


(1) Te is generated by four forms of weights —4, —2, 0, 0. 
(2) I, is generated by four forms of weights 12, 12, 14, 16. 
(3) Shey is ei by ten forms of weights —14, —12, —12 
—8, —6, —2,0. 
(4) Jo, is eee by ten forms of weights 12, 12, 12, 14, 14, 16, 16, 
16, 18, 18. 


; F , —12, —10, 


To prove the above theorem, we first use the differential operators ap- 
proach in” to estimate the minimal weight of weak Jacobi forms of a given 
index. Then we combine the arguments int™!! to construct generators 
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such as Hecke operator. We also construct one of the singular-weight gen- 
erators of ean as the t-th Fourier—Jacobi coefficient of Borcherds’ au- 
tomorphic form ®ı2 for the unimodular lattice H262 (see!®). The main 
difficulty of the proof is to calculate the Fourier expansions of generators, 
because Conway invariant Jacobi forms have unwieldy Fourier expansions 
in 25 variables. To overcome this difficulty, we write the Fourier expansion 
of a Jacobi form in terms of Conway orbits defined as the Coo-invariant 


exponential polynomials 


orb(v) = 5 gored), 


a€Coo /(Coo)» 


where v € A and (Coo), is the stabilizer of Coo with respect to v. The 
Conway orbits orb(v) of type $(v,v) < 16 are available in?°. In order to 
calculate the Fourier expansions of products of Jacobi forms, we have to 
know the decomposition of some products orb(v) orb(w) into linear combi- 
nations of Conway orbits. We determine such non-trivial decompositions 
by comparing the Fourier-Jacobi expansion of ®j2 and the Borcherds de- 
nominator formula for the fake monster Lie algebra (see 1419) 
these arguments together, we prove the theorem. 

The differential operator H also called heat operator used in the proof 
and the construction of generators is defined by the following lemma. 


. Combining 


Lemma 4.1. Given a Conway invariant weak Jacobi form of weight k and 
index t > 1 


Ady E orob: 
n=0 rEA/ Coo 


Then Hp(p) is a Conway invariant weak Jacobi form of weight k +2 and 
index t, where 


noed- D (n- G2) iane: obt) 


n€N rEA/ Coo 


and Eo(T) = 1 — 24 „>; 0(n)q” is the Eisenstein series of weight 2. 


Similar to Sakai’s A; and B; forms for the Eg case, we can also define 
analogously Conway invariant holomorphic Jacobi forms A; and B; for the 
Leech lattice which have weight 12 and 14 respectively. The other type 
of Conway invariant holomorphic Jacobi forms of singular weight 12 come 
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from the famous Borcherds automorphic form ®j given by the Weyl—Kac-— 
Borcherds denominator identity of the fake monster Lie algebra. In terms 


of Conway invariant Jacobi forms, one can express ®19 as follows (see +?:?') 


®12(Z) = A(T) - exp ( = 5 (ATAT aege) 


m=i 


[e.e] 
= X. Disa ™, 
m=0 
where T_(m) are Hecke operators. Then ®12 (7,3) is a Conway invariant 
holomorphic Jacobi forms of singular weight 12 and index m. The above 
operator and basic forms allow us to construct the explicit generators for 
the space of holomorphic forms at index 2 and 3. 


Theorem 4.2. The free C[E4, E¢|-module Jona is generated by Ag, ®12,2, 
Bə and H B2, which have weights 12, 12, 14, 16 respectively. 


Theorem 4.3. The free C[E4, E¢|-module Jo is generated by ten forms 
As, ®12,3, V123, Bs, Via3, HB3, HY14,3, Vie3, H?Bs, H?’Y14,3 
which have weights 12, 12, 12, 14, 14, 16, 16, 16, 18, 18 respectively. 


For the explicit construction of holomorphic forms Y12,3, Y14,3 and Vie 3, 
we refer to”. We also use these holomorphic generators to construct the 
weak generators”. 

As applications, we determine many product decompositions of Conway 
orbits by means of modular linear relations among Conway invariant holo- 
morphic Jacobi forms. These results are formulated in Appendix A of?. 
For example, 


O2 ® Oz = 19656009 © 460002 P 55203 © 4604 P 205 D 206p P Osa, 
O2 ® O3 = 4710402 © 1117803 © 204804 © 27505 Ð 24060 © O7 ® Ope. 


Here the Conway orbits follow the notion of the ALTAS?°. It would be 
very difficult to compute these product decompositions in a brutal way due 
to the huge size of the Conway orbits. We also classify Conway invariant 
holomorphic Jacobi forms of singular weight 12 and index t < 3 with non- 
trivial character and use the Fourier expansions of our Jacobi forms to 
determine all conjugate relations among Conway orbits of type s(v, v) < 16 
modulo 2A and 3A. Besides, we calculate the pullbacks of Conway invariant 
Jacobi forms and Conway orbits along Leech vectors of types 2, 3 and 4, 
which characterize the intersection behaviors of Leech vectors. 
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We know from Borcherds’ thesis?” that A/4A have 31 orbits with respect 
to Cop. We give an explicit description of the representative system of 
minimal length of A/4A in?. Borcherds’ result yields that the rank of Joo 
is 31. It seems very difficult to determine and construct the associated 31 
generators of index 4. 

The famous 2d Monster CFT has one single character x = O,(T)/n74 — 
24 where O;(T) = Oz(T,3 = 0). This can be seen as the index one case 
of the discussion in this section. It is interesting to consider whether there 
exist certain level two generalization of 2d Monster CFT which has 4 char- 
acters, and level three generalization which has 10 characters, and so on. 
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